[bzoj 1101--Poi2007]Zap

FGD正在破解一段密码，他需要回答很多类似的问题：对于给定的整数a,b和d，有多少正整数对x,y，满足x<=a，y<=b，并且gcd(x,y)=d。作为FGD的同学，FGD希望得到你的帮助。

莫比乌斯反演的处女题。其实还是较水的，首先先把a%d，b%d，这样问题便简化成求gcd(x,y)=1的对数了。f(i)表示的是gcd(x,y)=i的对数，F(i)表示的是gcd(x,y)%i=0的对数，那问题又简化成求f(1)。
熟悉莫比乌斯的同学明白在知道f与F的含义后，这题便可做了。

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<algorithm>
using namespace std;
int pr=0,prime[51000];
bool v[51000];
long long Mu[51000],sM[51000];
void get_Mu()
{
    memset(v,true,sizeof(v));
    Mu[1]=sM[1]=1;
    for(int i=2;i<=50000;i++)
    {
        if(v[i]==true)
        {
            prime[++pr]=i;
            Mu[i]=-1;
        }
        for(int j=1;(j<=pr && i*prime[j]<=50000);j++)
        {
            v[i*prime[j]]=false;
            if(i%prime[j]==0)
            {
                Mu[i*prime[j]]=0;
                break;
            }
            Mu[i*prime[j]]=-Mu[i];
        }
        sM[i]=sM[i-1]+Mu[i];
    }
}
long long find(int x,int y)
{
    if(x>y)swap(x,y);
    int last=0;long long s=0;
    for(int i=1;i<=x;i=last+1)
    {
        last=min(x/(x/i),y/(y/i));
        s+=(sM[last]-sM[i-1])*(x/i)*(y/i);
    }
    return s;
}
int main()
{
    get_Mu();
    int t;
    scanf("%d",&t);
    while(t--)
    {
        int a,b,d;
        scanf("%d%d%d",&a,&b,&d);
        if(d==0)
        {
            printf("0\n");
            continue;
        }
        long long ans=0;
        ans=find(a/d,b/d);
        printf("%lld\n",ans);
    }
    return 0;
}