登录
A107067号
系数为{0,1}并除以x^n-1的多项式的数目。
1, 2, 2, 4, 2, 6, 2, 8, 4, 6, 2, 17, 2, 6, 6, 16, 2, 18, 2, 17, 6, 6, 2, 48, 4, 6, 8, 17, 2, 36, 2, 32, 6, 6, 6, 69, 2, 6, 6, 47, 2, 36, 2, 17, 17, 6, 2, 136, 4, 18, 6, 17, 2, 54, 6, 47, 6, 6, 2, 176, 2, 6, 17, 64, 6, 36, 2, 17, 6, 36, 2, 257, 2, 6, 18, 17, 6, 36, 2, 131, 16, 6, 2, 177, 6, 6, 6, 47, 2, 183, 6, 17, 6, 6, 6, 389, 2, 18, 17, 70, 2, 36, 2, 47, 35
抵消
1,2
评论
发件人罗伯特·伊斯雷尔2017年5月22日:(开始)
这些多项式中的每一个都是k>1除以n的不同分圆多项式C_k(x)的乘积。
a(n)<=2^(A000005号(n) -1)。
如果n是素数,则a(n)=2。(结束)
链接
罗伯特·伊斯雷尔,n=1..719时的n,a(n)表(前359个术语来自Antti Karttunen)
MAPLE公司
f: =程序(n)局部t,C,x,S;
C: =映射(m->numtheory:-分圆(m,x),numtheori:-除数(n)减去{1});
t: =0:
S: =组合:-子集(C);
虽然不是S[完成]do
如果{coeffs(expand(convert(S[nextvalue](),`*`)),x)}={1},则
t: =t+1;
fi(菲涅耳)
od;
结束进程:
地图(f,[1..100]美元)#罗伯特·伊斯雷尔2017年5月22日
数学
a[n_]:=模块[{c,s},
c=分圆[#,x]和/@休息@除数[n] ;
s=系数列表[#,x]和/@(Times@@@子集[c]);
选择[s,AllTrue[#,#==0||#==1&]&]//长度];
表[a[n],{n,1,105}](*Jean-François Alcover公司2023年10月24日*)
黄体脂酮素
(PARI)对于(n=1100,m=0;p=x^n-1;nE=numdiv(n);P=系数(P);E=P[,2];P=P[,1];对于vec(v=矢量(nE,i,[0,E[i]]),divp=prod(k=1,nE,P[k]^v[k]);m++;对于(j=0,极度(divp),divpcof=polceoff(divp,j);if(divpcof<0||divpcoff>1,m-;break));print1(m,“,”)\\Herman Jamke(hermanjamke(AT)fastmail.fm),2006年11月15日
关键词
非n
作者
拉尔夫·斯蒂芬,根据的建议马克斯·阿列克塞耶夫2005年6月11日
扩展
来自Herman Jamke(hermanjamke(AT)fastmail.fm)的更多条款,2006年11月15日
数据段进一步扩展,并使用Jamke的PARI程序计算b文件安蒂·卡图恩2017年5月22日
状态
经核准的