|
terms[p_]：=系数规则[p，x/@最大范围[{0}，案例[p，x[t_]：>t，\[Infinity]]]；
循环符号术语循环Symm[n_]：=循环符号术语循环Symm[无]=模块[{c（c）=条款@CycleIndexPolynomial[对称群[n]，x/@范围[n]]},如果[头部[c（c）] ===Plus（加）,列表@@c（c）, {c（c）}]];]];
系数[Times[c_？NumericQ，_]]：=c；系数[c_？NumericQ]：=c；系数[_]：=1；
indets[e_]：=并集[Cases[{e}，x[t_]：>t，\[Infinity]]]；
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（v）[变量n个_,米_] :=扁平化术语总和[最后[一]最后[变量] =连环体b条]使用[{aa公司=弗斯特@一,英国广播公司=弗斯特@表b条},产品[总和[vb格式 英国广播公司[[vb格式]], {vb格式,交叉[常量数组约数[v（v）va公司],范围@米]}]^aa公司[[va公司]], {va公司,指数选择[变量范围@长度@aa公司,aa公司[[#]] !=0&]}]], {b条,x个循环Symm[v（v）]]], {v（v）米]}, {一,指数循环边缘[变量n个]}];
循环编辑[0]=1；
cycleindEdg[n_]：=cycleindEdg[n]=Sum[Module[{cycs=扁平项[t]}，系数[t]乘积[x[LCM@@p]^GCD@@p，{p，子集[cycs，{2}]}]乘积[Cf[OddQ[l]，x[l]^商[l，2]，x[l/2]x[l]^（l/2-1）]，{l，cycs}]]，{t，cycleindSymmTerms[n]}]；
varintoPower[n_，za_]：=模块[{indvars=indets[za]}，总和[coeff[beta]（za/.表[x[v]->1+z^v总和[d指数[beta，x[d]]，{d，除数[v]}]，{v，indvars}]），{beta，cycleindSymmTerms[n]}]；
vgf[n_，q_]：=vgf[n，q]=展开[varintoPower[n，cycleindEdg[q]]；
v（v）一【n】_,q个_] :=系数[vgf公司1+v（v）[n个, -1+n个,q个],z（z）,q个(q个-n个-1)/2];
a[n]：=1+v[-1+n（n-1）/2，n]；
表[a[n]，{n，2，7}](*安德烈·扎博洛茨基,八月二月 06 282024,之后_马尔科 2022*)里德尔_ *)
| 