行=10；
col[k_]：=col[k]=系数列表[乘积[（1-j^（k*j）*x^j），{j，1，行+3}]+O[x]^（行+3），x]；
A[n_，k_]：=列[k][[n+1]]；
（*或：*）
A[0，_]=1；A[n_，k_]：=A[n，k]=-（1/n）*和[DivisorSum[j，#^（1+k*j）&]*A[n-j，k]，{j，1，n}]；
表[A[n-k，k]，{n，0，行-1}，{k，n，0和-1}]//展平(*Jean-François Alcover公司2017年11月10日*）
