hdu2841(2009多校第三场) gcd(x,y)==1对数

关键是要能看出来是gcd==1个数===

然后就是素数分解容斥定理了。。

 1 #include<stdio.h>
 2 #include<math.h>
 3 #include<string.h>
 4 int k,n,num;
 5 int vis[100005],prime[100005],a[10];
 6 void dfs(int now,int flag,int bei)
 7 {
 8   int i;
 9   if (flag) k+=n/bei;
10   else k-=n/bei;
11   for (i=now+1;i<=num;i++)
12     dfs(i,1-flag,bei*a[i]);
13 }
14 int main()
15 {
16   int cnt,i,m,j,T,tmp,sqr;
17   long long ans;
18   memset(vis,0,sizeof(vis));
19   for (i=2,cnt=0,vis[1]=1;i<=100000;i++)
20     if (vis[i]==0)
21     {
22       prime[++cnt]=i;
23       for (j=2;i*j<=100000;j++) vis[i*j]=1;
24     }
25   scanf("%d",&T);
26   while (T--)
27   {
28     scanf("%d%d",&m,&n);
29     if (m>n) {tmp=m; m=n; n=tmp; }
30     ans=0;
31     for (i=1;i<=m;i++)
32     {
33       num=0; k=i; sqr=(int)sqrt(1.0*k);
34       for (j=1;j<=cnt&&prime[j]<=sqr&&prime[j]<=k;j++)
35         if (k%prime[j]==0)
36         {
37           a[++num]=prime[j];
38           while (k%prime[j]==0) k/=prime[j];
39         }
40       if (k!=1) a[++num]=k;
41       for (j=1,k=0;j<=num;j++)
42         dfs(j,1,a[j]);
43       ans=ans-k+n;
44     }
45     printf("%I64d\n",ans);
46   }
47   return 0;
48 }
View Code

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=2841

posted on 2014-12-10 13:31  xiao_xin  阅读(193)  评论(0编辑  收藏  举报

导航