%I#11 2018年9月26日16:24:28

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

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

%U 45,46,47,38,48,49,50,51,52,53,54,42,55,56,57,58,59,11,60,61,62,20,63,33,64,65,66,67,68,42,69,70,71,38,72,73,74,75,76,38,77,78,79,80,81,82,83,11,84,85,86,38

%N过滤序列，对于素数p，它在二进制表示法[A286622（p）]中为每个由1行长度组成的不同多集记录一个不同的值，并为所有其他数字指定一个唯一的数字。

%C当n是素数时，函数f的限制增长序列变换定义为f（n）=A278222（n），否则为-n。

%C在初始项3之后，Fermat素数（A019434）给出了这个序列中5的位置，而Mersenne素数（P000668）则分别被分配给它们自己的单粒子等价类。

%C对于所有i、j：

%C a（i）=a（j）=>A305900（i）=A305900，

%C a（i）=a（j）=>A286622（i）=A286621（j），

%C a（i）=a（j）=>A305795（i）=A305796（j）。

%H Antti Karttunen，n的表，n=1..100000的a（n）</a>

%o（PARI）

%o up_to=100000；

%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 A005940（n）={my（p=2，t=1）；n--；直到（！n\=2，if（（n%2），（t*=p），p=下一素数（p+1））；t}；

%o A046523（n）={my（f=vecsort（factor（n）[，2]，4），p）；prod（i=1，#f，（p=nextprime（p+1））^f[i]）；}；\\来自A046523

%o A278222（n）=A046523（A005940（1+n））；

%o A319705aux（n）=如果（i素数（n），A278222（n）和-n）；

%o v319705=rgs_transform（向量（up_to，n，A319705aux（n）））；

%o A319705（n）=v319705【n】；

%Y参见A000668、A019434、A278222、A286163、A286622、A305795、A305900。

%Y另请参阅A319704、A319706。

%K nonn公司

%O 1,2号机组

%A _Antti Karttunen，2018年9月26日