%I#28 2023年10月25日14:59:00
%S 1,0,2,1,4,3,8,7,15,15,27,29,48,53,82,94137160225265362430572,
%电话:68389210661370164020782487311737254624519679180929885,
%电话:117521426316922204162416729073425440921482135734567409
%N a(N)=和{k=0..N}(-1)^(N-k)*A000041(k)。
%C基本上是A024786的第一个差异(参见公式)。此外,a(n)是n+2分区集最后一段中的2个数(见A135010)_Omar E.Pol_,2008年9月10日
%H Vaclav Kotesovec,<a href=“/A0877787/b087787.txt”>n表,n=0..10000时的a(n)</a>
%F G.F.:1/(1+x)*1/产品{k>0}(1-x^k)。
%F a(n)=1/n*和{k=1..n}(σ(k)+(-1)^k)*a(n-k)。
%F a(n)=A024786(n+2)-A024786(n+1)_Omar E.Pol_,2008年9月10日
%F a(n)~exp(Pi*sqrt(2*n/3))/(8*sqert(3)*n)*(1+(11*Pi/(24*sqort(6))-sqrt(3/2)/Pi)/sqrt(n)-(11/16+(23*Pi^2)/6912)/n)_瓦茨拉夫·科特索维奇,2016年11月5日
%F a(n)=A000041(n)-a(n-1).-_乔恩·麦加(Jon Maiga),2019年8月29日
%t表[总和[(-1)^(n-k)*分区P[k],{k,0,n}],{n,0,50}](*_Vaclav Kotesovec_,2015年8月16日*)
%t(*更高效的程序*)sig=1;su=1;扁平[{1,表[sig=-sig;su=su+sig*PartitionsP[n];Abs[su],{n,1,50}]}](*_Vaclav Kotesovec_,2016年11月6日*)
%o(Python)
%o来自sympy导入npartitions
%o def A087787(n):返回和(-npartitions(k)if n-k&1 else npartitions(k)for k in range(n+1))#_Chai Wah Wu_,2023年10月25日
%Y参考A000041、A024786、A135010、A138121、A141285。
%K nonn公司
%0、3
%A _弗拉德塔·乔沃维奇,2003年10月7日
|