%I#6 2017年7月22日10:13:16

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

%电话：22,23,24,25,26,27,28,29,23,30,31,32,33,34,20,35,36,37,30,39,40,41，

%U 42,19,43,44,45,37,46,47,23,48,49,50,19,51,52,53,54,34,55,56,57,58,59,60,62,30,63,64,65,66,42,67,68,34,69,30,70,71,68,72,61

%N限制性生长序列转换A289625（A003961（N））。

%C对于所有i，j:a（i）=a（j）=>A003972（i）=A003971（j）。

%H Antti Karttunen，<a href=“/A290082/b290082.txt”>n的表，a（n）表示n=1.16384</a>

%o（PARI）

%o rgs_transform（invec）={my（occurrences=Map（），outvec=vector（length（invec）），u=1）；对于（i=1，length，invec，if（mapisdefined（occurements，invec[i]），my（pp=mapget（occursions，invec[i]；

%o写入to_b文件（start_offset，vec，bfilename）={对于（n=1，长度（vec），写入（bfilename，（n+start_offset）-1，“”，vec[n]）；}

%o A003961（n）={my（f=因子（n））；对于（i=1，#f~，f[i，1]=下一素数（f[i、1]+1））；因子回退（f）；}；\\此函数来自_Michel Marcus_

%o A289625（n）={my（m=1，p=2，v=znstar（n）[2]）；对于（i=1，长度（v），m*=p^v[i]；p=nextprime（p+1））；（m）；}；

%o写入to_b文件（1，rgs_transform（矢量（16384，n，A289625（A003961（n））），“b290082_upto16384.txt”）；

%Y参见A003961、A003972、A289625、A290076。

%K nonn公司

%O 1,2号机组

%2017年7月20日，安蒂·卡图内