简单的容斥原理

容斥原理：奇加偶减

hdu2588
本以为只要求一下GCD就行了。。。结果果然是我天真！
由题可知两个数里面有一个已经被固定了，那么他的gcd一定是n的一个因子，所以我们可以把n的所有因子全部找出来，排序以后直接找大于m的所有因子，由于可能出现因子之间有倍数关系，因此再去重，然后直接使用容斥原理就行了！！

容斥原理：使用一个二进制数字1表示使用过，0表示不用，然后枚举每一位是否为1直接可得每种方案。

#include<stdio.h>
#include<math.h>
#include<algorithm>
using namespace std;
const int N=2e5+5;
int main()
{
    int n,m,i,j,k,t;
    int a[N],cnt=0,b[N],cnt1=0;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d %d",&n,&m);
        cnt=0,cnt1=0;
        int num=sqrt(n)+0.5;
        for(i=1;i<=num;i++)
        {
            if(n%i==0)
            {
                //printf("%ddd\n",i);
                a[++cnt]=i;
                if(n/i!=i) a[++cnt]=n/i;
            }
        }
        sort(a+1,a+1+cnt);
        int x=lower_bound(a,a+1+cnt,m)-a;
        for(i=x;i<=cnt;i++)
        {
            if(i==1||a[i]!=1) b[++cnt1]=a[i];
            else continue;
            for(j=i+1;j<=cnt;j++)
            {
                if(a[j]%a[i]==0) a[j]=1;
            }
        }
        int ans=0;
        for(i=1;i<(1<<(cnt1));i++)
        {
            int s=0;
            long long x=1;
            for(j=0;j<cnt1;j++)
            {
                if(i&(1<<j))
                {
                    x=(x*b[j+1])/__gcd(x,(long long)b[j+1]);//找最小公倍数
                    s++;
                }
            }
            if(s&1) ans+=n/x;
            else ans-=n/x;
        }
        printf("%d\n",ans);
    }
}