A123157 Sum of digits of the squares of prime numbers.

4, 9, 7, 13, 4, 16, 19, 10, 16, 13, 16, 19, 16, 22, 13, 19, 16, 13, 25, 10, 19, 13, 31, 19, 22, 4, 16, 19, 19, 25, 19, 16, 31, 16, 7, 13, 25, 28, 34, 31, 10, 19, 22, 25, 28, 19, 16, 31, 22, 16, 28, 16, 22, 10, 25, 31, 19, 19, 31, 31, 25, 34, 28, 25, 40, 22, 22, 25, 16, 13, 22, 28

Offset: 1

Author

Herman Jamke (hermanjamke(AT)fastmail.fm), Oct 01 2006

Keywords

• nonn
• base

Links

• Albert Frank, International Contest Of Logical Sequences, 2002 - 2003. Item 11
• Albert Frank, Solutions of International Contest Of Logical Sequences, 2002 - 2003.

Programs

Maple

read("transforms"):
A123157 := proc(n)
    digsum(ithprime(n)^2) ;
end proc: # R. J. Mathar, Jul 05 2013

Mathematica

Total[IntegerDigits[#]]&/@(Prime[Range[80]]^2) (* Harvey P. Dale, Nov 04 2012 *)

PARI

for(n=1,100,dp2=eval(Vec(Str(prime(n)^2)));print1(sum(i=1,length(dp2),dp2[i]),","))