（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）））
