（平价）
解剖ModCyclicRooted（v）={my（n=#v）；
my（q=向量（n））；q[1]=序列反转（x和（i=3，#v，x^i*v[i]）/x+O（x*x^n））；
对于（i=2，n，q[i]=q[i-1]*q[1]）；
my（vars=变量（q[1]））；
my（u（m，r）=substvec（q[r]+O（x^（n\m+1）），变量，应用（t->t^m，变量））；
my（p=O（x*x^n）+总和（i=3，#v，my（c=v[i]））；如果（c，c*sumdiv（i，d，eulerphi（d）*u（d，i/d）））；
向量（n，i，polceoff（p，i））}
{my（v=解剖ModCyclicRooted（应用（i->if（i>=3&i<=4，y^（i-3）+O（y^2）），[1..25]））；应用（p->polcoff（p，1），v[4..#v]）}\\安德鲁·霍罗伊德2017年11月22日
