[SDOI2015]约数个数和

题目链接：Click here

Solution:

首先，我们转化式子

[sum_{i=1}^nsum_{j=1}^m d(ij)\ sum_{i=1}^nsum_{j=1}^m sum_{x|i} sum_{y|j}[gcd(x,y)=1]\ ]

我们把(x,y)给提前

[sum_{x=1}^nsum_{y=1}^m lfloor frac{n}{x} floor lfloor frac{m}{y} floor [gcd(x,y)=1] ]

我们把(gcd(x,y))提前，(x,y)不太好看，再给他换个名字(i,j)

[sum_{i=1}^n sum_{j=1}^m lfloor frac{n}{i} floor lfloor frac{m}{j} floor sum_{d|i} sum_{d|j} mu(d)\ ]

我们把(d)提前

[sum_{d=1}^n mu(d) sum_{i=1}^{lfloor frac{n}{d} floor} sum_{j=1}^{lfloor frac{m}{d} floor} lfloor frac{n}{di} floor lfloor frac{m}{dj} floor ]

我们设一个函数(g(n)=sum_{i=1}^n lfloor frac{n}{i} floor)，则有

[sum_{d=1}^n mu(d) g(lfloor frac{n}{d} floor) g(lfloor frac{m}{d} floor) ]

这个数论分块即可，考虑怎么筛(g(n))，我们可以发现(g(n)=sum_{i=1}^n d(i))，则我们筛(d)后做前缀和即可

Code:

#include<bits/stdc++.h>
#define int long long
using namespace std;
const int N=5e4+11;
int n,m,u[N],p[N],vis[N];
int ans,cnt,g[N],num[N];
int read(){
	int x=0,f=1;char ch=getchar();
	while(!isdigit(ch)){if(ch=='-')f=-f;ch=getchar();}
	while(isdigit(ch)){x=x*10+ch-48;ch=getchar();}
	return x*f;
}
void prepare(){
    u[1]=1;g[1]=1;
    for(int i=2;i<N;i++){
        if(!vis[i]) p[++cnt]=i,u[i]=-1,g[i]=2,num[i]=1;
        for(int j=1;j<=cnt&&i*p[j]<N;j++){
            vis[i*p[j]]=1;
            if(i%p[j]==0){
                u[i*p[j]]=0;
                g[i*p[j]]=g[i]/(num[i]+1)*(num[i]+2);
                num[i*p[j]]=num[i]+1;
                break;
            }
            u[i*p[j]]=-u[i];
            g[i*p[j]]=g[i]*2;
            num[i*p[j]]=1;
        }
    }
    for(int i=1;i<N;i++)
        u[i]=u[i]+u[i-1],g[i]=g[i]+g[i-1]; 
}
void solve(){
    n=read(),m=read();
    ans=0;
    for(int i=1,j;i<=min(n,m);i=j+1){
        j=min(n/(n/i),m/(m/i));
        ans=ans+(u[j]-u[i-1])*g[n/i]*g[m/i];
    }
    printf("%lld
",ans);
}
signed main(){
    prepare();
    int t=read();
    while(t--) solve();
    return 0;
}