%I#21 2017年11月10日05:55:10
%S 1,1,-1,1,-1,-1,1,-1,-4,0,1,-1-,-16,-23,0,1.,-1,-64,-713,-229,1,1,-1,
%电话-256、-19619、-64807、-2761、0,1、-1、-1024、-531185、-16757533、-9688425、,
%U-42615、1、1、-1、-4096、-14347883、-4294435855、-304995411977、-2165979799、-758499、0
%N平方数组A(N,k),N>=0,k>=0(反对偶),其中k列是乘积{j>=1}(1-j^(k*j)*x^j)的x次幂展开式。
%H Seiichi Manyama,反对角线n=0..52,扁平</a>
%F A(0,k)=1和A(n,k)=-(1/n)*求和{j=1..n}(求和{d|j}d^(1+k*j))*A(n-j,k)对于n>0。
%e方阵开始:
%e 1,1,1,1。。。
%e-1,-1,-1,-1-,-1。。。
%e-1,-4,-16,-64,-256。。。
%e 0,-23,-713,-19619,-531185。。。
%电子邮箱0,-229,-64807,-16757533,-4294435855。。。
%t行=10;
%t col[k_]:=col[k]=系数列表[乘积[(1-j^(k*j)*x^j),{j,1,行+3}]+O[x]^(行+3),x];
%t A[n_,k_]:=列[k][[n+1]];
%t(*或:*)
%tA[0,_]=1;A[n_,k_]:=A[n,k]=-(1/n)*和[DivisorSum[j,#^(1+k*j)&]*A[n-j,k],{j,1,n}];
%t表[A[n-k,k],{n,0,rows-1},{k,n,0-1}]//扁平(*_Jean-François Alcover_,2017年11月10日*)
%Y列k=0..1表示A010815、A292312。
%Y行n=0..2表示A000012,(-1)*A000012、(-1)*A000302。
%Y参见A283675、A294758。
%K符号,表格
%0、9
%A _Seiichi Manyama,2017年11月6日
|