bzoj3994 [SDOI2015]约数个数和

题目

反演的经典套路题目。

对于这种题，套路无非就是先用狄利克雷卷积一下，然后交换求和顺序，最后就能化为一个 $\sqrt n$的求和了。

对于交换顺序这个事，只要注意一下各个变量的依赖关系就好了。

对于这道题来说：

$\sum_{i=1}^n \sum_{j=1}^md(i,j)=\sum_{i=1}^n \sum_{j=1}^m [\frac{n}{i}][\frac{m}{j}][gcd(i,j)==1]$
$\sum_{i=1}^n \sum_{j=1}^m [\frac{n}{i}][\frac{m}{j}][gcd(i,j)==1]=\sum_{i=1}^n\sum_{j=1}^m\sum_{d|gcd(i,j)}[\frac{n}{i}][\frac{m}{j}]\mu(d)$
$\sum_{i=1}^n\sum_{j=1}^m\sum_{d|gcd(i,j)}[\frac{n}{i}][\frac{m}{j}]\mu(d)=\sum_{d=1}^{min(m,n)}\sum_{i=1}^{[\frac{n}{d}]}\sum_{j=1}^{[\frac{m}{d}]}\mu(d)[\frac{n}{i}][\frac{m}{j}]$
$\sum_{d=1}^{min(m,n)}\sum_{i=1}^{[\frac{n}{d}]}\sum_{j=1}^{[\frac{m}{d}]}\mu(d)[\frac{n}{i}][\frac{m}{j}]=\sum_{d=1}^{min(m,n)}\mu(d)\sum_{i=1}^{[\frac{n}{d}]}[\frac{n}{i}]\sum_{j=1}^{[\frac{m}{d}]}[\frac{m}{j}]$
$\sum_{d=1}^{min(m,n)}\mu(d)\sum_{i=1}^{[\frac{n}{d}]}[\frac{n}{i}]\sum_{j=1}^{[\frac{m}{d}]}[\frac{m}{j}]=\sum_{d=1}^{min(m,n)}\mu[d]f([\frac{n}{d}])f([\frac{m}{d}])$

所以，最终答案就是$\sum_{d=1}^{min(m,n)}\mu[d]f([\frac{n}{d}])f([\frac{m}{d}])$,其中，f(n)表示n的约数个数。

之后就又是套路的求前缀和了。

#include<bits/stdc++.h> 
#define N 50000
using namespace std;
int miu[N+5],f[N+5],d[N+5];
int prime[N+5],P[N+5],cnt;
int T,m,n;
long long ans;
inline char nc()
{
    static char buf[100000],*p1=buf,*p2=buf;
    return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}
inline int read()
{
    int x=0,b=1;
    char c=nc();
    for(;!(c<='9'&&c>='0');c=nc())if(c=='-')b=-1;
    for(;c<='9'&&c>='0';c=nc())x=x*10+c-'0';
    return x*b;
}
inline void init()
{
    miu[1]=f[1]=1;P[1]=true;
    for(int i=1;i<=N;i++)
    {
        if(!P[i])prime[++cnt]=i,f[i]=2,miu[i]=-1,d[i]=1;
        for(int j=1;j<=cnt&&i*prime[j]<=N;j++)
        {
            P[i*prime[j]]=true;
            if(i%prime[j]==0)
            {
                miu[i*prime[j]]=0;
                d[i*prime[j]]=d[i]+1;
                f[i*prime[j]]=f[i]/(d[i]+1)*(d[i]+2);
                break;
            }
            miu[i*prime[j]]=-miu[i];
            f[i*prime[j]]=2*f[i];
            d[i*prime[j]]=1;
        }
    }
    for(int i=1;i<=N;i++)f[i]+=f[i-1],miu[i]+=miu[i-1];
}
inline void write(long long x)
{
    if(x==0)putchar('0');
    else
    {
        char buf[15];
        int len=0;
        if(x<0)putchar('-'),x=-x;
        while(x)buf[++len]=x%10+'0',x/=10;
        for(int i=len;i>=1;i--)putchar(buf[i]);
    }
    putchar('\n');
}
int main()
{ 
    freopen("in.txt","r",stdin);
    init();
    T=read();
    while(T--)
    {
        n=read(),m=read();
        if(n>m)swap(m,n);
        ans=0;
        for(int i=1,j;i<=n;i=j+1)
        {
            j=min(n/(n/i),m/(m/i));
            ans+=1LL*(miu[j]-miu[i-1])*f[n/i]*f[m/i];
        }
        write(ans);
    }
    return 0;
}

哇，这个公式把我恶心到了，太难写了。。。