|
terms[p_]：=系数规则[p，x/@最大范围[{0}，案例[p，x[t_]：>t，\[Infinity]]]；
cycleindSymm[n_]：=循环符号[n]=术语@CycleIndexPolynomial[对称组[n]，x/@范围[n]]；
cycleindEdge[0|1]=1；
cycleindEdge[n_]：=循环边缘[n]=术语@总和[Last[t]使用[{tt=First[t]}，使用[{ind=Select[范围@n，tt[[#]]！=0 &]},
乘积[x[LCM@@p]^（GCD@@p倍@@tt[[p]]），{p，子集[ind，{2}]}]
乘积[With[{e=tt[[k]]}，x[k]^（ke（e-1）/2+商[k-1，2]e）If[EvenQ[k]，x[k/2]^e，1]]，{k，ind}]
]]，{t，循环符号[n]}]；
v[1，_]=v[_，1]=1；
v[n_，m_]：=[Last[a]Last[b]与[{aa=之和第一个@a，bb=第一个@b}，乘积[Sum[vb bb[[vb]]，{vb，交集[Divisors[va]，范围@m]}]^aa[[va]]，｛va，选择[范围@长度@aa，aa[[#]]！=0&]}]]，{b，循环符号[m]}，{a，循环边[n]}]；
a[n]：=1+v[n，-1+n（n-1）/2]；
| 