%I#28 2021年5月10日02:07:37
%S 3,7,11,13,31,43,61,73127151157211241307331421463521547,
%电话:601683757109312123148317232551273128012971330735413907,
%电话:44234561483151135419570160076163648180118191909199011031131
%N整数k>2、|m|>1和k!=的素数Phi(k,m)2^j(j>1)
%C Phi(k,m)表示分圆多项式数分圆(k,m)。
%C这些是A206944的主要条款。
%C A059055是该序列的子序列。
%Mathematica程序可以生成任意边界maxdata的序列,而无需用户选择参数。A206864中解释了该程序的参数确定部分。
%e仅采用A206944中的主要术语:
%e A206944(1)=3是素数,所以a(1)=3。。。
%t phiinv[n_,pl_]:=模块[{i,p,e,pe,val},如果[pl=={},返回[If[n==1,{1},{}]];val={};p=最后[pl];对于[e=0;pe=1,e==0||Mod[n,(p-1)pe/p]==0,e++;pe*=p,val=Join[val,pe*phiinv[If[e==0,n,n*p/pe/(p-1)],Drop[pl,-1]]];排序[val]];phiinv[n_]:=phiinv[n,选择[1+除数[n],素数Q]];最大数据=12000;最大值=天花板[(1+Sqrt[1+4*(maxdata-1)])/4]*2;eb=2*楼层[(Log[2,maxdata])/2+0.5];而[eg=phiinv[eb];lu=长度[eg];lu==0,eb=eb+2];t=选择[Range[eg[[Length[eg]]],((EulerPhi[#]<=eb)&&((!IntegerQ[Log[2,#]])||(#<=2)))&];t=排序依据[t,分圆[#,2]&];a={};Do[i=2;While[i++;cc=分圆[t[i]],m];cc<=最大数据,如果[PrimeQ[cc],a=Append[a,cc]]],{m,2,max}];联盟[a]
%Y参见A059055、A206944、A20694、A206864、A194712。
%K nonn公司
%O 1,1号机组
%A _雷州_,2012年2月13日
|