Description
FGD正在破解一段密码,他需要回答很多类似的问题:对于给定的整数a,b和d,有多少正整数对x,y,满足x<=a
,y<=b,并且gcd(x,y)=d。作为FGD的同学,FGD希望得到你的帮助。
Input
第一行包含一个正整数n,表示一共有n组询问。(1<=n<= 50000)接下来n行,每行表示一个询问,每行三个
正整数,分别为a,b,d。(1<=d<=a,b<=50000)
Output
对于每组询问,输出到输出文件zap.out一个正整数,表示满足条件的整数对数。
Sample Input
2
4 5 2
6 4 3
4 5 2
6 4 3
Sample Output
3
2
//对于第一组询问,满足条件的整数对有(2,2),(2,4),(4,2)。对于第二组询问,满足条件的整数对有(
6,3),(3,3)。
2
//对于第一组询问,满足条件的整数对有(2,2),(2,4),(4,2)。对于第二组询问,满足条件的整数对有(
6,3),(3,3)。
反演板子。。刷通过数用的
代码:
1 #include<iostream> 2 #include<cstdio> 3 #include<cstring> 4 #define M 50010 5 using namespace std; 6 7 int n,tot; 8 int P[M],mul[M],sum[M]; 9 bool prime[M]; 10 11 void get_mul() 12 { 13 mul[1]=1; 14 for(int i=2;i<=50000;i++) 15 { 16 if(!prime[i]) 17 { 18 P[++tot]=i; 19 mul[i]=-1; 20 } 21 for(int j=1;j<=tot&&P[j]*i<=50000;j++) 22 { 23 prime[P[j]*i]=true; 24 if(i%P[j]==0) 25 { 26 mul[P[j]*i]=0; 27 break; 28 } 29 else mul[P[j]*i]=-mul[i]; 30 } 31 } 32 for(int i=1;i<=50000;i++) sum[i]=sum[i-1]+mul[i]; 33 } 34 35 int cal(int n,int m,int k) 36 { 37 int ans=0; 38 n/=k; m/=k; 39 for(int l=1,r;l<=n&&l<=m;l=r+1) 40 { 41 r=min(n/(n/l),m/(m/l)); 42 ans+=(n/l)*(m/l)*(sum[r]-sum[l-1]); 43 } 44 return ans; 45 } 46 47 int main() 48 { 49 get_mul(); 50 scanf("%d",&n); 51 for(int i=1;i<=n;i++) 52 { 53 int a,b,c,d,k; 54 scanf("%d%d%d",&a,&b,&k); 55 printf("%d\n",cal(a,b,k)); 56 } 57 return 0; 58 }