|
（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）=子向量（r+O（x^（n\d+1）），变量，应用（t->t^d，vars））；
my（q=x*deriv（p）/p）；
my（T=和（k=1，#v，my（T=v[k]）；如果（T，x^k*T*sumdiv（k，d，eulerphi（d）*varpow（p，d）^（k/d））/k））；
形式上的（1+p+d）和（2）
}
Vec（NCPartitionsModCyclic（向量（40，k，isprime（k）））
| 