多多色-多人伦交性欧美在线观看-多人伦精品一区二区三区视频-多色视频-免费黄色视屏网站-免费黄色在线

國內最全IT社區平臺 聯系我們 | 收藏本站
阿里云優惠2
您當前位置:首頁 > php開源 > 綜合技術 > BZOJ 2818 gcd(莫比烏斯反演)

BZOJ 2818 gcd(莫比烏斯反演)

來源:程序員人生   發布時間:2015-07-27 07:53:01 閱讀次數:3298次
Gcd
Time Limit:10000MS     Memory Limit:262144KB     64bit IO Format:%lld & %llu
Submit Status

Description

給定整數N,求1<=x,y<=N且Gcd(x,y)為素數的
數對(x,y)有多少對.

Input

1個整數N

Output

如題

Sample Input

4

Sample Output

4


#include<cstdio> #include<algorithm> #include<iostream> #include<cmath> #include<cstring> #include<vector> #define ll long long #define N 10000010 using namespace std; bool check[N+10]; int prime[N]; int mu[N+10]; int tot; int n; void Moblus() { memset(check,0,sizeof check); mu[1]=1; tot=0; for(int i=2; i<N; i++) { if(!check[i]) { prime[tot++]=i; mu[i]=⑴; } for(int j=0; j<tot; j++) { if(i*prime[j]>N)break; check[i*prime[j]]=1; if(i%prime[j]==0) { mu[i*prime[j]]=0; break; } else mu[i*prime[j]]=-mu[i]; } } } int main() { Moblus(); while(scanf("%d",&n)!=EOF) { ll ans=0; for(int i=0; i<tot; i++) { int x=n/prime[i]; if(x==0)break; for(int j=1; j<=x; j++) ans+=(ll)mu[j]*(x/j)*(x/j); } printf("%lld ",ans); } return 0; }


生活不易,碼農辛苦
如果您覺得本網站對您的學習有所幫助,可以手機掃描二維碼進行捐贈
程序員人生
------分隔線----------------------------
分享到:
------分隔線----------------------------
關閉
程序員人生
主站蜘蛛池模板: 在线一区观看 | 国产精品欧美一区二区三区不卡 | 不卡视频一区二区三区 | 亚洲视频在线播放 | 国产精品一区欧美日韩制服 | a级片在线免费看 | 亚洲天堂久久 | 日本一级高清不卡视频在线 | 一本久草 | 亚洲欧美日韩第一页 | 亚洲第成色999久久网站 | 欧美福利一区二区三区 | 伊人中文在线 | www.性欧美| 在线中文字幕网站 | 中文字幕人成不卡一区 | 国产一区二区三区欧美精品 | 成人在线视频国产 | 一级毛片女人18水真多 | 国产欧美日韩在线一区二区不卡 | 午夜大片免费男女爽爽影院久久 | 狠狠躁天天躁夜夜躁婷婷 | 亚洲欧美视频一区二区三区 | 免费观看的男女羞羞视频 | 操你啦网站 | 成人在线不卡视频 | 美国一级黄色毛片 | 日本一区二区三区在线网 | 国产成人一区二区三区影院免费 | xxxxxxx国产精品视频 | 亚洲国产精品成 | 日本免费中文字幕在线看 | 国产区成人精品视频 | 另类综合图片 | 国产真人毛片一级视频 | 成人不卡 | 国产美女福利在线 | 日本美女影院 | 欧美激情bbbbbxxxxⅹ | 国产欧美自拍 | 亚洲第一成人在线 |