bzoj 3561: DZY Loves Math VI

最新推荐文章于 2022-02-12 10:50:10 发布

OI界第一麻瓜

最新推荐文章于 2022-02-12 10:50:10 发布

阅读量245

点赞数

分类专栏：数论

本文链接：https://blog.csdn.net/qq_36797743/article/details/79690765

版权

数论专栏收录该内容

23 篇文章 0 订阅

订阅专栏

题意

bzoj3561

题解

今天听课的时候听到一个一个gcd在下面的版本。。
然后积分算复杂度，好吧，我还不会积分
于是来做一下这题
推了一下式子，发现很简单啊
m开500000，那我 $mlogm$ 不就好了。。
调和级数啊！
然后就没什么了。。
式子很简单的，自己推一下吧。。
懒）

CODE:

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
using namespace std;
typedef long long LL;
const int MOD=1e9+7;
const int N=500005;
int n,m;
int mu[N],pri[N],tot;
bool ok[N];
void get_mu ()
{
    memset(ok,true,sizeof(ok));
    mu[1]=1;
    for (int u=2;u<=500000;u++)
    {
        if (ok[u])  {pri[++tot]=u;mu[u]=-1;}
        for (int i=1;i<=tot;i++)
        {
            int j=pri[i];
            if (u*j>500000) break;
            ok[u*j]=false;
            if (u%j!=0) mu[u*j]=mu[u]*mu[j];
            else {mu[u*j]=0;break;}
        }
    }
}
LL pow (LL x,int y)
{
    if (y==1) return x;
    LL lalal=pow(x,y>>1);
    lalal=lalal*lalal%MOD;
    if (y&1) lalal=lalal*x%MOD;
    return lalal;
}
int t[N];
LL o[N];
int main()
{
    get_mu();
    scanf("%d%d",&n,&m);
    if (n>m) swap(n,m);
    LL ans=0;
    for (int u=1;u<=m;u++) o[u]=1;
    for (int d=1;d<=n;d++)
    {
        for (int u=1;u<=m/d;u++) 
        {
            o[u]=o[u]*u%MOD;
            t[u]=(t[u-1]+o[u])%MOD;
        }
        LL lalal=pow(d,d);
        int sum=0;
        for (int u=1;u<=n/d;u++)    
            sum=(sum+mu[u]*pow(u,2*d)%MOD*t[n/d/u]%MOD*t[m/d/u]%MOD)%MOD;
        ans=(ans+sum*lalal%MOD)%MOD;
    }
    ans=(ans+MOD)%MOD;
    printf("%lld\n",ans);
    return 0;
}