A078815 a(n) = A078215(n+1)/A078215(n).

1, 25, 121, 16, 121, 1089, 121, 1089, 1089, 1156, 1089, 1089, 1089, 1089, 1024, 1089, 10404, 1089, 1024, 1089, 1024, 10404, 10404, 10404, 10404, 10404, 10404, 10404, 1024, 10404, 10404, 1024, 10404, 1024, 1024, 10404, 1024, 1024, 1024, 1024, 1024

Offset

1,2

Links

Table of n, a(n) for n=1..41.

Prog

(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)}

Crossrefs

Cf. A078215, A078216.

Sequence in context: A360511 A124953 A126412 * A270368 A206472 A036057

Adjacent sequences: A078812 A078813 A078814 * A078816 A078817 A078818

Keyword

base,nonn

Author

Klaus Brockhaus, Dec 06 2002

Status

approved