%I#51 2024年2月12日13:28:56
%S 0,1,2521618493648131190426666895260600944976578435364121,
%电话:758578289296485619490565545012396225109386936710041,
%电话:9352258807734007959040299927611577358868780506251348595071513294176115300880662371655692274550391912122916411944673157896228824166263789956099627581699900092104464129
%N数N,使三角形(N)+三角形(2*N)为正方形。
%前10个术语中的C8的格式为x^y。两个例外是a(7)=311904=2^5*3^3*19^2和a(10)=449765784=2^3*3^5*13^2*37^2。
%C对应的平方由A075873(2*n-1)^2给出。例如,三角形(a(10))+三角形(2*a(10))=711142146^2=A075873(19)^2。
%C A147875中正方形的位置,相当于解丢番图方程n*(5*n+3)=2*s^2。-_R.J.Mathar,2013年4月19日
%H Max Alekseyev,n的表,n=1..1000时的a(n)</a>
%H<a href=“/index/Rec#order_07”>具有常系数的线性重复出现的索引条目,签名(1,01442,-1442,0,-1,1)。
%F a(n)=(A228209(2*n-1)-3)/10.-_Max Alekseyev_,2013年9月4日
%传真:x^2*(x+1)*(x^4+23*x^3+168*x^2+23*x1)/(x^6-1442*x^3+1)/(1-x)_Max Alekseyev_,2013年9月4日
%t线性递归[{1,01442,-1442,0,-1,1},{0,1,25216184936481311904},30](*哈维·P·戴尔,2015年1月23日*)
%o(Python)
%o导入数学
%o对于i在范围内(1L<<30):
%o s=i*(i+1)/2+i*(2*i+1)
%o t=int(数学.sqrt(s))
%o如果s==t*t:print i,
%Y参见A000217、A075873。
%Y参考A220186(数字n使得三角形(2*n)-三角形(n)是正方形)。
%K nonn公司
%氧1,3
%A _Alex Ratushnyak,2013年4月18日
%E 2013年9月4日最大Alekseyev之后的第a(11)条
|