%I#16 2021年5月9日09:51:41

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

%电话18,19,12,20,21,22,13,23,24,25,14,26,27,28,15,16,29,30,17,18,19,31,20，

%U 32,21,33,22,34,35,36,23,37,38,24,25,39,40,41,26,42,43,44,27,45,46,28,29,47,48,49,30,31,50,51,32,52,53,54,33,55,34,56,35,58,59,36

%N a（1）=0；对于n>1，如果n不是平方自由的，则a（n）=A057627（n），否则a（n）=A013928（n）。

%C每个数字n>0在这个序列中正好出现两次，位置A005117（1+n）和A013929（n）。

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

%F a（1）=0；对于n>1，如果A008683（n）为0[当n不为平方自由]，则a（n）=A057627（n），否则a（n）=A013928（n）。

%o（方案）（定义（A284584 n）（条件（（=1 n）0）（（零？（A008683 n））（A057627 n）））（否则（A013928 n））

%o（Python）

%o来自sympy import mobius

%o来自sympy.theory.factor导入核心

%o定义a057627（n）：返回n-和（[mobius（k）**2表示范围（1，n+1）中的k）]

%o def a013928（n）：返回和（[1代表范围（1，n）中的i，如果核心（i）==i]）

%o定义a（n）：

%o如果n==1：返回0

%o如果核心（n）==n：返回a013928（n）

%o其他：返回a057627（n）

%o打印（[a（n）代表范围（121）中的n）]#_Indranil Ghosh，2017年4月17日

%Y参见A005117、A008683、A013928、A013929、A057627、A243343、A243354、A24334、A24335。

%Y参考A066136（类似序列）。

%K nonn公司

%氧1,3

%2017年4月17日