HDU2841 （队列容斥）

#include<bits/stdc++.h>
using namespace std;
__int64 n;
__int64 a,b;
__int64  slove( __int64 m,__int64 gg)
{
    __int64  que[1000];
    __int64  a[1000];
    memset(que,0,sizeof(que));
    memset(a,0,sizeof(a));
    __int64 sum=0;
    __int64 t=0;
    __int64 ss=0;
    for(__int64 i=2;i*i<=m;i++)
    {
        if(m%i==0)
        {
            que[ss++]=i;
            while(m%i==0)
                m=m/i;
        }
    }
    if(m>1)
    que[ss++]=m;
    a[t++]=-1;
    for(__int64 i=0;i<ss;i++)
    {
         int k=t;
         for(__int64 j=0;j<k;j++)
         {
             a[t++]=que[i]*a[j]*(-1);
         }
    }
    for(__int64 i=1;i<t;i++)
       sum=sum+gg/a[i];
    return sum;
}
int main()
{
    while(scanf("%I64d",&n)!=EOF)
    {
        for(__int64 i=1;i<=n;i++)
        {
            __int64 re=0;
            scanf("%I64d %I64d",&a,&b);
            for(__int64 j=1;j<=a;j++)
                re+=slove(j,b);
            printf("%I64d\n",a*b-re);
        }
    }
    return 0;
}