OFFSET
1,1
EXAMPLE
An eight-digit term is 11224121, a ten-digit term is 1111111843.
11353 belongs to the sequence because 1^2+1^2+3^2+5^2+3^2=45=1*1*3*5*3
PROG
(PARI) forprime(p=2, 6e6, n=p; sd=0; pd=1; while(n>0, d=divrem(n, 10); n=d[1]; sd=sd+d[2]*d[2]; pd=pd*d[2]); if(sd==pd, print1(p, ", ")))
CROSSREFS
KEYWORD
easy,nonn,base
AUTHOR
Felice Russo, Feb 06 2002
EXTENSIONS
Edited and extended by Klaus Brockhaus Feb 11 2002
STATUS
approved