4407: 于神之怒加强版 莫比乌斯反演+线性筛积性函数

好虚。。不会做啊
Ans=ni=1mj=1gcd(i,j)k ,我们考虑枚举公约数 d ,然后令f(d)表示 gcd(i,j)=d 的数对 (i,j) 的数量,然后莫比乌斯反演可得 f(d)=ndx=1μ(x)ndxmdx ,所以有

Ans=d=1ndk×x=1ndμ(x)ndxmdx

我们令 T=dx ,
Ans=T=1nnTmTdTdkμ(Td)

g(T)=dTdkμ(Td) ,求出 g(T) 就可以分块辣。
可以发现 g(T) 显然是一个积性函数,然后我们可以进行线性筛。
is_prime(x) ,则 g(x)=xk1
否则,若 gcd(x,p)=1 ,则 g(x×p)=g(x)×g(p)
否则 g(x×p)=g(x)×pk
前两个是显然的,我们来YY一下第三个。
我们加入一个已经存在的质数 p ,由于含pi(i>1)的数的莫比乌斯函数为0,所以加入一个数后式子的项数是不变的,而每个数都多乘了一个 pk ,所以将原来的和乘 pk 就好辣。

#include<iostream>
#include<cstdio>
#define P 1000000007
#define ll long long
#define MAXN 5000005
using namespace std;
int prime[MAXN];
ll f[MAXN],p[MAXN];
bool flag[MAXN];
int k;
inline int read()
{
    int a=0,f=1; char c=getchar();
    while (c<'0'||c>'9') {if (c=='-') f=-1; c=getchar();}
    while (c>='0'&&c<='9') {a=a*10+c-'0'; c=getchar();}
    return a*f;
}
ll quick_power(ll a,ll b)
{
    ll ans=1;
    for (ll i=b;i;i>>=1,a=a*a%P)
        if (i&1) ans=ans*a%P;
    return ans;
}
inline void prepare()
{
    f[1]=1;
    for (int i=2;i<MAXN;i++)
    {
        if (!flag[i]) prime[++prime[0]]=i,p[prime[0]]=quick_power(i,k),f[i]=p[prime[0]]-1;
        for (int j=1;j<=prime[0]&&i*prime[j]<MAXN;j++)
        {
            flag[i*prime[j]]=1;
            if (i%prime[j]==0)
            {
                f[i*prime[j]]=f[i]*p[j]%P;
                break;
            }
            else f[i*prime[j]]=f[i]*f[prime[j]]%P;
        }
    }

    for (int i=1;i<MAXN;i++) (f[i]+=f[i-1])%=P;
}
int main()
{
    int testcase=read(); k=read();
    prepare();
    while (testcase--)
    {
        ll ans=0;
        int n=read(),m=read();
        if (n>m) swap(n,m);
        for (int i=1,pos;i<=n;i=pos+1)
        {
            pos=min(n/(n/i),m/(m/i));
            (ans+=(f[pos]-f[i-1]+P)%P*(n/i)%P*(m/i)%P)%=P;
        }
        printf("%lld\n",ans);
    }
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值