首页 > 其他 > 详细

bzoj 2005

时间:2016-12-10 08:59:11      阅读:128      评论:0      收藏:0      [点我收藏+]

    裸的2D gcd。ans=(Σ(d<=n)phi[d]*(n/d)*(m/d))*2-n*m;

    

 1 #include<iostream>
 2 #include<cstdio>
 3 #include<algorithm>
 4 #include<cstring>
 5 #define int long long
 6 using namespace std;
 7 int n,m;
 8 int pr[100005],su[100005],cnt,phi[100005];
 9 void shai()
10 {
11     phi[1]=1;
12     for(int i=2;i<=100000;i++)
13     {
14         if(!phi[i])su[++cnt]=i,pr[i]=i,phi[i]=i-1;
15         for(int j=1;su[j]<=pr[i]&&su[j]*i<=100000&&j<=cnt;j++)
16         {
17               pr[su[j]*i]=su[j];
18               if(su[j]==pr[i])phi[su[j]*i]=phi[i]*su[j];
19               else phi[su[j]*i]=phi[i]*(su[j]-1);
20         }
21     }
22 }
23 signed main()
24 {
25    shai();
26    scanf("%lld%lld",&n,&m);int ans=0;
27    for(int i=1;i<=min(n,m);i++)ans+=phi[i]*(n/i)*(m/i);
28    ans*=2;ans-=(n*m);
29    printf("%lld\n",ans);
30    return 0;
31 }

 

     

bzoj 2005

原文:http://www.cnblogs.com/ezyzy/p/6154769.html

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!