|
%I #3 Mar 30 2012 17:27:40
%S 1,25,121,16,121,1089,121,1089,1089,1156,1089,1089,1089,1089,1024,
%T 1089,10404,1089,1024,1089,1024,10404,10404,10404,10404,10404,10404,
%U 10404,1024,10404,10404,1024,10404,1024,1024,10404,1024,1024,1024,1024,1024
%N a(n) = A078215(n+1)/A078215(n).
%o (PARI) {print1(a=1,","); for(n=2,41,k=floor(log(a)/log(10))+1; s=n; t=s+1; m=floor(log(s)/log(10))+1; d=k-m; s=s*10^d; t=t*10^d; b=1; while(b>0,q=floor(s/a); while(b>0&&(p=q*a)<t,if(p>=s&&issquare(p),print1(q,","); b=0,q++)); s=10*s; t=10*t); a=p)}
%Y Cf. A078215, A078216.
%K base,nonn
%O 1,2
%A _Klaus Brockhaus_, Dec 06 2002
|
