BZOJ2154 Crash的数字表格

原题链接：https://www.lydsy.com/JudgeOnline/problem.php?id=2154

Crash的数字表格

Description

今天的数学课上，Crash小朋友学习了最小公倍数(Least Common Multiple)。对于两个正整数a和b，LCM(a, b)表示能同时被a和b整除的最小正整数。例如，LCM(6, 8) = 24。回到家后，Crash还在想着课上学的东西，为了研究最小公倍数，他画了一张N*M的表格。每个格子里写了一个数字，其中第i行第j列的那个格子里写着数为LCM(i, j)。一个4*5的表格如下： 1 2 3 4 5 2 2 6 4 10 3 6 3 12 15 4 4 12 4 20 看着这个表格，Crash想到了很多可以思考的问题。不过他最想解决的问题却是一个十分简单的问题：这个表格中所有数的和是多少。当N和M很大时，Crash就束手无策了，因此他找到了聪明的你用程序帮他解决这个问题。由于最终结果可能会很大，Crash只想知道表格里所有数的和mod 20101009的值。

Input

输入的第一行包含两个正整数，分别表示N和M。

Output

输出一个正整数，表示表格中所有数的和mod 20101009的值。

Sample Input

4 5

Sample Output

122

【数据规模和约定】

100%的数据满足N, M ≤ 10^7。

题解

初始式子：

$$ans=\sum_{i=1}^n\sum_{j=1}^mlcm(i,j)$$

显然有$lcm(a,b)=\frac{ab}{gcd(a,b)}$：

$$\sum_{i=1}^n\sum_{j=1}^m\frac{ij}{gcd(i,j)}$$

基本操作之更该枚举项+提取公因数：

$$\begin{align*} \sum_{i=1}^n\sum_{j=1}^m\frac{ij}{gcd(i,j)}&=\sum_{d=1}^{min(n,m)}\sum_{i=1}^n\sum_{j=1}^m\frac{ij}{d}[d=gcd(i,j)]\\ &=\sum_{d=1}^{min(n,m)}d^2\sum_{i=1}^{\lfloor\frac{n}{d}\rfloor}\sum_{j=1}^{\lfloor\frac{m} {d}\rfloor}\frac{ij}{d}[1=gcd(i,j)]\\ &=\sum_{d=1}^{min(n,m)}d\sum_{i=1}^{\lfloor\frac{n}{d}\rfloor}\sum_{j=1}^{\lfloor\frac{m}{d}\rfloor}ij[1=gcd(i,j)] \end{align*}$$

我们把$\sum_{i=1}^{\lfloor\frac{n}{d}\rfloor}\sum_{j=1}^{\lfloor\frac{m}{d}\rfloor}ij[1=gcd(i,j)]$提出来单独调♂教，设$F(x,y)=\sum_{i=1}^{x}\sum_{j=1}^{y}ij[1=gcd(i,j)]$：

$$\begin{align*} F(x,y)&=\sum_{i=1}^{x}\sum_{j=1}^{y}ij[1=gcd(i,j)]\\ &=\sum_{i=1}^{x}\sum_{j=1}^{y}ij\sum_{d|gcd(i,j)}\mu(d)\\ &=\sum_{d=1}^{min(x,y)}\mu(d)\sum_{i=1}^{x}\sum_{j=1}^{y}ij[d|gcd(i,j)]\\ &=\sum_{d=1}^{min(x,y)}\mu(d)d^2\sum_{i=1}^{\lfloor\frac{x}{d}\rfloor}\sum_{j=1}^{\lfloor\frac{y}{d}\rfloor}ij\\ &=\sum_{d=1}^{min(x,y)}\mu(d)d^2\sum_{i=1}^{\lfloor\frac{x}{d}\rfloor}i\sum_{j=1}^{\lfloor\frac{y}{d}\rfloor}j\\ &=\sum_{d=1}^{min(x,y)}\mu(d)d^2\frac{(1+\lfloor\frac{x}{d}\rfloor)\lfloor\frac{x}{d}\rfloor}{2}\frac{(1+\lfloor\frac{y}{d}\rfloor)\lfloor\frac{y}{d}\rfloor}{2}\\ \end{align*}$$

代回原式：

$$\begin{align*} \sum_{d=1}^{min(n,m)}d\times F(\lfloor\frac{n}{d}\rfloor,\lfloor\frac{m}{d}\rfloor)&=\sum_{d=1}^{min(n,m)}d\sum_{d'=1}^{\lfloor\frac{min(n,m)}{d}\rfloor}\mu(d')d'^2\frac{(1+\lfloor\frac{x}{dd'}\rfloor)\lfloor\frac{x}{dd'}\rfloor}{2}\frac{(1+\lfloor\frac{y}{dd'}\rfloor)\lfloor\frac{y}{dd'}\rfloor}{2} \end{align*}$$

自然想到将$dd'$替换为$T$并枚举$T$：

$$=\sum_{T=1}^{min(n,m)}\frac{(1+\lfloor\frac{x}{T}\rfloor)\lfloor\frac{x}{T}\rfloor}{2}\frac{(1+\lfloor\frac{y}{T}\rfloor)\lfloor\frac{y}{T}\rfloor}{2}T\sum_{d'|T}\mu(d')d'$$

设$T=i\times p[j],f(T)=T\sum_{d'|T}\mu(d')d',g(T)=\sum_{d'|T}\mu(d')d'$，我们单独拿出$g(T)$讨论如何线筛：

1.当$T$为素数时，$d'$的取值只能为$1,T$，此时$g(T)=1-T$；

2.当$T$拥有多个最小质因子时，由于没有新增质因子，所以不影响我们的枚举，此时$g(T)=g(i)$；

3.当$T$的最小质因子只有一个时，在原来枚举的所有因数的基础上，我们还多了一堆“原来的因数$\times p[j]$”，因为多了一个质因数之后，$\mu(d'\times p[j])=-\mu(d')$，所以这个最小质因子的贡献为$-p[j]\times g(i)$，此时$g(T)=g(i)\times(1-p[j])$。

在筛完$g(T)$后，依次乘以$T$即可得到$f(T)$。

代码实现时，对$f(T)$求前缀和，套上下底分块即可实现$O(\sqrt{n})$求解。

代码

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int M=1e7+5,mod=20101009;
int n,m,p[M/3];
ll f[M],ans,s1,s2;
bool isp[M];
void pre()
{
    f[1]=1;int t;
    for(int i=2;i<=n;++i)
    {
        if(!isp[i])p[++p[0]]=i,f[i]=1-i;
        for(int j=1;j<=p[0];++j)
        {
            t=i*p[j];if(t>n)break;
            isp[t]=1;
            if(i%p[j]==0){f[t]=f[i];break;}
            f[t]=f[i]*(1-p[j]);
        }
        f[i]=(f[i-1]+f[i]*i%mod)%mod;
    }
}
void in(){scanf("%d%d",&n,&m);}
void ac()
{
    if(n>m)swap(n,m);pre();
    for(int l=1,r;l<=n;l=r+1)
    {r=min(n/(n/l),m/(m/l));s1=1ll*(1+n/l)*(n/l)/2%mod;s2=1ll*(1+m/l)*(m/l)/2%mod;(ans+=(f[r]-f[l-1])*((s1*s2)%mod)%mod)%=mod;}
    printf("%lld",(ans+mod)%mod);
}
int main(){in();ac();}