BZOJ2820: YY的GCD

最新推荐文章于 2019-07-04 22:58:56 发布

CR1SceNT

最新推荐文章于 2019-07-04 22:58:56 发布

阅读量337

点赞数

分类专栏：莫比乌斯反演

本文链接：https://blog.csdn.net/Ep1C_HeReT1c/article/details/72743241

版权

莫比乌斯反演专栏收录该内容

9 篇文章 0 订阅

订阅专栏

BZOJ2820

直接枚举质数像1101那样做肯定会T…然而我还是不信邪的T了一发才开心。。

根据前面的题的经验。
这道题的式子可以先化成

\sum k = 1 p [k] < = m i n (n, m) \sum d = 1 m i n (n, m) μ (d) * ⌊ n p k * d ⌋ ⌊ m p k * d ⌋

$\sum_{k=1}^{p[k]<=min(n,m)}{\sum_{d=1}^{min(n,m)}}μ(d)*\lfloor{\frac{n}{{p_k*d}}}\rfloor\lfloor{\frac{m}{{p_k*d}}}\rfloor$

将式子进一步推演：
令 $T=p*d$

\sum T = 1 m i n (n, m) (⌊ n T ⌋ ⌊ m T ⌋ * \sum d | T μ (T d))

$\sum_{T=1}^{min(n,m)}(\lfloor{\frac{n}{{T}}}\rfloor\lfloor{\frac{m}{{T}}}\rfloor*\sum_{d|T}μ(\frac{T}{d}))$

预 处 理 出 \sum d | T μ (T d)

$预处理出\sum_{d|T}μ(\frac{T}{d})$ 直接枚举质数算就好了，更厉害一点可以线性筛求。

最后就分块求就好了。

【代码】


#include <cstdio>
#include <iostream>
#include <algorithm>
#define N 10000005
using namespace std;
typedef long long ll;
int read()
{
    int x=0,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-') f=-1;ch=getchar();}
    while(isdigit(ch)){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
    return x*f;
}

int T,n,m;
int Miu[N],p[664580],sum[N];
bool Not_Prime[N];

void Get_Miu()
{
    Miu[1]=1;
    for(register int i=2;i<N;i++) {
        if(!Not_Prime[i]) p[++p[0]]=i,Miu[i]=-1;
        for(register int j=1;j<=p[0]&&i*p[j]<N;j++)
        {
            Not_Prime[i*p[j]]=1;
            if(i%p[j]==0) {Miu[i*p[j]]=0;break;}
            Miu[i*p[j]]=-Miu[i];
        }
    }
    for(register int i=1;i<=p[0];i++)
    for(register int j=p[i];j<N;j+=p[i]) sum[j]+=Miu[j/p[i]];
    for(register int i=1;i<N;i++) sum[i]+=sum[i-1];
}

void Calc()
{
    ll ans=0;int pos=1;
    for(int i=1;i<=n;i=pos+1)
    {
        pos=min(n/(n/i),m/(m/i));
        ans+=1LL*(sum[pos]-sum[i-1])*(n/i)*(m/i);
    }
    printf("%lld\n",ans);
}

int main()
{
    Get_Miu();
    T=read();
    while(T--)
    {
        n=read(),m=read();
        if(n>m) swap(n,m);
        Calc();
    }
    return 0;
}