%I#45 2017年3月27日10:25:16

%S 1,0,0，-1,0,0,0，

%T-54,0,0,0，-5，0,0,10,30,0,0-0,00,0.0，0,0114，-1，-2,0,,0,6,0126,0,00,0，

%U 0,0,0，-4470,0,0,1,0,,0,0252,0,0-0,00,0_0,0

%N序列由Sum_{d|N}a（d）^（N/d）=mu（N）^2的条件定义，其中mu（N）是Möbius函数。

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

%F a（1）=1。

%如果p>2，a（p）=0，a（p^2）=-1，a（p^n）=0表示n>2。

%F a（p1*p2*..*pn）=0，a（2*p1*…*pn。

%F如果p1，。。。，pn很奇怪，似乎：

%F a（p1^2*p2^2*…*pn^2）=（-1）^n，

%F a（p1^k1*p2^k2*…*pn^kn）=0，如果k1，。。。，kn>2，

%F a（2*p1^k1*p2^k2*…*pn^kn）>0，如果k1，。。。，kn>1。

%F a（2^n）=A264610（n）。

%e对于素数p，a（p）^1+a（1）^p=mu。

%e对于n=6，a（1）^6+a（2）^3+a（3）^2+a（6）^1=mu。

%p a:=proc（n）选项记忆；numtheory：-mobius（n）^2-add（procname（n/d）^d，d=“minus”（numtheory：-除数（n），{1}））结束proc；a（1）：=1；La:=seq（a（i），i=1。。100）

%t a[n_]：=如果[n<2，1，MoebiusMu[n]^2-总和[如果[d==1，0，a[n/d]^d]，{d，除数[n]}]；表[a[n]，{n，100}]（*_Indranil Ghosh_，2017年3月26日*）

%o（PARI）a（n）=如果（n==1，1，moebius（n）^2-sumdiv（n，d，如果（d==1，0，a（n/d）^d））；

%o（哈斯克尔）

%o a263774 1=1

%o a263774 n=foldl（-）（a008966 n）$zipWith（^）（映射a'$reversed ds）ds

%o其中a'x=如果x==n，则0为a263774 x

%o ds=a027750_低n

%o--_Reinhard Zumkeller_，2015年12月6日

%Y参考A008683，A264610。

%Y参见A008966、A027750。

%K符号

%O 1,8型

%2015年11月28日，亚美尼亚