[bzoj2226] LCMSum 数学+筛法

题目链接：bzoj2226 .

—————————————-

概述

题目大意如下：
给定正整数$n$，求

$$\sum_{i=1}^nlcm(i,n).$$

题目有多组数据。$(T\leq 300000,\ n\leq 1000000)$

$lcm(x,y)$即$x,y$的最小公倍数。

—————————————-

分析

对于$lcm(i,n)$，我们很显然可以把它转化成如下形式：

$$lcm(i,n)=\frac{i\times n}{gcd(i,n)}.$$

对于某个$i$，我们令$d=gcd(i,n)$，$n=d\times x$，$i=d\times y$，那么$x$与$y$是互质的。

所以，我们不妨枚举$d$，再枚举所有与$x$互质的$y$并统计答案，其中$x=\cfrac{n}{d}.$

我们得到下式：

$$Ans\ =\ \sum_{d|n}\sum_{y=1}^{\frac{n}{d}}\frac{n}{d}\times y\times d\times [gcd(\frac{n}{d}, y)=1]\ =\ n\sum_{d|n}\sum_{y=1}^{\frac{n}{d}}y\times [gcd(\frac{n}{d},y)=1].$$

化到这里，问题就变成了：对于$n$所有的约数$d$，计算小于$\cfrac{n}{d}$且与其互质的数的和，将结果$\times n$就是答案了。

那么关键来了，所有与$\cfrac{n}{d}$互质的数的和怎么求呢？

——————–

结论： $$\sum_{i=1}^ni\times [gcd(i,n)=1]\ =\ \cfrac{n\times \varphi(n)}{2}.$$

特别地，$n=1$时值为1.

证明：

引理： 对于一个与$n$互质的数$x$，$n-x$也一定与$n$互质。

引理证明：设$gcd(n-x,n)=d$，那么$\cfrac{n-x}{d}$一定是整数，即$\cfrac{n}{d}-\cfrac{x}{d}$是整数。由于$d$是$n$的约数，$x$与$n$互质，所以$\cfrac{x}{d}$是整数当且仅当$d=1$，所以$gcd(n-x,n)=d=1$。引理证毕。

我们要求的是与$n$互质的所有数的和，即：

$$\sum_{i=1}^{n}i\times [gcd(i,n)=1].$$

根据引理，我们又有：

$$\sum_{i=1}^{n}i\times [gcd(i,n)=1]\ =\ \sum_{i=1}^{n\over 2}(i+(n-i))\times [gcd(i,n)=1]\ =\ n\sum_{i=1}^{n\over 2}[gcd(i,n)=1]\ =\ \cfrac{n\times \varphi(n)}{2}.$$

结论证毕。

——————–

有了这个结论就好办了，我们原来要求的是：

$$Ans=n\sum_{d|n}\sum_{y=1}^{\frac{n}{d}}y\times [gcd(\frac{n}{d},y)=1]$$

现在我们可以将它转化为：

$$Ans\ =\ n(1+\sum_{d|n}\cfrac{\frac{n}{d}\times \varphi(\frac{n}{d})}{2})\ =\ \cfrac{n}{2}(1+\sum_{d|n}d\times \varphi(d)).$$

其中，加$1$是因为当$d=1$时本应该有值为$1$，然而若按上面的式子计算$d=1$时为0，所以要加$1$.

至此，有2种做法，一个是80分，一个是100分。

---

80分：

待求：

$$Ans\ =\ \cfrac{n}{2}(1+\sum_{d|n}d\times \varphi(d)).$$

我们可以线性筛出所有的$\varphi(d)$，对于每个询问$O(\sqrt n)$枚举$n$所有约数$d$并统计答案。总复杂度$O(T\sqrt n).$

代码

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cmath>
#define ll long long
#define For(i,j,k) for(register int i=j; i<=(int)k; ++i)
#define Forr(i,j,k) for(register int i=j; i>=(int)k; --i)
#define INF 0x3f3f3f3f
using namespace std;

const int maxn = 1000000+5;

int T, n, tot;
int pri[maxn/5], phi[maxn];
bool vis[maxn];

inline void read(int &x){
    x = 0;
    char c = getchar();
    while(!isdigit(c))    c = getchar();
    while(isdigit(c))    x = (x<<3)+(x<<1)+(c^48),  c = getchar();
}

inline void init(){
    phi[1] = 1;
    For(i, 2, maxn-1){
        if(!vis[i]){
            pri[++tot] = i;
            phi[i] = i-1;
        }
        for(int j=1,x; j<=tot&&(x=i*pri[j])<maxn; ++j){
            vis[x] = true;
            if(i%pri[j] == 0){
                phi[x] = phi[i]*pri[j];
                break;
            }
            else    phi[x] = phi[i]*(pri[j]-1);
        }
    }
}

inline ll calc(int n){
    ll sum, back;
    sum = 1,  back = 1ll * n;
    for(int i = 1; i*i<=n; ++i)
        if(n%i == 0){
            sum += 1ll * phi[i] * i >> 1;
            if(i*i == n)    break;
            sum += 1ll * phi[n/i] * (n/i) >> 1;
        }
    back *= sum;
    return back;
}

int main(){
    init();
    read(T);
    while(T --){
        read(n);
        printf("%lld\n", calc(n));
    }
    return 0;
}

---

100分

待求：

$$Ans\ =\ \cfrac{n}{2}(1+\sum_{d|n}d\times \varphi(d)).$$

我们令$h(n)=\sum_{d|n}d\times \varphi(d)$，那么答案即是:

$$Ans=\cfrac{n}{2}(1+h(n)).$$

我们观察一下$h(n)=\sum_{d|n}d\times \varphi(d)$，由于$d$与$\varphi(d)$都是积性函数，所以$h(n)$也是积性函数。

既然$h(n)$是积性函数，那么就可以通过线性筛筛出来，我们思考一下怎么筛。

设$p$是一个质数，那么有$h(p)=p\times (p-1)+1.$
那么我们接下来要解决的问题是$h(a\times p)$如何筛得。

总共分2种情况：

1. $gcd(a,p)=1$，即$a,p$互质。这种情况很好处理： $$h(a\times p)=h(a)\times h(p).$$

2. $gcd(a,p)\neq 1$.此时又分两种小情况：

   （1）$a=x\times p^r.$ 那么，$h(a\times p)=h(x)\times h(p^{r+1}).$

   （2）$a=p^r.$ 这种情况下要稍微推一下式子：

   $$已知：h(a)=h(p^r)=1+\sum_{i=1}^rp^i\times \varphi(p^i).$$

   我们对它进行进一步转化：

   $$h(a)=1+\sum_{i=1}^rp^i\times (p^i-p^{i-1})=1+(p-1)\sum_{i=1}^rp^{2i-1}.$$

   那么同理：

   $$h(a\times p)=h(p^{r+1})=1+(p-1)\sum_{i=1}^{r+1}p^{2i-1}.$$

   对于$h(p^r)$，我们可以发现如下的关系：

   $$\qquad \qquad \qquad \qquad \qquad \qquad h(p^r)-1=(p-1)\sum_{i=1}^{r}p^{2i-1}$$

$$\qquad \qquad \qquad \qquad \qquad \quad (h(p^r)-1)\times p^2=(p-1)\sum_{i=2}^{r+1}p^{2i-1}$$

$$\qquad \qquad \quad \ (h(p^r)-1)\times p^2+(p-1)\times p=(p-1)\sum_{i=1}^{r+1}p^{2i-1}$$

$$1+(h(p^r)-1)\times p^2+(p-1)\times p=h(p^{r+1})$$

$$\qquad \qquad \qquad \ \ \ 1+h(p^r)\times p^2-p=h(p^{r+1})$$

至此我们得到了$h(p^r)$与$h(p^{r+1})$的关系，也就是$h(a)$与$h(a\times p)$的关系。

整理一下：

$$h(a\times p)=\begin{cases} h(a)\times h(p), & \text{gcd(a,p)=1} \\ h(x)\times h(p^{r+1}), & \text{gcd(a,p)$\neq$1, a=$x\times p^r$}\\ h(a)\times p^2-p+1, &\text{gcd{a,p}$\neq$1, a=$p^r$} \end{cases}$$

至此，线性筛筛出所有的$h(n)$，对于每一个询问，$Ans=\cfrac{(h(n)+1)\times n}{2}$可以$O(1)$计算。总复杂度$O(n+T).$

代码

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cmath>
#define ll long long
#define For(i,j,k) for(register ll i=j; i<=(ll)k; ++i)
#define Forr(i,j,k) for(register ll i=j; i>=(ll)k; --i)
#define INF 0x3f3f3f3f
using namespace std;

const ll maxn = 1000000+5;

ll T, n, tot;
ll pri[maxn/5], phi[maxn];
bool vis[maxn];

inline void read(ll &x){
    x = 0;
    char c = getchar();
    while(!isdigit(c))    c = getchar();
    while(isdigit(c))    x = (x<<3)+(x<<1)+(c^48),  c = getchar();
}

inline void init(){
    phi[1] = 1;
    For(i, 2, maxn-1){
        if(!vis[i]){
            pri[++tot] = i;
            phi[i] = (i-1) * i + 1;
        }
        for(ll j=1,x; j<=tot&&(x=i*pri[j])<maxn; ++j){
            vis[x] = true;
            if(i%pri[j] == 0){
                ll ii = i,  jj = pri[j];
                while(ii%pri[j] == 0){
                    ii /= pri[j];
                    jj *= pri[j];
                }
                if(ii > 1)    phi[x] = phi[ii] * phi[jj];
                else    phi[x] = phi[i] * pri[j]*pri[j] - pri[j]+1;
                break;
            }
            else    phi[x] = phi[i]*phi[pri[j]];
        }
    }
}

int main(){
    init();
    read(T);
    while(T --){
        read(n);
        printf("%lld\n", 1ll * (phi[n]+1 >> 1) * n);
    }
    return 0;
}

—————————————-

小结

这一题的解决关键在于中间那个欧拉函数的结论，利用那个结论可以把式子化成能够用数论函数解决的形式。之后的重点是推导$h(n)$如何通过线性筛筛出来，情况的考虑并不复杂，但是式子不是很好推。总之，这是一道练习推式子以及筛法的好题。

—————————————-

$\mathscr ——wrote\ by\ miraclejzd.$