%I#20 2017年12月17日08:02:57

%S 9、12、17、19、24、25、28、33、40、49、51、52、57、60、64、67、72、73、79、81、84、88、89，

%电话96,97,99103105108112115116121124129134136144145148156，

%U 161163168169172177180184192193199电话：

%基数，其中ab和ba是不同的正方形，对于某些a和b。

%C自2016年3月14日以色列罗巴特起：（开始）

%C不允许使用前导0。

%C猜想：除1以外的所有奇数平方（A016754）都是序列项。（结束）

%CN=（2n+1）^2，a=N^2，b=4n^2+2n+1表明（2n+1）^2是一个项，所以这个序列是无限的_Michael R Peake，2017年3月21日

%H Robert Israel，n的表，n=1..10000的a（n）</a>

%e 17_9和71_9是正方形。14_12和41_12是正方形。

%p滤波器：=proc（n）局部x，a，b，R；

%x的p从ceil（sqrt（n））到n-1 do

%pa:=x^2模块；

%p如果a=0，则下一个fi；

%p b：=（x^2-a）/n；

%p如果赋值（R[b，a]），则返回真fi；

%p R[a，b]：=1；

%p od；

%p假

%p端程序：

%p选择（过滤器，[1..1000]）；#_罗伯特·伊斯雷尔（Robert Israel），2016年3月14日

%o（MATLAB）

%o匹配=零（1100）；

%o对于N=2:200，十=零（1，N-1）；单位=零（1，N-1）；对于a=N-1:-1:sqrt（N），c=a^2；十（a）=楼层（c/N）；单位（a）=rem（c，N）；结束；对于a=N-1:-1:sqrt（N），h=find（（单位==十（a））&（[1:N-1]~=a））；如果长度（h），匹配=任意（单位（a）=十（h））；如果匹配，Sol（N）=Sol（N）+1；结束；结束；结束；结束；

%o查找（匹配>0）

%Y参考A016754。

%K基础，简单，无

%O 1,1号机组

%A _迈克尔·R·派克，2010年3月21日

%E MATLAB程序由_Robert Israel_修订，2016年3月14日