日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 > 编程资源 > 编程问答 >内容正文

编程问答

P2231-[HNOI2002]跳蚤【容斥】

發布時間:2023/12/3 编程问答 45 豆豆
生活随笔 收集整理的這篇文章主要介紹了 P2231-[HNOI2002]跳蚤【容斥】 小編覺得挺不錯的,現在分享給大家,幫大家做個參考.

正題

題目鏈接:https://www.luogu.com.cn/problem/P2231


題目大意

求一個由[1,m][1,m][1,m]的整數組成的長度為nnn的序列使得他們的gcdgcdgcdmmm互質。


解題思路

考慮容斥減去不合法的答案。那就是要求序列的gcdgcdgcdmmm不互質的個數,那么我們依舊需要容斥計算這個問題,顯然一個約數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]跳蚤【容斥】的全部內容,希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯,歡迎將生活随笔推薦給好友。