HDU1695 GCD 莫比乌斯反演

题目地址：http://acm.hdu.edu.cn/showproblem.php?pid=1695
题目大意就是在(1,b)中取x，在(1,d)中取y，求gcd(x,y)=k的个数。
本题使用的算法为莫比乌斯反演。
思路如下：
gcd(x,y)=k---->gcd(x/k,y/k)=1,x/k<=b/k;y/k<=d/k;

我们选择写法二，然后看下面 　　　当时一直以为是写法一，研究了半天，也没有研究出来，哭

上面这个图片是从我同学的PPT上拿下来的，侵权即删
这道题的话还有一个点，就是gcd(5,7)和gcd(7,5)，他视作是一个数，所以我们需要把重复的个数去掉。gcd(x,y)=d，分为x<min(n,n)&&y<min(n,m)和x>min(n,m)||y>min(n,m)两大部分，前一部分会有重合，后一部分没有。大家可以在纸上写写看看
重复的数可以理解为在1<=x<=min(n,m),1<=y<=min(n,m)相同部分的f(d)，我们再求一遍然后减去他的二分之一即可,即f(d)(n,m)-0.5*(d)(min(n,m),min(n,m)) //f(d)(x,y)，x，y分别为gcd(a,b)=k中a和b的取值上限
为什么要减1/2呢 f(d)(min(n,m),min(n,m))是奇数 ，有gcd(d,d)和 x个gcd(n1,n2) gcd（n2,n1)构成 ，即0.5*f(d)(min(n,m),min(n,m)) =x.

#include <iostream>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#define  maxn 100009
using namespace std;
int prime[maxn],cnt,mo[maxn],sum[maxn];///mo为莫比乌斯函数
bool f[maxn];
void data()///线性筛莫比乌斯函数
{
    cnt=0;
    memset(f,true,sizeof(f));
    f[0]=f[1]=false;
    mo[1]=1;
    for(int i=2;i<maxn;i++)
    {
        if(f[i])
        {
            prime[cnt++]=i;
            mo[i]=-1;
        }
        for(int j=0;j<cnt&&i*prime[j]<maxn;j++)
        {
            f[i*prime[j]]=false;
            if(i%prime[j])
            mo[i*prime[j]]=-1*mo[i];
            else
            {
                mo[i*prime[j]]=0;
                break;
            }
        }
    }
}
int main()
{
    data();
    int t,a,b,c,d,k,x,min1,max1;
    scanf("%d",&t);
   //printf("%d***\n",t);
        for(x=1;x<=t;x++)
        {
            scanf("%d%d%d%d%d",&a,&b,&c,&d,&k);
            //printf("%d%d%d%d%d\n",a,b,c,d,k);
            printf("Case %d: ",x);
            if(k==0)///不要漏了k为零
            {
            cout<<"0"<<endl;
            continue;
           }
            b=b/k;d=d/k;
            min1=min(b,d);
             // printf("%d %d %d***\n",b,d,min1);
            long long sum1,sum2;
            sum1=sum2=0;
            for(int i=1;i<=min1;i++)
            {
                sum1+=(long long)mo[i]*(b/i)*(d/i);
            }
            for(int i=1;i<=min1;i++)
            {
                sum2+=(long long)mo[i]*(min1/i)*(min1/i);
            }
           // printf("%lld\n",sum1);
            sum1=sum1-sum2/2;
            printf("%lld\n",sum1);
        }

    return 0;
}