OFFSET
1,2
PROG
(PARI) {print1(a=1, ", "); for(n=2, 17, 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(sqrtint(p), ", "); b=0, q++)); s=10*s; t=10*t); a=p)}
CROSSREFS
KEYWORD
base,nonn
AUTHOR
Amarnath Murthy, Nov 22 2002
EXTENSIONS
Edited and extended by Klaus Brockhaus, Dec 06 2002
STATUS
approved