Gcd反应堆 (pgcd)

Gcd反应堆 (pgcd)

题目描述

 

不知什么时候起,TA突然对gcd产生了浓厚的兴趣,于是他为此编写了个程序,输入分别不大于m,n (1<m,n<=10^7)的两个数,就能得出gcd(m,n)!TA兴奋之余,给这个程序起了个酷(zhuang)雅(bi)的名字:gcd反应堆。

现在TA要做一项伟大的科学实验,需要大量素数,他打算生成素数这种任务就交给他的反应堆实现。现在TA问你,多少对在他反应堆承载范围内的整数对i,j (1<i<=m,1<=j<=n),生成的数是他所需要的原料呢?

 

 

输入

 

本题包含多组数据,第一行包含一个整数t,表示数据组数。

接下来T行,每行两个整数m,n,意义为题目所述。

 

 

输出

 

输出T行,每行一个整数,为满足(1<i<=m,1<=j<=n)且gcd(i,j)为素数的整数对(i,j)的个数。

 

【数据规模】

对于30%的数据,m,n<=10^3

对于70%的数据,m,n<=10^5

对于100%的数据额,m,n<=10^7,T<=10


solution

题目求

\sum_{i=1}^{n}\sum_{j=1}^{m} [ gcd(i,j)==prime]

先枚举质数p

\sum_{p} \sum_{i=1}^{n}\sum_{j=1}^{m} [ gcd(i,j)==p]

除掉p

\sum_{p} \sum_{i=1}^{n/p}\sum_{j=1}^{m/p} [ gcd(i,j)==1]

套用反演的式子

\sum_{p} \sum_{i=1}^{n/p}\sum_{j=1}^{m/p} \sum_{d|gcd(i,j)} \mu (d)

先枚举gcd

\sum_{p} \sum_{d} \mu (d) \sum_{i=1}^{n/pd}\sum_{j=1}^{m/pd}

后面可以O(1)算,前面O(n)预处理mu

效率O(nlogn)

#include<cstdio>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
int T,n,m,mu[10000007],pri[10000007],tot;
bool flag[10000007];
int main()
{
    int Max=10000000;mu[1]=1;
    for(int i=2;i<=Max;i++){
        if(!flag[i]){mu[i]=-1;pri[++tot]=i;}
        for(int j=1;j<=tot&&i<=Max/pri[j];j++){
            flag[i*pri[j]]=1;mu[i*pri[j]]=-mu[i];
            if(i%pri[j]==0){mu[i*pri[j]]=0;break;}
        }
    }
    cin>>T;
    while(T--){
        scanf("%d%d",&n,&m);
        if(n<m)swap(n,m);
        long long ans=0;
        for(int i=1;i<=tot&&pri[i]<=m;i++){
            int p=pri[i];
            for(int d=1;p<=m/d;d++){
                ans=ans+(1LL)*mu[d]*(n/(p*d))*(m/(p*d));
            }
        }
        cout<<ans<<endl;
    }
    return 0;
}
 

 

 

 

posted @ 2018-11-30 19:21  liankewei123456  阅读(449)  评论(0编辑  收藏  举报