%I#18 2018年3月13日22:11:12
%S 0,1,1,2,1,4,1,6,2,4,1,60,1,4,42,1,60,4,4,14620,2,4,6,60,1,
%电话:1000,1546,4,421780,1,4,44620,11000,1,60,60,4,11021020,2,60,4,
%U 60,14620,44620,4,416897000,14,6012558,41000,160,41000:175162780,14,60,41000,11021020,42,4,1
%N a(1)=0;对于n>1,a(n)=Product{d|n,1<d<n}素数(A101296(d)-1)。
%C a(n)=由素数与n的真除数的素数签名索引而成的乘积。
%这个序列的限制增长序列变换是A101296,因为从n的真除数的素数签名集总是可以确定n本身的素数,反之亦然,从n的素数签字,我们可以形成其所有真除数素数签名的集合。
%C对于所有i,j:a(i)=a(j)<=>A101296(i)=A101296(j)。
%H Antti Karttunen,n的表,n的a(n)=1..8192</a>
%H<a href=“/index/Eu#epf”>根据n的因式分解中的指数计算序列的索引项</a>
%e对于n=12,其真除数>1为2,3,4,6,其从A101296起的素数签名秩为:2,2,3和4。我们从每一个中减去一个,得到乘积素数(1)*素数(一)*素值(二)*素元(三)=2*2*3*5=60,即a(12)的值。
%t块[{nn=83,s},s=Map[#1->#2&@@#&,Transpose@{Values@#,Keys@#}]&@PositionIndex@Table[Times@@MapIndexed[Prime[First@#2]^#1&,Sort[FactorInteger[n][[All,-1]],Greater]]-Boole[n==1],{n,nn}];表[If[n==1,0,Times@@Map[Prime[FirstPosition[Keys@s,#][[1]-1]&,Most@Rest@Divisors@n]],{n,nn}]](*_Michael De Vlieger_,2018年3月13日*)
%o(PARI)
%o up_to=8192;
%o rgs_transform(invec)={my(om=Map(),outvec=vector(length(invec)),u=1);对于(i=1,长度(invesc),如果(mapisdefined(om,invec[i]),my(pp=mapget(om,invec[i];
%o写入to_b文件(start_offset,vec,bfilename)={对于(n=1,长度(vec),写入(bfilename,(n+start_offset)-1,“”,vec[n]);}
%o A046523(n)={my(f=vecsort(因子(n)[,2],4),p);prod(i=1,#f,(p=nextprime(p+1))^f[i]);};\\来自A046523
%o v101296=rgs_transform(向量(up_to,n,A046523(n)));
%o A101296(n)=v101296【n】;
%o A300716(n)={my(m=1);如果(1==n,0,fordiv(n,d,if(d>1)&(d<n),m*=素数(A101296(d)-1));(m));};
%o表示(n=1,up_to,写入(“b300716.txt”,n,“”,A300716(n));
%Y参见A046523、A101296、A297174、A300250。
%Y另请参阅A292258、A294897、A300827。
%K nonn公司
%O 1,4型
%A _Antti Karttunen,2018年3月13日
|