%I#29 2016年12月14日11:06:24
%S 2,3,6,7,8,11,12,15,18,19,22,23,24,25,26,28,29,31,32,33,41,43,
%第44、45、46、47、48、53、54、60、61、67、70、72、74、76、79、82、84、87、90、92、93、96页,
%电话:105106107108111112114117122128133139141148159
%N表示k,使A276516(k)=0。
%C这与A001422不同,第一个差异:a(14)=25,A001421(14)=27。
%C推测:对于k>7169,这个序列中没有更多的项(测试k<10000000)。
%H Vaclav Kotesovec,<a href=“/A276517/b276517.txt”>n表,n=1..173的a(n)</a>
%e3在序列中,因为A276516(3)=0
%e4不在序列中,因为A276516(4)=-1
%e 4222在序列中,因为A276516(4222)=0
%e 7169在序列中,因为A276516(7169)=0
%t nn=100;A276516=静止[系数列表[系列[乘积[(1-x^(k^2))),{k,nn}],{x,0,nn^2}],x]];选择[Range[nn^2],A276516[[#]]==0&]
%t nmax=10000;nn=楼层[Sqrt[nmax]]+1;poly=常量数组[0,nn^2+1];聚[1]]=1;聚[[2]]=-1;poly[[3]]=0;Do[Do[poly[[j+1]]-=多边形[[j-k^2+1]],{j,nn^2,k^2,-1}];,{k,2,nn}];A276516=取[poly,{2,nmax+1}];选择[Range[nmax],A276516[[#]]==0&]
%Y参见A001422、A001661、A276516、A279486和A279487。
%K nonn公司
%O 1,1号机组
%A _ Vaclav Kotesovec_,2016年12月12日
|