(PARI)单(n)=n=eval(Vec(Str(n)));对于(i=2,#n,如果(n[i]<n[i-1],返回(0));1
对于(n=1,1e5,if(单声道(n^2),打印1(n“,”))\\查尔斯·格里特豪斯四世2011年8月22日
(Python)
从itertools导入组合with_replacement
从gmpy2导入issquare,isqrt
A028819号_list=[0]+[int(isqrt(n))for n in(int(''.join(i))for l in range(1,11)for i in combinations_with_replacement('123456789',l))if is_square(n)]#柴华武2015年12月7日
|