登录
词法上最早的无限序列,使得a(i)=a(j)=>f(i)=f(j)对于所有i,j>=1,其中f(1)=0,对于n>1,f(n)=A334870型(n) ●●●●。
4

%I#7 2020年6月9日22:11:27

%S 1,2,3,3,4,5,6,7,5,8,9,10,11,12,13,7,14,15,16,17,18,19,20,21,8,22,23,

%第24、25、26、27、28、29、30、31、41、32、33、34、35、36、37、38、39、40、41、42、43、12、44页,

%U 45、46、47、48、49、50、51、52、53、54、55、56、57、10、58、59、60、61、62、63、64、65、66、67、68、69、70、71、72、73、15、74、75、76、77、78、79、80、81、82、83、84、85、86、87、88、89、90、91、13、92

%N词法上最早的无限序列,使得a(i)=a(j)=>f(i)=f(j)对于所有i,j>=1,其中f(1)=0,对于N>1,f(N)=A334870(N)。

%H Antti Karttunen,n的表格,n=1..65537的a(n)</a>

%o(PARI)

%o up_to=65537;

%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 A334870(n)=如果(issquare(n),sqrtint(n)),my(c=芯(n)、m=n);对于素数(p=2,如果(!(c%p),m/=p;断裂,m*=p));(m) );

%o A334869aux(n)=如果(1==n,0,A334870(n));

%o v334869=rgs_transform(向量(up_to,n,A334869aux(n)));

%o A334869(n)=v334869【n】;

%Y参见A334868、A334870和A334871。

%K nonn公司

%O 1,2号机组

%2020年6月8日,安蒂·卡图宁