http://acm.hdu.edu.cn/showproblem.php?pid=2841 题意很简单

感觉不错的一个题目 又学习了一点知识。
//注意数据:
//求n*m的矩形内,有多少对互质的数(x,y)
#include<stdio.h>
#include<string.h>
#define max 100001
__int64 Prim[max],sum[max],num[max][20];
void init()
{
    //打表求出1--100000中的所有数进行质因数分解
    memset(Prim,0,sizeof(Prim));
    memset(sum,0,sizeof(sum));
    for(int i=2; i<=100000; i++)
        if(Prim[i]==0)
        {
            num[i][0]=i;
            sum[i]++;
            for(int j=i*2; j<=100000; j+=i)
            {
                num[j][sum[j]]=i;
                sum[j]++;
                Prim[j]=1;
            }
        }
}
int Noprim(int m,int n,int k)
{
    int i;
    __int64 ans=0;
    if(m==1) return 0;
    for(i=k;i<sum[m];i++)
    ans+=n/num[m][i]-Noprim(m,n/num[m][i],i+1);
    return ans;
}
int main()
{
    int Stask;
    int m,n;
    __int64 ans;
    init();
    scanf("%d",&Stask);
    for(int i=0; i<Stask; i++)
    {
        scanf("%d%d",&m,&n);
        if(m>n)
        {
            n=m^n;
            m=m^n;
            n=m^n;
        }
       ans=n;
       for(int i=2;i<=m;i++)
         ans+=n-Noprim(i,n,0);
         printf("%I64d\n",ans);
    }
    return 0;
}

你可能感兴趣的:(http://acm.hdu.edu.cn/showproblem.php?pid=2841 题意很简单)