根sof1（p，n）={本地（a）；/*奇素数p的a^（p-1）%p^n==1的本原解*/a=向量（（p-3）/2，i，i+1）；对于（j=2，n，a=向量（#a，i，升力（Mod（升力（Mod（a[i]，p^j）^（p-1）））\p^（j-1），p）/Mod；)；向量排序（concat（[a，向量（#a，i，p^n-a[i]），[1，p^n-1]]）}A125609（n，p=3）={局部（v，i）；v=根sof1（p，n）；i=0；而（1，对于（j=1，#v，如果是（i素数（i*p^n+v[j]），则返回（i*p ^n+v[j））；)；i++；)；}