%I#8 2017年10月18日00:56:12

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

%电话20,2,21,22,23,24,25,2,26,27,28,2,29,2,30,31,32,2,33,34,35,12,36,2,37，

%U 38,39,40,41,2,42,2,43,44,45,46,47,2,48,49,50,2,51,252,53,54,55,56,2,57,58,59,2,60,61,62,63,64,2,65,66,67,68,69,70,71,2,72

%A293225的限制生长序列变换，这是一种组合两个乘积的滤波器，另一个乘积是由N的真除数的三元展开式中的1位数字（A293221）和2位数字（A2 93222）形成的。

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

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

%o（PARI）

%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写入to_b文件（start_offset，vec，bfilename）={对于（n=1，长度（vec），写入（bfilename，（n+start_offset）-1，“”，vec[n]）；}

%o A019565（n）={my（j，v）；因子回复（Mat（向量（如果（n，#n=向量提取（二进制（n），“-1..1”）），j，[素数（j），n[j]]）~））}；\\此函数来自_M.F.Hasler_

%o A289813（n）={my（d=数字（n，3））；从数字（向量（#d，i，if（d[i]==1，1，0）），2）；}；

%o A289814（n）={my（d=数字（n，3））；从数字（向量（#d，i，如果（d[i]==2，1，0）），2）；}；

%o A293221（n）={my（m=1）；对于div（n，d，如果（d<n，m*=A019565（A289813（d）））；m；}；

%o A293222（n）={my（m=1）；对于div（n，d，如果（d<n，m*=A019565（A289814（d）））；m；}；

%o Anot_submited（n）=（1/2）*（2+（（A293222（n）+A293221（n））^2）-A293222（n）-3*A293221n））；\\等式等级等于A293225。

%o写入to_b文件（1，rgs_transform（矢量（19683，n，Anot_submited（n））），“b293226.txt”）；

%Y参见A001065、A293221、A29322、A293233、A293244、A293225。

%Y另请参阅A290093、A290094、A293215、A29321、A293232。

%K nonn公司

%O 1,2号机组

%2017年10月3日，安蒂·卡图内