BZOJ2154,BZOJ2693 (Mobius反演)

BZOJ2154:Crash的数字表格
BZOJ2693:jzptab

• BZOJ 2154
  给$n,m$,求$\sum\limits_{i=1}^n\sum\limits_{j=1}^n[lcm(i,j)]$

• 将原式适当变形得：

  $$\begin{matrix}\sum_{i=1}^n\sum_{j=1}^n[lcm(i,j)]&=&\sum_{i=1}^n\sum_{j=1}^n\frac{i·j}{gcd(i,j)}\\&=&\sum_d\sum_{i=1}^{\lfloor\frac{n}{d}\rfloor}\sum_{i=1}^{\lfloor\frac{m}{d}\rfloor}\frac{d^2·i·j}{d}[gcd(i,j)==1]\\&=&\sum_d d\sum_{i=1}^{\lfloor\frac{n}{d}\rfloor}\sum_{i=1}^{\lfloor\frac{m}{d}\rfloor}i·j·\sum_{x|i\land x|j}\mu(x)\\&=&\sum_d d\sum_x\mu(x)·x^2·\frac{\lfloor\frac{n}{dx}\rfloor·(\lfloor\frac{n}{dx}\rfloor+1)}{2}·\frac{\lfloor\frac{m}{dx}\rfloor·(\lfloor\frac{m}{dx}\rfloor+1)}{2}\end{matrix}$$

有了上述公式，通过预处理$\mu(x)·x^2$前缀和便可做到$O(\sqrt{n}·\sqrt{n})=O(n)$的单次询问优秀的复杂度。

• BZOJ2693
  给多组n,m,求$\sum\limits_{i=1}^n\sum\limits_{j=1}^n[lcm(i,j)]$

• 明显，再使用刚才的算法会达到$O(n^2)$的复杂度,无法接受。

不妨设$D=dx$.
再次变形得

$$\begin{align}&\sum_d d\sum_x\mu(x)·x^2·\frac{\lfloor\frac{n}{dx}\rfloor·(\lfloor\frac{n}{dx}\rfloor+1)}{2}·\frac{\lfloor\frac{m}{dx}\rfloor·(\lfloor\frac{m}{dx}\rfloor+1)}{2}\\=&\sum_D\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}\sum_{x|D}\frac{D}{i}·\mu(x)·x^2\end{align}$$

不难发现，如果求出$F(D)=\sum_{x|D}\frac{D}{i}·\mu(x)·x^2$的前缀和，单次询问将会在$O(\sqrt{n})$时间内完成。

对于$F(x)$函数，可以通过两种方式实现。
1.枚举法。对于每一个约数枚举其倍数。根据调和数列的性质可以知道该方法预处理时间复杂度为$O(nlog\,n)$。

2.线筛法。
首先有一个结论。
积性函数的约数和为积性函数

因为积性函数的卷积也是积性函数。

此外，当$D=d_1*q,q为d_1约数$时
枚举$x$,当$x$包含$q^2时\mu(x)=0$此时函数无意义。
对于所有有意义的数，其$\frac{D}{x}$项必然扩大$q$倍,$F(D)=F(d_1)*q$

其余的情况都可以转化而来，便得到了线性筛法。

• Code 2693

#include<bits/stdc++.h>
using namespace std;
const int Maxn=1e7;
const int Mod=1e8+9;
typedef long long ll;

inline int read()
{
    char ch=getchar();int i=0,f=1;
    while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
    while(isdigit(ch)){i=(i<<1)+(i<<3)+ch-'0';ch=getchar();}
    return i*f;
}

int Prime[Maxn+50],pr[Maxn+50],tot;
ll F[Maxn+50];

inline void sieve()
{
    F[1]=1;
    for(ll i=2;i<=Maxn;i++)
    {
        if(!pr[i])
        {
            pr[i]=i;
            F[i]=(i-(i*i%Mod))%Mod;
            Prime[++tot]=i;
        }
        for(int j=1;j<=tot;j++)
        {
            int k=i*Prime[j];
            if(k>Maxn)break;
            pr[k]=Prime[j];
            if(i%Prime[j]==0)
            {
                F[k]=F[i]*Prime[j]%Mod;
                break;
            }
            F[k]=F[i]*F[Prime[j]]%Mod;
        }
    }
    for(int i=1;i<=Maxn;i++)(F[i]+=F[i-1])%=Mod;
}

inline int sum(int x,int y)
{
    return (1ll*x*(x+1)/2)%Mod*((1ll*y*(y+1)/2)%Mod)%Mod;
}

int main()
{
    sieve();
    int T=read();
    while(T--)
    {
        int n=read(),m=read();
        if(n>m)swap(n,m);
        int pos,ans=0;
        for(int bg=1;bg<=n;bg=pos+1)
        {
            pos=min(n/(n/bg),m/(m/bg));
            ans=(1ll*ans+1ll*sum(n/bg,m/bg)*(F[pos]-F[bg-1])%Mod)%Mod;
        }
        printf("%d\n",(ans+Mod)%Mod);
    }
}