%I#6 2019年9月5日20:49:15
%S 1,2,2,2,2,3,2,2,2,4,2,5,2,4,6,2,7,2,8,2,42,2,9,2,2,5,2,6,4,
%T 2,4,2,4,2,2,7,2,5,10,4,2,5,2,4,12,4,4,7,2,1,1,2,4,2,7,2,
%U 4,6,4,2,13,2,4,2,5,2,7,2,42,4,4,2,2,5,14,15,5,2,4,16,9,2,46,8,2,72,4,6
%N词汇学上最早的序列,对于所有i,j,a(i)=a(j)=>f(i)=f(j),其中f。
%C函数f的限制增长序列变换,当n是酉乘完美数(A327158)之一时定义为f(n)=-A323166(n)=-n,否则定义为f
%C对于所有i、j:
%C A305800(i)=A305800。
%H Antti Karttunen,n的表,n的a(n)=1..87360</a>
%o(PARI)
%o上限=87360;
%o rgs_transform(invec)={my(om=Map(),outvec=vector(length(invec)),u=1);对于(i=1,length(invec),if(mapisdefined(om,invec[i]),my(pp=mapget(om,invec[i]));outvec[i]=outvec[pp],mapput(om,invec[i],i);outvec[i]=u;u++);outvec;};
%o A034448(n)={my(f=因子(n));prod(k=1,#f~,1+(f[k,1]^f[k,2]));};\\A034448中的后代码
%o A323166(n)=gcd(n,A034448(n));
%o Aux327163(n)={my(u=A323166(n));u*((-1)^(u==n));
%o v327163=rgs_transform(向量(up_to,n,Aux327163(n)));
%o A327163(n)=v327163【n】;
%Y参见A034448、A305800、A323166、A327158和A327164。
%Y另请参阅A326193。
%K nonn公司
%O 1,2号机组
%A _Antti Karttunen_,2019年8月28日
|