已知n,dn,dn,d,它们2个同阶,求 ∑k=0d−1(1+ωdk)n\sum_{k=0}^{d-1}(1+\omega_d^k)^n∑k=0d−1(1+ωdk)n 的值,模一个素数p,如果你想的话可以是NTT模数,我想找一个O(n)O(\sqrt n)O(n)(或更好)的算法,目前已有O(nlog2(n))O(\sqrt n\log^2(n))O(nlog2(n))的方法,做了好几天了qwq