(PARI)块大小受限的分区数
NCPartitionsModCyclic(v)={my(n=#v);
my(p=serreverse(x/(1+总和(k=1,#v,x^k*v[k]))+O(x^2*x^n))/x);
my(vars=变量(p));
my(varpow(r,d)=substvec(r+O(x^(n\d+1)),变量,应用(t->t^d,变量));
my(q=x*deriv(p)/p);
my(T=总和(k=1,#v,my(T=v[k]));if(t,x^k*t*sumdiv(k,d,eulerphi(d)*varpow(p,d)^(k/d))/k));
T+2+整数(和(d=1,n,eulerphi(d)*varpow(q,d))/x)-p
}
Vec(NCPartitionsModCyclic(向量(40,k,isprime(k)))