OFFSET
1,2
LINKS
Jinyuan Wang, Table of n, a(n) for n = 1..10000
EXAMPLE
15 is a term because the 15th positive nonprime in A018252 is 24, and 1 + 5 = 2 + 4 = 6. - Bernard Schott, Feb 04 2019
MAPLE
A007953 := proc(n) add(d, d=convert(n, base, 10)) ; end proc:
A018252 := proc(n) option remember ; if n = 1 then 1; else for a from procname(n-1)+1 do if not isprime(a) then return a; end if; end do: end if; end proc:
for n from 1 to 900 do if A007953(n) = A007953(A018252(n)) then printf("%d, ", n) ; end if; end do: # R. J. Mathar, Apr 20 2010
PROG
(PARI) k=0; for(n=1, 684, k++; while(isprime(k), k++); if(sumdigits(n)==sumdigits(k), print1(n, ", "))) \\ Jinyuan Wang, Feb 04 2019
CROSSREFS
KEYWORD
nonn,base,less
AUTHOR
Juri-Stepan Gerasimov, Apr 14 2010
EXTENSIONS
More terms from R. J. Mathar, Apr 20 2010
STATUS
approved