%I#12 2019年4月14日07:51:31
%S 0,1,1,-1,-2,2,6,-5,-18,15,59,-54,-215199813,-744,-3135289012394,
%电话:-11538,-5001746806204893,-192451,-8496818009743560927,-3367656,
%电话:-150584781427942664171736,-60992032,-2753046652621990501188070488,-1133572891,-5153913606
%N G.f.A(x)满足:A(x)=x*exp(-A(-x)+A(-x^2)/2-A(-x^3)/3+A(-x^4)/4-A(-x^5)/5+…)。
%F G.F.:A(x)=和{n>=1}A(n)*x^n=x*产品{n>=1}1/(1+x^n)^(-1)^n*A(n))。
%F递归:a(n+1)=(1/n)*和{k=1..n}(和{d|k}(-1)^(k/d+d)*d*a(d))*a(n-k+1)。
%通用公式:A(x)=x+x^2-x^3-2*x^4+2*x^5+6*x^6-5*x^7-18*x^8+15*xs^9+59*x^10-54*x^11-215*x^12+。。。
%t项=36;A[_]=0;Do[A[x_]=x Exp[Sum[(-1)^k A[-x^k]/k,{k,1,terms}]]+O[x]^(terms+1)//正常,terms+1];系数列表[A[x],x]
%t a[n_]:=a[n]=系列系数[x乘积[1/(1+x^k)^(-1)^k a[k]),{k,1,n-1}],{x,0,n}];a[0]=0;表[a[n],{n,0,36}]
%Y参考A000081、A004111、A045648、A049075、A073075、A115593、A307365、A307366、A307538。
%K符号
%O 0.5
%A_Ilya Gutkovskiy_,2019年4月14日