|
%I #7 Mar 30 2012 18:35:54
%S 3,6,7,11,20,29,30,38,56,63,89,114,124,131,139,140,176,184,192,195,
%T 198,208,214,215,223,243,300,301,331,337,345,359,372,377,388,393,407,
%U 411,451,455,459,461,476,482,483,487,506,508,516,523,536
%N Numbers n such that the sum of the square of digits of n^n is a prime.
%e 6 is in the sequence because 6^6 = 46656 and 4^2 + 6^2 + 6^2 + 5^2 + 6^2 =
%e 149 is a prime.
%p with(numtheory):for n from 1 to 540 do: l:=length(n^n) : n0:=n^n: s:=0:for
%p m from 1 to l do: q:=n0:u:=irem(q,10):v:=iquo(q,10): n0:=v :s:=s+u^2: od:if
%p type(s,prime)=true then printf(`%d, `, n):else fi:od:
%K nonn,base
%O 1,1
%A _Michel Lagneau_, Nov 17 2010
|
