A058196 Sum of digits of first 10^n primes.
2, 57, 1111, 17024, 224155, 2674957, 31352962, 363506969, 4127383480, 45630284509, 502160313778
Offset: 0
Programs
-
Mathematica
NextPrime[ n_Integer ] := Module[ {k}, k = n + 1; While[ ! PrimeQ[ k ], k++ ]; k ]; c = d = p = q = 0; Do[ While[ d++; d <= 10^n, q = NextPrime[ q ]; p = p + Apply[ Plus, RealDigits[ q ][ [ 1 ] ] ]; If[ PrimeQ[ p ], c++ ] ]; d--; Print[ p ], {n, 0, 10}]
-
PARI
a(n) = {my(k = 0, s = 0, pow = 10^n); forprime(p = 2, , k++; s += sumdigits(p); if(k == pow, break)); s;} \\ Amiram Eldar, Jun 03 2024
Extensions
a(9)-a(10) from Amiram Eldar, Jun 03 2024