BZOJ2154 Crash的数字表格及其拓展BZOJ2693 jzptab（Mobius反演）

题目：BZOJ2154.
题目大意：给定$n$和$m$，求：
$$\sum _{i=1}^n\sum _{j=1}^{m}lcm(i,j)$$

$n,m\le 10^7$.

很容易发现：
$$\sum _{i=1}^n\sum _{j=1}^{m}lcm(i,j)=\sum _{i=1}^n\sum _{j=1}^m\frac{ij}{gcd(i,j)}$$

我们钦定$n\le m$，那么：
$$\begin{gathered} \sum _{i=1}^n\sum _{j=1}^{m}lcm(i,j) \\ =\sum _{k=1}^n\frac{1}{k}\sum _{i=1}^n\sum _{j=1}^{m}ij[gcd(i,j)=k] \end{gathered}$$

设$a=\lfloor \frac{n}{k}\rfloor ,b=\lfloor \frac{m}{k}\rfloor$，那么：
$$\begin{gathered} \sum _{i=1}^n\sum _{j=1}^{m}lcm(i,j) \\ =\sum _{k=1}^n\frac{1}{k}\sum _{i=1}^a\sum _{j=1}^{b}ij{k}^2\ast [gcd(i,j)=1] \\ =\sum _{k=1}^{n}k\sum _{i=1}^a\sum _{j=1}^{b}ij\sum _{d|gcd(i,j)}\mu (d) \\ =\sum _{k=1}^{n}k\sum _{d=1}^a\mu (d)\sum _{i=1}^a\sum _{j=1}^{b}ij[d|i][d|j] \end{gathered}$$

设$x=\lfloor \frac{a}{d}\rfloor ,y=\lfloor \frac{b}{d}\rfloor$，那么：
$$\begin{gathered} \sum _{i=1}^n\sum _{j=1}^{m}lcm(i,j) \\ =\sum _{k=1}^{n}k\sum _{d=1}^a{d}^2\mu (d)\sum _{i=1}^x\sum _{j=1}^{y}ij \\ =\sum _{k=1}^{n}k\sum _{d=1}^a{d}^2\mu (d)\ast \frac{x(x+1)}{2}\ast \frac{y(y+1)}{2} \end{gathered}$$

设函数$f(x,y)=\frac{x(x+1)}{2}\ast \frac{y(y+1)}{2}$，那么原式就变为：
$$\begin{gathered} \sum _{i=1}^n\sum _{j=1}^{m}lcm(i,j) \\ =\sum _{k=1}^{n}k\sum _{d=1}^a{d}^2\mu (d)f(\lfloor \frac{a}{d}\rfloor ,\lfloor \frac{b}{d}\rfloor ) \end{gathered}$$

设函数（钦定$x\le y$）$g(x,y)={\sum }_{d=1}^x{d}^2\mu (d)f(\lfloor \frac{x}{d}\rfloor ,\lfloor \frac{y}{d}\rfloor )$，那么原式变为：
$$\sum _{i=1}^n\sum _{j=1}^{m}lcm(i,j)=\sum _{k=1}^{n}kg(\lfloor \frac{n}{k}\rfloor ,\lfloor \frac{m}{k}\rfloor )$$

发现这个式子$g$函数的计算可以用除法分块优化，而外面式子也可以用除法分块优化，计算这个式子的时间复杂度即为$O((\sqrt{n}+\sqrt{m}{)}^2)=O(n+m)$.

代码如下：

#include<bits/stdc++.h>
  using namespace std;
 
#define Abigail inline void
typedef long long LL;
 
const int N=10000000;
const LL mod=20101009;
 
int pr[N+9],tp,b[N+9];
LL mu[N+9],div2;
 
LL M(LL x){
  x=x%mod+mod;
  return (x>=mod)?x-=mod:x; 
}
 
void sieve(int n){
  for (int i=2;i<=n;++i) b[i]=1;
  mu[1]=1;
  for (int i=2;i<=n;++i){
    if (b[i]) pr[++tp]=i,mu[i]=-1;
    for (int j=1;j<=tp&&pr[j]*i<=n;++j){
      b[i*pr[j]]=0;
      if (i%pr[j]) mu[i*pr[j]]=-mu[i];
      else {mu[i*pr[j]]=0;break;}
    }
  }
  for (int i=1;i<=n;++i)
    mu[i]=M(mu[i]*i*i+mu[i-1]);
}
 
LL sum(int l,int r){return M(LL(l+r)*LL(r-l+1)>>1);}
LL f(int x,int y){return M(sum(1,x)*sum(1,y));}
 
LL g(int n,int m){
  LL ans=0;
  if (n>m) swap(n,m);
  for (int l=1,r;l<=n;l=r+1){
    r=min(n/(n/l),m/(m/l));
    ans=M(ans+(mu[r]-mu[l-1])*f(n/l,m/l));
  }
  return ans;
}
 
LL Query(int n,int m){
  LL ans=0;
  if (n>m) swap(n,m);
  for (int l=1,r;l<=n;l=r+1){
    r=min(n/(n/l),m/(m/l));
    ans=M(ans+sum(l,r)*g(n/l,m/l));
  }
  return ans;
}
 
Abigail work(){
  sieve(N);
}
 
Abigail getans(){
  int n,m;
  scanf("%d%d",&n,&m);
  printf("%lld\n",Query(n,m));
}
 
int main(){
  work();
  getans();
  return 0;
}

我们再来看一道题，是上面这道题多组询问的形式.

题目：BZOJ2693.
题目大意：多组询问，给定$n$和$m$，求：
$$\sum _{i=1}^n\sum _{j=1}^{m}lcm(i,j)$$

多组询问了，可是$n$和$m$的数据范围还没变，怎么办呢？

我们考虑这样的问题应该是要预处理然后在低于$O(n)$的时间复杂度内回答询问.

观察下式：
$$\begin{gathered} \sum _{k=1}^{n}k\sum _{d=1}^a\mu (d)d^2\ast \frac{x(x+1)}{2}\ast \frac{y(y+1)}{x} \\ =\sum _{k=1}^{n}k\sum _{d=1}^a\mu (d)d^2\ast \frac{\lfloor \frac{n}{dk}\rfloor (\lfloor \frac{n}{dk}\rfloor +1)}{2}\ast \frac{\lfloor \frac{m}{dk}\rfloor (\lfloor \frac{m}{dk}\rfloor +1)}{2} \end{gathered}$$

我们让$t=dk$，那么：
$$\begin{gathered} \sum _{k=1}^{n}k\sum _{d=1}^a\mu (d)d^2\ast \frac{x(x+1)}{2}\ast \frac{y(y+1)}{x} \\ =\sum _{t=1}^n\sum _{d|t}\mu (d)d^2\ast \frac{t}{d}\ast \frac{\lfloor \frac{n}{t}\rfloor (\lfloor \frac{n}{t}\rfloor +1)}{2}\ast \frac{\lfloor \frac{m}{t}\rfloor (\lfloor \frac{m}{t}\rfloor +1)}{2} \\ =\sum _{t=1}^{n}f(\lfloor \frac{n}{t}\rfloor ,\lfloor \frac{m}{t}\rfloor )\ast t\sum _{d|t}\mu (d)\ast d \end{gathered}$$

我们设一个函数$G(t)=t{\sum }_{d|t}\mu (d)d$，可以证明这个函数是积性的：
当$gcd(a,b)=1$时：
$$\begin{gathered} G(a)G(b) \\ =a\sum _{d_1|a}\mu (d_1)d_1\ast b\sum _{d_2|b}\mu (d_2)d_2 \\ =ab\sum _{d_1|a}\sum _{d_2|b}\mu (d_1)\mu (d_2)d_1{d}_2 \\ =ab\sum _{d_1|a}\sum _{d_2|b}\mu (d_1{d}_2)d_1{d}_2 \\ =ab\sum _{d|ab}\mu (d)d \\ =G(ab) \end{gathered}$$

证毕.

然后我们就可以用线性筛筛这个函数：
1.当$i$为素数时，$G(i)=i(\mu (1)+\mu (i)i)=i+i^2\mu (i)=i-i^2$.
2.当$gcd(i,pr[j])=1$时，根据积性函数的定义，$G(i\ast pr[j])=G(i)G(pr[j])$.
3.当$pr[j]|i$时，可以推导出：
$$G(i\ast pr[j])=i\ast pr[j]\sum _{d|i\ast pr[j]}\mu (d)d$$

考虑将式子拆成两部分，一部分为$d|i$的情况，另一部分为$gcd(d,i)̸=d$的情况，那么容易发现后者的$\mu$值均为$0$，所以：
$$G(i\ast pr[j])=pr[j]\ast i\sum _{d|i}\mu (d)d=pr[j]G(i)$$

那么我们将$G$函数线性筛出来之后，就可以将原始转变为：
$$\sum _{k=1}^{n}k\sum _{d=1}^a\mu (d)d^2\ast \frac{x(x+1)}{2}\ast \frac{y(y+1)}{x}=\sum _{t=1}^{n}f(\lfloor \frac{n}{t}\rfloor ,\lfloor \frac{m}{t}\rfloor )G(t)$$

发现后面的式子在知道$G$函数的取值时是可以用除法分块做到$O(\sqrt{n})$回答一次询问的.

那么我们只需要对$G$函数进行线性筛并前缀和就可以做到$O(n+T\sqrt{n})$解决这个问题了.

代码如下：

#include<bits/stdc++.h>
  using namespace std;
 
#define Abigail inline void
typedef long long LL;
 
const LL mod=100000009;      //把模数看成1e9+9死了一次...
const int N=10000000;
 
LL g[N+9];
int pr[N+9],b[N+9],tp;
 
LL M(LL x){
  x=x%mod+mod;
  return x>=mod?x-mod:x;
}
 
LL f(int x,int y){return M(M(LL(x+1)*x>>1)*M(LL(y+1)*y>>1));}
 
void sieve(int n){
  for (int i=2;i<=n;++i) b[i]=1;
  g[1]=1LL;
  for (int i=2;i<=n;++i){
    if (b[i]) pr[++tp]=i,g[i]=M(LL(1-i)*i);
    for (int j=1;j<=tp&&pr[j]*i<=n;++j){
      b[i*pr[j]]=0;
      if (i%pr[j]) g[i*pr[j]]=M(g[i]*g[pr[j]]);
      else {g[i*pr[j]]=M(g[i]*pr[j]);break;}
    }
  }
  for (int i=1;i<=n;++i)
    g[i]=M(g[i]+g[i-1]);
}
 
LL Query(int n,int m){
  if (n>m) swap(n,m);
  LL ans=0LL;
  for (int l=1,r;l<=n;l=r+1){
    r=min(n/(n/l),m/(m/l));
    ans=M(ans+f(n/l,m/l)*(g[r]-g[l-1]));
  }
  return ans;
} 
 
Abigail work(){
  sieve(N);
}
 
Abigail getans(){
  int T,n,m;
  scanf("%d",&T);
  while (T--){
    scanf("%d%d",&n,&m);
    printf("%lld\n",Query(n,m));
  }
}
 
int main(){
  work();
  getans();
  return 0;
}