%I#10 2018年6月15日09:26:45
%S 1,2,3,4,5,6,7,8,4,4,4],9,4,10,4,11,12,13,14,15,16,17,18,19,20,21,22,
%电话23,4,4,24,25,4,26,4,27,28,4,29,30,4,31,4,32,33,34,35,4,36,37,38,39,
%U 40,41,4,42,4,4,43,44,45,46,47,48,49,50,51,4,52,53,54,55,56,57,58,59,60,61,62,63,64,4,65,66,67,4,68,4,69,70,71,72,73,74,4
%N限制性生长序列转换A305900(A048673(N))。
%C对于所有i、j:
%C a(i)=a(j)=>A278224(i)=A278222(j)。
%C a(i)=a(j)=>A286583(i)=A286582(j)。
%C a(i)=a(j)=>A292251(i)=A29225一(j)。
%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,长度(invesc),如果(mapisdefined(om,invec[i]),my(pp=mapget(om,invec[i];
%o A003961(n)=我的(f=系数(n));对于(i=1,#f~,f[i,1]=下一素数(f[i、1]+1));factorback(f);\\来自A003961
%o A048673(n)=(A003961(n)+1)/2;
%o A305900(n)=如果(n<=5,n,if(i素数(n),5,3+n-素数(n)));
%o v305902=rgs_transform(向量(up_to,n,A305900(A048673(n)));
%o A305902(n)=v305902【n】;
%Y参考A048673,A305900。
%Y另请参阅A305901。
%K nonn公司
%O 1,2号机组
%A _Antti Karttunen,2018年6月14日
|