%I#19 2019年11月2日06:53:47
%S 7,17,15,19,5,43,19,517,7167,1,11,3,19,1,67,1,17,17,7,5211,1,7,
%电话:11,13,3139,1,31,9,1,5109,1,1,3,85,3,61,11,11,11,7,1211,11,5,3,
%U 31,5,31,1,13,1353,1,1,9,31,3,71,1,5,3,19,1317,1,5,1,1,31167,7,5
%N a(N)是最小的m,因此2*N*m是一个不明显的值;即2*n*a(n)在A005277中。
%H Amiram Eldar,n的表,a(n)表示n=1..10000</a>
%e n=5:invphi(10k)中的项数为2,5,2,9,0,9,。。。对于k=1,2,3,。。。;a(5)=5,因为0出现在第5位。
%p带有(numtheory);a:=proc(n)局部m;对于1中的m,如果(invphi(2*n*m)=[]),则返回m end end
%t invphi[n_,plist_]:=模块[{i,p,e,pe,val},如果[plist=={},返回[If[n==1,{1},{}]];val={};p=最后[plist];对于[e=0;pe=1,e==0||Mod[n,(p-1)pe/p]==0,e++;pe*=p,val=Join[val,pe*invphi[If[e==0,n,n*p/pe/(p-1)],Drop[plist,-1]]];排序[val]];invphi[n_]:=invphi[n,选择[1+除数[n],素数Q]];a[n_]:=对于[m=1,True,m++,如果[invphi[2n*m]=={},返回[m]](*invphi[n,plist]是x的列表,其中phi(x)=n和plist中x的所有素数除数。*)
%Y参考A0000010、A005277、A0002202、A071616。
%K nonn公司
%O 1,1号机组
%A _Labos Elemer,2002年5月27日
%E由_Robert G.Wilson v_编辑和扩展,2002年5月28日
%E由_Dean Hickerson编辑和扩展,2002年6月4日