%I#19 2016年4月1日10:39:45

%S 4,83634845445464642522138312123,0,442323231213686，

%电话：531876781354494427224642171,0525,88575464522122562426，

%电话：36963252464,0178379973871291080192211241145325696332010293639,06568588

%N a（N）=最小回文k>N，使得k/N是正方形；如果不存在解，则a（n）＝0。

%C如果n是10的倍数，那么a（n）=0，因为没有回文以0结尾。

%C最多200个术语目前只有3个未知，a（125）>5.2*10^28，a（177）>3.5*10^27和a（185）>4.5*10^27。有关已知值的表格，请参见链接。-_Giovanni Resta_，2015年8月5日

%C如果a（125）>0，则a（125。证明：设m^2=a（125）/125。那么m是奇数，否则125*m^2是10的倍数，不是回文。因为m是奇数，所以m^2==1模8，因此125*m^2==125模1000_Chai Wah Wu_，2016年3月31日

%H Giovanni Resta，已知值小于等于a（200）</a>

%e a（3）=363，因为363/3=11^2。363*3=1089，这也是一个正方形。

%e a（15）=53187678135，因为531876781 35/15=59547^2和53187678 135*15=893205^2。

%p ispali：=proc（n）局部L；五十： =换算（n，基数，10）；ListTools:反向（L）=L结束进程：

%p f：=proc（n）局部m；

%p如果n mod 10=0，则返回0 fi；

%对于从2到10^6的m，p如果是ispali（m^2*n），则返回m^2*n fiod：

%p-1#信号超时

%p端程序：

%p序列（f（n），n=1..50）；#_罗伯特·伊斯雷尔（Robert Israel），2015年8月21日

%t palQ[n_]：=块[{d=整数位数@n}，d==反向@d]；a[n_]：=如果[Mod[n，10]==0，0，块[{q=2}，而[！palQ[q^2*n]，q++]；q^2*n]]；阵列[a，42]（*Giovanni Resta_，2015年8月18日*）

%o（Python）

%o对于范围（1150）内的k：

%哦。。。。c=2

%哦。。。。而c<10**8：

%o。。。。。。。。kk=k*c**2

%o。。。。。。。。如果kk==int（str（kk）[：：-1]）

%o。。。。。。。。。。。。打印（k，kk）

%o。。。。。。。。。。。。c=10**9

%o。。。。。。。。c=c+1

%Y参考A002113、A000290、A023108、A061563。

%K nonn，基础

%O 1,1号机组

%A _邮政信箱_，2015年7月30日

%E 2015年8月5日，_Giovanni Resta中缺少a（13）