bzoj2154: Crash的数字表格

最新推荐文章于 2021-04-10 12:46:11 发布

Miao_zc

最新推荐文章于 2021-04-10 12:46:11 发布

阅读量643

点赞数 1

分类专栏： bzoj 数论文章标签： bzoj 莫比乌斯反演

本文链接：https://blog.csdn.net/Miao_zc/article/details/50899408

版权

bzoj 同时被 2 个专栏收录

69 篇文章 0 订阅

订阅专栏

数论

9 篇文章 0 订阅

订阅专栏

依题意： $d=gcd(i,j)$

$Ans=\sum_{i=1}^{n}\sum_{j=1}^{m}Lcm(i,j)=\sum_{i=1}^{n}\sum_{j=1}^{m}\frac{i}{d}\frac{j}{d}d$

$\Rightarrow Ans=\sum_{d=1}^{min(n,m)}d\sum_{i=1}^{\left \lfloor \frac{n}{d} \right \rfloor} \sum_{j=1}^{\left \lfloor \frac{m}{d} \right \rfloor}ij|gcd(i,j)=1$

$\Rightarrow Ans=\sum_{d=1}^{min(n,m)}d\sum_{i=1}^{\left \lfloor \frac{n}{d} \right \rfloor} \sum_{j=1}^{\left \lfloor \frac{m}{d} \right \rfloor}ij*\sum_{p|i,p|j}\mu(p)$

$\Rightarrow Ans=\sum_{d=1}^{min(n,m)}d\sum_{p=1}^{\left \lfloor min(\frac{n}{d},\frac{m}{d})\right \rfloor}\mu(p)p^{2}\sum_{i=1}^{\left \lfloor \frac{n}{dp} \right \rfloor} \sum_{j=1}^{\left \lfloor \frac{m}{dp} \right \rfloor}ij$

枚举d，对一段d有 $\left \lfloor \frac{n}{d} \right \rfloor 和\left \lfloor \frac{m}{d} \right \rfloor$ 不变， $O(\sqrt n)$ 。同理枚举p，总复杂度 $O(n)$ 。

感谢lych

#include<iostream>
#include<cstdio>
#define N 10000005
#define mo 20101009
#define ll long long
using namespace std;
ll n,m,Ans,mu[N],b[N],f[N];
void getMob()
{
    mu[1]=1;
    ll l=0;
    for (ll i=2;i<=m;i++)
    {
        if (f[i]==0) mu[i]=-1,b[++l]=i;
        for (ll j=1;j<=l;j++)
        {
            ll t=i*b[j];
            if (t>m)break;
            f[t]=1;
            if (i%b[j]==0) {mu[t]=0;break;}
            mu[t]=-mu[i];
        }
    }
    for (ll i=1;i<=m;i++)
        mu[i]=(mu[i-1]+mu[i]*i*i)%mo;
}
ll Get(ll a,ll b)
{
    a=a*(1+a)/2%mo;
    b=b*(1+b)/2%mo;
    return a*b%mo;
}
int main()
{
    scanf("%lld%lld",&n,&m);
    if (n>m) swap(n,m);
    getMob();
    for (ll d=1,lastd;d<=n;d=lastd+1)
    {
        lastd=min(n/(n/d),m/(m/d));
        ll sum=0,m1=m/d,n1=n/d;
        for (ll p=1,lastp;p<=n1;p=lastp+1)
        {
            lastp=min(n1/(n1/p),m1/(m1/p));
            sum+=(mu[lastp]-mu[p-1])*Get(n1/p,m1/p)%mo;
            sum=sum%mo;
        }
        Ans=(Ans+((d+lastd)*(lastd-d+1)/2)%mo*sum%mo)%mo;
    }
    printf("%lld\n",(Ans%mo+mo)%mo);
    return 0; 
}