BZOJ 2693 jzptab

http://www.lydsy.com/JudgeOnline/problem.php?id=2693

题解：

考虑把lcm转化成gcd
那答案就是

然后神奇的设：

就有：

一样可以枚举

的取值，这是O(√n)的；

然后求f(x,y)；

大概证明了一下= =

线性筛之后也可以O(√n)求出f(x,y)
总复杂度O(n)，常数略大。。

这题显然是卡O(n)过不了呗
那就还得优化

预处理这玩意

然后O(√n)就搞出来啦！

设

“积性函数的约数和也是积性函数”  ->好像比较显然?
所以g(D)是积性函数
线性筛裸上就好

#include<cstdio>
#include<cstring>
#include<cmath>
#include<iostream>
#include<algorithm>
#define ll long long
const ll Mod=;
ll g[],p[],sum[];
bool mark[];
ll read(){
    ll t=,f=;char ch=getchar();
    while (ch<''||ch>''){if (ch=='-') f=-;ch=getchar();}
    while (''<=ch&&ch<=''){t=t*+ch-'';ch=getchar();}
    return t*f;
}
void init(){
    g[]=sum[]=;
    for (ll i=;i<;i++){
        if (!mark[i]){
            p[++p[]]=i;
            g[i]=(ll)((i-i*i)%Mod+Mod)%Mod;
        }
        for (ll j=;j<=p[]&&p[j]*i<=;j++){
            mark[i*p[j]]=;
            if (i%p[j]==){
                g[i*p[j]]=g[i]*(p[j])%Mod;
                break;
            }else
            g[i*p[j]]=g[i]*g[p[j]]%Mod;
        }
        sum[i]=sum[i-]+g[i];
    }
}
ll F(ll x,ll y){
    return (((x*(x+)/2LL)%Mod)*((y*(y+)/2LL)%Mod))%Mod;
}
int main(){
    init();int T=read();
    while (T--){
        ll n=read(),m=read();
        if (n>m) std::swap(n,m);
        ll j;ll ans=;
        for (ll i=;i<=n;i=j+){
             j=std::min(n/(n/i),m/(m/i));
             ans+=((sum[j]-sum[i-]%Mod+Mod)%Mod)*F(n/i,m/i);
             ans%=Mod;
        }
        printf("%lld\n",ans);
    }
}