---

題目

---

?

1101: [POI2007]Zap

Time Limit:?10 Sec??Memory Limit:?162 MB

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

Sample Output

3
2

HINT

?

對于第一組詢問，滿足條件的整數對有(2,2)，（2,4），（4，2）。對于第二組詢問，滿足條件的整數對有（6,3），（3,3）。

?

---

題解

---

由數論推知->

且只會有種取值，所以我們暴力就好了！

---

?

代碼

---

?

/*Author:WNJXYK*/ #include<cstdio> #include<iostream> #include<cstring> #include<string> #include<algorithm> #include<queue> #include<set> #include<map> using namespace std;#define LL long long #define Inf 2147483647 #define InfL 10000000000LLinline void swap(int &x,int &y){int tmp=x;x=y;y=tmp;} inline void swap(LL &x,LL &y){LL tmp=x;x=y;y=tmp;} inline int remin(int a,int b){if (a<b) return a;return b;} inline int remax(int a,int b){if (a>b) return a;return b;} inline LL remin(LL a,LL b){if (a<b) return a;return b;} inline LL remax(LL a,LL b){if (a>b) return a;return b;} inline void read(int &x){x=0;int f=1;char ch=getchar();while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}x=x*f;} inline void read(LL &x){x=0;LL f=1;char ch=getchar();while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}x=x*f;} inline void read(int &x,int &y){read(x);read(y);} inline void read(LL &x,LL &y){read(x);read(y);} inline void read(int &x,int &y,int &z){read(x,y);read(z);} inline void read(int &x,int &y,int &n,int &m){read(x,y);read(n,m);} inline void read(LL &x,LL &y,LL &z){read(x,y);read(z);} inline void read(LL &x,LL &y,LL &n,LL &m){read(x,y);read(n,m);}const int Maxn=50000; LL miu[Maxn+10]; inline void getMiu(){for (int i=1;i<=Maxn;i++){int target=i==1?1:0;int delta=target-miu[i];miu[i]=delta;for (int j=i+i;j<=Maxn;j+=i) miu[j]+=delta;} for (int i=2;i<=Maxn;i++) miu[i]+=miu[i-1]; }inline LL getAns(int n,int m){if (n>m) swap(n,m);LL Ans=0;int pos;for (int i=1;i<=n;i=pos+1){pos=remin(n/(n/i),m/(m/i));Ans+=(miu[pos]-miu[i-1])*(n/i)*(m/i);} return Ans; }int T; LL a,b,d; int main(){read(T);getMiu();for (;T--;){read(a,b,d);printf("%lld\n",getAns(a/d,b/d));}return 0; }

轉載于:https://www.cnblogs.com/WNJXYK/p/4063912.html

總結

以上是生活随笔為你收集整理的BZOJ 1101: [POI2007]Zap的全部內容，希望文章能夠幫你解決所遇到的問題。

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