login
a(n) = A078215(n+1)/A078215(n).
2

%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