%I#23 2020年6月26日23:44:19

%S 0,1,0,2,0,1,0,3,0,1,4,0,12,4,0,1,2,0,1,6,2,1,0,5,0,1,10,0,3,1,5,1，

%温度4,2,0,1,2,7,0,3,0,18,4,1,0,6,0,12,34,0,4,11,0,8,0,1,8,6,4,0，

%U 66,2,5,0,3,0,1,2130,8,3,0，8,0,12,4,1,2,19,0,5,8258,2,1,4,7,0,16,2,0,3,15,6

%N a（1）=0；对于k>=2，a（素数（k））=0，a（k^2）=2*a（k）；否则a（n）=a（A334870（n））+1。

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

%H<a href=“/index/Pri#prime_indices”>根据素因式分解中的索引计算出的序列的索引项</a>

%F替代定义：（开始）

%F a（1）=0，a（2）=1；否则，对于n=k*m^2，k平方自由：

%F如果m=1，a（n）=A048675（A052126（k））；

%如果m>1，则a（n）=A048675（k）+2*a（m）。

%F（结束）

%F对于n=4*A122132（k），a（n）=A048675（n）。

%F更一般地说，当且仅当n在A335738中时，a（n）=A048675（n）。

%F a（n）=A335426（A225546（n））。

%F a（A003961（2k+1））=2*a（2k/1）。

%F如果n在A036554中，a（n）=a（n/2）+1；否则，对于n<>3，a（n）=2*a（A019565（k/2）*m^2）-a（m^2），其中n=A019565（k）*m^2。

%o（PARI）

%o A334870（n）=如果（issquare（n），sqrtint（n）），my（c=芯（n）、m=n）；对于素数（p=2，如果（！（c%p），m/=p；断裂，m*=p））；（m） ）；

%o A335427（n）=如果（n<=2，n-1，if（i质数（n），0，if）（发行方（n）、2*A335427-平方（n））、1+A335427.（A334870（n）；

%o（PARI）

%o A048675（n）={my（f=因子（n））；和（k=1，#f~，f[k，2]*2^素数（f[k、1]））/2；}；

%o A052126（n）=如果（1==n，n，（n/vecmax（系数（n）[，1]））；

%o A335427（n）=如果（n<=2，n-1，如果（无平方（n），A048675（A052126（n）），my（k=核心（n））；A048675（k）+2*A335427（平方（n/k））；

%Y A052126、A225546、A334870、A335426用于定义此序列的公式中。

%Y相关全加性序列：A048675。

%Y参见A062090（零值指数）、A003159（偶数值指数）和A036554（奇值指数）。

%Y参考A005117、A122132、A334872、A335738。

%Y A003961、A019565用于表示此序列项之间的关系。

%K nonn公司

%O 1,4个

%2020年6月15日，安提·卡图宁和佩特·穆恩