P2231-[HNOI2002]跳蚤【容斥】
生活随笔
收集整理的這篇文章主要介紹了
P2231-[HNOI2002]跳蚤【容斥】
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
正題
題目鏈接:https://www.luogu.com.cn/problem/P2231
題目大意
求一個由[1,m][1,m][1,m]的整數組成的長度為nnn的序列使得他們的gcdgcdgcd和mmm互質。
解題思路
考慮容斥減去不合法的答案。那就是要求序列的gcdgcdgcd和mmm不互質的個數,那么我們依舊需要容斥計算這個問題,顯然一個約數xxx的貢獻就是(mx)n(\frac{m}{x})^n(xm?)n,設fif_ifi?表示第iii個約數的容斥系數為?∑dj∣difi-\sum_{d_j|d_i}f_i?∑dj?∣di??fi?。
這樣理論時間復雜度就是約數個數的平方,加上快速冪,也就是接近于O(m2log?n)O(\sqrt m^2\log n)O(m?2logn)。然而時間復雜度很小。
codecodecode
#include<cstdio> #include<cstring> #include<algorithm> #define ll long long using namespace std; ll n,m,q[1100000],cnt,ans,f[1100000]; ll power(ll x,ll b){ll ans=1;while(b){if(b&1)ans=ans*x;x=x*x;b>>=1;}return ans; } int main() {scanf("%lld%lld",&n,&m);for(ll i=2;i*i<=m;i++){if(m%i==0){q[++cnt]=i;if(i*i!=m)q[++cnt]=m/i;}}q[++cnt]=m; sort(q+1,q+1+cnt);for(ll i=1;i<=cnt;i++){f[i]=1;for(ll j=1;j<i;j++)if(q[i]%q[j]==0)f[i]-=f[j];ans+=power(m/q[i],n)*f[i];}printf("%lld",power(m,n)-ans); }總結
以上是生活随笔為你收集整理的P2231-[HNOI2002]跳蚤【容斥】的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 微信号怎么改 微信号改变步骤
- 下一篇: 什么是pid调节器