登录
序列由Sum_{d|n}a(d)^(n/d)=mu(n)^2的条件定义,其中mu(n)是Möbius函数。
1

%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,

%如果k1中的一个,F a(p1^k1*p2^k2*…*pn^kn)=0,。..,kn>2,

%如果k1中的一个,F a(2*p1^k1*p2^k2*…*pn^kn)>0,。..,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`(numtheoriy:-divisors(n),{1}))结束进程;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-总和(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型

%A _Gevorg Hmayakyan,2015年11月28日