%I#21 2023年11月5日11:35:46

%S 1,1,2,2,3,2,5,4,6,3,7,4,11,5,6,8,13,6,17,6,10,7,19,8,15,11,18,10,23，

%电话：6,29,16,14,13,15,12,31,17,22,12,37,10,41,14,18,19,43,16,35,15,26,22，

%U 47,18,21,20,34,23,53,12,59,29,30,32,33,14,61,26,38,15,67,24,71,31,30,34

%N素数级联：用a（素数（m）^k）=素数（m-1）*素数（米）^（k-1）进行乘法。

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

%F a（A000040（n））=A000040；

%F a（A000079（n））=A000079；

%F a（A002110（n））=A002110。

%F a（n）=A003557（n）*A064989（A007947（n_Antti Karttunen_，2017年12月31日

%F和{k=1..n}a（k）~c*n^2/2，其中c=Product_{p素数}（1-（p-q（p））/p^2）=0.52621951…，其中q（2）=1，且q（p

%e a（63）=a（3^2*7^1）=a（3|2）*a（7|1）=（2*3^1）*（5*7^0）=30。

%t f[p_，e_]：=如果[p==2，1，NextPrime[p，-1]]*p^（e-1）；a[1]=1；a[n_]：=倍@@f@@FactorInteger[n]；阵列[a，100]（*_Amiram Eldar_，2023年11月2日*）

%o（PARI）

%o A003557（n）={my（f=系数（n））；对于（i=1，#f~，f[i，2]=最大值（0，f[i，2]-1））；系数回退（f）；}；

%o A065769（n）={my（f=因子（n>>估值（n，2）））[，1]~）；（A003557（n）*因子回复（向量（#f，i，预素数（f[i]-1）））；}；\\_Antti Karttunen_，2017年12月31日

%o（方案）（定义（A065769 n）（*（A003557 n）（A064989（A007947 n）））_Antti Karttunen_，2017年12月31日

%Y参见A000040、A000079、A003557、A007947、A064989、A065770、A151799。

%K mult、nonn、easy

%氧1,3

%A _Henry Bottomley_，2001年11月19日