bzoj2154: Crash的数字表格

最新推荐文章于 2017-06-18 21:29:31 发布

zhouyuyang233

最新推荐文章于 2017-06-18 21:29:31 发布

阅读量328

点赞数

分类专栏：数论

本文链接：https://blog.csdn.net/zhouyuyang233/article/details/70048677

版权

数论专栏收录该内容

49 篇文章 0 订阅

订阅专栏

传送门
首先 $lcm(i,j)=\dfrac{i*j}{gcd(i,j)}$
所以 $答案=\sum_{i=1}^{n}\sum_{j=1}^m \dfrac{i*j}{gcd(i,j)}$

设 $d=gcd(i,j)$
则

答 案 = \sum d = 1 \sum i = 1 n d \sum j = 1 m d i * j (g c d (i, j) = = 1)

$答案=\sum_{d=1}\sum_{i=1}^{\dfrac{n}{d}}\sum_{j=1}^{\dfrac{m}{d}} i*j(gcd(i,j)==1)$

= \sum d = 1 \sum i = 1 n d i \sum j = 1 m d j (g c d (i, j) = = 1)

$=\sum_{d=1}\sum_{i=1}^{\dfrac{n}{d}}i\sum_{j=1}^{\dfrac{m}{d}} j(gcd(i,j)==1)$
显然

(gcd(i,j)==1) $(gcd(i,j)==1)$ 是可以用莫比乌斯反演在sqrt(i,j)的时间内求出来的。
再对

d <script type="math/tex" id="MathJax-Element-35">d</script>进行分块就行了。

#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<iostream>
#include<algorithm>
#define ll long long
using namespace std;
const ll mo=20101009;
ll n,m,l,r,ans,mou[10000005];
ll sum(ll x,ll y){return (x*(x+1)/2%mo)*(y*(y+1)/2%mo)%mo;}
ll F(ll n,ll m){
    ll ans=0,l,r;
    for (l=1;l<=n;l=r+1){
        r=min(n/(n/l),m/(m/l));
        ans+=(mou[r]-mou[l-1])%mo*sum(n/l,m/l)%mo;
        ans%=mo;
    }
    return ans;
}
int main(){
    scanf("%lld%lld",&n,&m);
    if (n>m) swap(n,m);
    for (ll i=2;i<=m;i++)
        if (!mou[i]){
            mou[i]=i;
            for (ll j=i*i;j<=m;j+=i) mou[j]=i;
        }
    mou[1]=1;
    for (ll i=2;i<=m;i++)
        if (i/mou[i]%mou[i]==0) mou[i]=0;
        else mou[i]=-mou[i/mou[i]];
    for (ll i=1;i<=m;i++)
        mou[i]=(mou[i-1]+mou[i]*i%mo*i%mo)%mo;
    for (l=1;l<=n;l=r+1){
        r=min(n/(n/l),m/(m/l));
        ans+=(l+r)*(r-l+1)/2%mo*F(n/l,m/l)%mo;
        ans%=mo;
    }
    if (ans<0) ans+=mo;
    printf("%lld",ans);
}