bzoj 2820-莫比乌斯函数反演+分块优化 +换元优化

BZOJ2820 YGY的GCD

•题目大意：求有多少数对(x,y)(1<=x<=n,1<=y<=m)满足gcd(x,y)为质数

T = 10000

N, M <= 10000000

参考popoqqqPPT：

首先，根据http://blog.csdn.net/viphong/article/details/52451540

 

f(x)表示gcd(i,j)==x的ij对数，F(x)是包含f(x)的

那么对于一个质数其贡献为 f(i)= 

直接计算，光枚举p就已经TLE了。

我们可以发现，p是所有的质数，f(p)枚举了每个质数的所有倍数，这其实不就是把所有的数都枚举了吗：

也就是从 

变成了

而我们发现，这个跟又跟那个分块优化的套路非常相似了

我们只需要预处理出 的前缀和，就可以sqrt(n)求到答案了

这个东西什么意思？ 就是对于一个d， 把mu(d的约数)求和。

而这个我们可以预处理，枚举所有质数，去更新这个质数的所有倍数。

假设我们枚举的是n个数，那么每一个数的倍数有n/i个，一共枚举的数量就是=n* ，调和级数嘛，也就是n*lgn咯，似乎可以证明是均摊logn，而质数的个数约为O（n/lgn），因此

枚举所有质数的倍数的复杂度约为，o(logn*n/logn)=o(n)

再求个前缀和，所以预处理这一部分的复杂度是o(n)

综上复杂度为sqrt(n)

参考代码：

#include<iostream>
#include<cstdio>
#include<algorithm>
using namespace std;
typedef long long  ll;
const ll p =1000000007;
const int  N=10000000;
bool is_prime[N+500];
int prime[N+50];
int mu[N+50];
ll sum[N+50];
ll tot;

void Moblus()
{
    tot = 0;
    mu[1] = 1;
    for(ll i = 2; i < N; i++)
    {
        if(!is_prime[i])
        {
            prime[tot++] = i;
            mu[i] = -1;
        }
        for(ll j = 0; j < tot && i*prime[j] < N; j++)
        {
            is_prime[i*prime[j]] = 1;
            if(i % prime[j])
            {
                mu[i*prime[j]] = -mu[i];
            }
            else
            {
                mu[i*prime[j]] = 0;
                break;
            }
        }
    }
}

ll ud[N+50];
void pre()
{
     for (int i=0;i<tot;i++)
    {
        for (ll j=prime[i];j<=N;j+=prime[i])
        {
            ud[j]+=mu[j/prime[i]];
        }
    }
}

//找[1,n],[1,m]内gcd(i,j)为质数的对数
ll solve(ll n,ll m )
{
    if (n>m)swap(n,m);
    ll ret=0;
    for (int i=1,last;i<=n;i=last+1)
    {
        last=min(n/(n/i),m/(m/i));
        ret+=(sum[last]-sum[i-1])*(n/i)*(m/i);
    }
    return ret;
}
int main()
{
    Moblus();
    pre();
    for (int i=1;i<N;i++)
        sum[i]=sum[i-1]+ud[i];
    freopen("YYnoGCD.in","r",stdin);
    freopen("YYnoGCD.out","w",stdout);
    int t;cin>>t;
    while(t--)
    {
        ll n,m;
        scanf("%lld%lld",&n,&m);
        ll ans=solve(n,m);

        printf("%lld\n",ans);
    }


    return 0;
}