标签:gcd LuoguP4464 sum mid mu JZPKIL ans aligned 国家集训队
瞎jb推一波式子,推出
\[ans=n^{x+y}\sum_{j=0}^{y+1}a_j\sum_{i|n}\sum_{d_1|n/i}i^{j-x}d_1^{y-x}\mu(d_1)
\]
其中 \(a_j\) 是自然数幂和的多项式系数。
然后后面显然是个积性函数,随便求求就好了。
\[\begin{aligned}
\frac{ans}{n^y}&=\sum_{d\mid n}d^{x-y}\sum_{i=1}^{n}[gcd(i,n)=d]i^y\\
&=\sum_{d\mid n}d^{x-y}\sum_{j=1}^{n/d}[\gcd(j,n/d)=1](jd)^y\\
&=\sum_{d\mid n}d^{x}\sum_{j=1}^{n/d}[\gcd(j,n/d)=1]j^y\\
&=\sum_{d\mid n}d^{x}f(n/d)\\
\end{aligned}
\]
\[\begin{aligned}
f(n)&=\sum_{i=1}^{n}[\gcd(i,n)=1]i^y\\
&=\sum_{d|n}\mu(d)d^y\sum_{i=1}^{n/d}i^y
\end{aligned}
\]
\[\begin{aligned}
ans&=n^y\sum_{d_1\mid n}d_1^{x}\sum_{d_2\mid n/d_1}\mu(d_2)d_2^y\sum_{i=1}^{n/d_1/d_2}i^y\\
\end{aligned}
\]
\[\frac{ans}{n^y}=\sum_{i|n}^{n}g(i,y)\sum_{d1\times d2=n/i}d_1^x\mu(d_2)d_2^y=\sum_{i|n}\sum_{j=0}^{y}a_ji^j\sum_{d_1|n/i}^{n/i}(n/i/d_1)^x\mu(d_1)d_1^y=\sum_{j=0}^{y}a_j\sum_{i|n}i^j\sum_{d_1|n/i}\frac{n^x}{(id_1)^x}\mu(d_1)d_1^y
\]
\[g(n,k)=\sum_{i|n}\sum_{j|n/i}i^{k-x}j^{y-x}\mu(j)\\
g(p^c, k)=\sum_{i|p^c}\sum_{j|p^c/i}j^{k-x}\mu(i)i^{y-x}=\sum_{i|p^c}i^{k-x}-p^{y-x}\sum_{i|p^{c-1}}i^{k-i}
\]
标签:gcd,LuoguP4464,sum,mid,mu,JZPKIL,ans,aligned,国家集训队
来源: https://www.cnblogs.com/zcr-blog/p/14614670.html
本站声明:
1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。