OFFSET
1,1
LINKS
Eric M. Schmidt, Table of n, a(n) for n = 1..1000
EXAMPLE
a(6) = 19, sum of digits is 10; a(7) = 29, sum of digits is 11 and 11 > 10.
MATHEMATICA
t = {s = 2}; Do[If[(y = Total[IntegerDigits[x = Prime[n]]]) > s, AppendTo[t, x]; s = y], {n, 2, 750000}]; t (* Jayanta Basu, Aug 09 2013 *)
PROG
(Sage)
def A061248(nterms, b=10) :
res = []; n_list = [2]; n = 2; dsum = 0
while len(res) < nterms :
while not (sum(n_list) >= dsum and n.is_prime()) :
i = next((j for j in range(len(n_list)) if n_list[j] < b-1), len(n_list))
if i == len(n_list) : n_list.append(0)
n_list[i] += 1
r = dsum - sum(n_list[i:])
for j in range(i) :
n_list[j] = min(r, b-1)
r -= n_list[j]
n = sum(n_list[i]*b^i for i in range(len(n_list)))
res.append(n); dsum = sum(n_list)+1
return res
# Eric M. Schmidt, Oct 08 2013
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Amarnath Murthy, Apr 23 2001
EXTENSIONS
More terms from Patrick De Geest, Jun 05 2001
STATUS
approved