bzoj2154: Crash的数字表格

链接

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

题解

　　这道题当时觉得太难了不敢做，今天在$\_\_stdcall$大爷的帮助下学会了。
　　其实我就是不懂套路，莫比乌斯题的套路有两种，一种是用$\sum_{d|x}\mu(d)$，即莫比乌斯函数的性质。还有一种就是这道题的莫比乌斯反演。
　　下面推式子。
　　

$$ans=\sum_{i=1}^{n}\sum_{j=1}^mlcm(i,j)\\ 　　=\sum_{d=1}^n\frac{1}{d}\sum_{d|i}^n\sum_{d|j}^m[gcd(i,j)==d]$$
　　注意这里千万别化成 $[gcd(i,j)==1]$的形式，否则后面就化不动了。
　　令 $$f(n,m,d)=\sum_{d|i}^n\sum_{d|j}^m[gcd(i,j)==d] 　　\\s(n,m,d)=\sum_{d|i}^n\sum_{d|j}^mij 　　\\=\sum_{i=1}^{\lfloor\frac{n}{d}\rfloor}\sum_{j=1}^{\lfloor\frac{m}{d}\rfloor}idjd\\ 　　=d^2\sum_{i=1}^{\lfloor\frac{n}{d}\rfloor}\sum_{j=1}^{\lfloor\frac{m}{d}\rfloor}ij\\ 　　=d^2\frac{\lfloor\frac{n}{d}\rfloor(\lfloor\frac{n}{d}\rfloor+1)}{2}\frac{\lfloor\frac{m}{d}\rfloor(\lfloor\frac{m}{d}\rfloor+1)}{2}$$
　　显然
　　 $$s(n,m,d)=\sum_{d|x}^nf(n,m,x)$$
　　莫比乌斯反演
　　 $$f(n,m,d)=\sum_{d|x}^n\mu(\frac{x}{d})s(n,m,x)\\ 　　=\sum_{i=1}^{\lfloor\frac{n}{d}\rfloor}\mu(i)s(n,m,id)\\ 　　=\sum_{i=1}^{\lfloor\frac{n}{d}\rfloor}\mu(i)i^2d^2\frac{\lfloor\frac{n}{id}\rfloor(\lfloor\frac{n}{id}\rfloor+1)}{2}\frac{\lfloor\frac{m}{id}\rfloor(\lfloor\frac{m}{id}\rfloor+1)}{2}$$
　　那么
　　 $$ans=\sum_{d=1}^n\frac{1}{d}f(n,m,d)\\ 　　=\sum_{d=1}^n\frac{1}{d}\sum_{i=1}^{\lfloor\frac{n}{d}\rfloor}\mu(i)i^2d^2\frac{\lfloor\frac{n}{id}\rfloor(\lfloor\frac{n}{id}\rfloor+1)}{2}\frac{\lfloor\frac{m}{id}\rfloor(\lfloor\frac{m}{id}\rfloor+1)}{2}\\ 　　=\sum_{d=1}^nd\sum_{i=1}^{\lfloor\frac{n}{d}\rfloor}\mu(i)i^2\frac{\lfloor\frac{n}{id}\rfloor(\lfloor\frac{n}{id}\rfloor+1)}{2}\frac{\lfloor\frac{m}{id}\rfloor(\lfloor\frac{m}{id}\rfloor+1)}{2}$$
　　显然 $\mu(i)i^2$可以前缀和预处理。最后的时间复杂度是 $O(n)$(这个是我找 $will7101$用微积分算出来的)

代码

//莫比乌斯反演
#include <cstdio>
#include <algorithm>
#define maxn 10000001
#define mod 20101009
#define ll long long
using namespace std;
int mu[maxn], prime[maxn/10], N, M, s[maxn];
bool mark[maxn];
ll f(ll x){return x*(x+1)/2%mod;}
void shai()
{
    int i, j;
    mu[1]=1;
    for(i=2;i<=N;i++)
    {
        if(!mark[i])prime[++prime[0]]=i,mu[i]=-1;
        for(j=1;j<=prime[0] and i*prime[j]<=N;j++)
        {
            mark[i*prime[j]]=1;
            if(i%prime[j]==0){mu[i*prime[j]]=0;break;}
            mu[i*prime[j]]=-mu[i];
        }
    }
    for(i=1;i<=N;i++)s[i]=((ll)s[i-1]+(ll)mu[i]*i%mod*i)%mod;
}
void work()
{
    ll ans=0, t;
    int d, i, last, n, m, dd;
    for(d=1;d<=N;d=dd+1)
    {
        dd=min(N/(N/d),M/(M/d));
        n=N/d, m=M/d;
        for(i=1,t=0;i<=n;i=last+1)
        {
            last=min(n/(n/i),m/(m/i));
            t=(t+(s[last]-s[i-1])*f(n/i)%mod*f(m/i))%mod;
        }
        ans=(ans+(f(dd)-f(d-1))*t)%mod;
    }
    printf("%lld\n",(ans+mod)%mod);
}
int main()
{
    scanf("%d%d",&N,&M);
    if(N>M)swap(N,M);
    shai();
    work();
    return 0;
}