BZOJ 4407: 于神之怒加强版|莫比乌斯反演

不会搞数学公式很苦恼！！
flag：会写数学公式之后一定好好写一发题解
非常感谢龙爷（sd第一男选手！！可惜神犇都不写blog）提供线性筛做法
2.16————————————————-
一下均设 n<=m

Ans=∑i=1n∑j=1mgcd(i,j)k

设 f(d) 为 gcd(x,y)=d 的 (x,y) 数对的数量

g(d)=∑i=1⌊nd⌋f(i∗d)=⌊nd⌋∗⌊md⌋

f(d)=∑i=1⌊nd⌋u(i)∗g(i∗d)=∑i=1⌊nd⌋u(i)∗⌊md∗i⌋∗⌊nd∗i⌋

Ans=∑d=1ndk∑i=1⌊nd⌋u(i)∗⌊md∗i⌋∗⌊nd∗i⌋

此时两次分块可以做到单次讯问 O(n) 但是还不够
继续等价变换令 T=d∗i

Ans=∑T=1n⌊nT⌋∗⌊mT⌋∑d|Tdk∗u(Td)

h(T)=∑d|Tdk∗u(Td)

只需要搞出 h(T) 的前缀和再用分块处理就可以做到单次讯问 n−√ 的复杂度
然后就是鬼畜的线性筛，还要感谢龙爷

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define N 5000005
#define ll long long
#define R 1000000007ll
using namespace std;
int sc()
{
    int i=0,f=1; char c=getchar();
    while(c>'9'||c<'0'){if(f=='-')f=-1;c=getchar();}
    while(c>='0'&&c<='9')i=i*10+c-'0',c=getchar();
    return i;
}
ll f[N],K;
int prime[N],low[N];
int a[N],n;
ll cal(ll x)
{
    ll ans=1,y=K;
    while(y)
    {
        if(y&1)ans=ans*x%R;
        x=x*x%R;
        y>>=1;
    }
    return (ans-1+R)%R;
}
void pre()
{
    int top=0; f[1]=low[1]=1;
    for(int i=2;iif(!a[i])
        {
            low[i]=prime[++top]=i;
            f[i]=cal(i);
        }
        for(int j=1;i*prime[j]*prime[j]]=1;
            if(i%prime[j]==0)
            {
                low[i*prime[j]]=low[i]*prime[j];
                if(low[i]==i) 
                    f[i*prime[j]]=f[i]*(f[prime[j]]+1)%R;
                else
                    f[i*prime[j]]=f[i/low[i]]*f[low[i]*prime[j]]%R;
                break;
            }
            low[i*prime[j]]=prime[j];
            f[i*prime[j]]=f[i]*f[prime[j]]%R;
        }
    }
    for(int i=2;i1])%R;
}
int main()
{
    n=sc(),K=sc();
    pre();
    while(n--)
    {
        ll ans=0;
        int n=sc(),m=sc();
        if(n>m)swap(n,m);
        for(int i=1,last;i<=n;i=last+1)
        {
            last=min(n/(n/i),m/(m/i));
            ans=(ans+((ll)(n/i)*(m/i)%R)*(R+f[last]-f[i-1])%R)%R;
        }
        printf("%lld\n",ans);
    }
    return 0;
}