【noi】【2010】【能量收集】【莫比乌斯函数】

---

最近终于是把莫比乌斯这个一直搁置着的知识点重拾起来了 ， 而关于由莫比于斯函数μ的特性而产生的二级结论也是略知一二了 ， 那么由此题入手便是再适合不过的了。。。

---

能量收集

栋栋有一块长方形的地，他在地上种了一种能量植物，这种植物可以采集太阳光的能量。在这些植物采集能量后，
栋栋再使用一个能量汇集机器把这些植物采集到的能量汇集到一起。 栋栋的植物种得非常整齐，一共有n列，每列
有m棵，植物的横竖间距都一样，因此对于每一棵植物，栋栋可以用一个坐标(x, y)来表示，其中x的范围是1至n，
表示是在第x列，y的范围是1至m，表示是在第x列的第y棵。 由于能量汇集机器较大，不便移动，栋栋将它放在了
一个角上，坐标正好是(0, 0)。 能量汇集机器在汇集的过程中有一定的能量损失。如果一棵植物与能量汇集机器
连接而成的线段上有k棵植物，则能量的损失为2k + 1。例如，当能量汇集机器收集坐标为(2, 4)的植物时，由于
连接线段上存在一棵植物(1, 2)，会产生3的能量损失。注意，如果一棵植物与能量汇集机器连接的线段上没有植
物，则能量损失为1。现在要计算总的能量损失。 下面给出了一个能量采集的例子，其中n = 5，m = 4，一共有20
棵植物，在每棵植物上标明了能量汇集机器收集它的能量时产生的能量损失。 在这个例子中，总共产生了36的能
量损失。

【Input】

仅包含一行，为两个整数n和m。

【Output】

仅包含一个整数，表示总共产生的能量损失。

【Sample Input】

【样例输入1】
5 4
【样例输入2】
3 4

【Sample Output】

【样例输出1】
36
【样例输出2】
20
对于100%的数据：1 ≤ n, m ≤ 100,000。

---

这道题如果用莫比乌斯函数来做的话似乎不是最快的（233）但是可以进一步熟练莫比乌斯函数的运用和一些二级结论，此处用到的第一个耳机结论便是求：
当 a∈（s_a ， e_a）, b∈（s_b , e_b）时，满足 gcd（a ， b） = 1 的对数个数；
即 ：

∑[i]（1 -> e_a）∑[j]（1 -> e_b）[gcd(i , j) == 1]

由于μ有个性质为：

那么可以把gcd（i ，j） == 1 变为：

Σ[d](d | i && d | j) μ(d)

上面的式子的意思为：
只有当i， j的公因数为一的情况下才对答案贡献 ， 其余的都不会时答案增加：
那么原始便可以变为：

∑[i]（1 -> e_a）∑[j]（1 -> e_b） Σ[d](d | i && d | j) μ(d)

那么就可以化简为：

Σ[d] μ(d) * (e_a / d) * (e_b / d)

,那么就有了最初始的代码：
而本题就是枚举gcd的值，再转化为上面的情况：（代码如下）

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define M 300001
using namespace std;
int n , m;
long long ans = 0;
int miu[M] , ch[M] , qu[M] , tot;
int main(){
    scanf("%d%d",&n,&m);
    if(n > m)swap(n , m);
    miu[1] = 1;
    for(int i = 2 ; i <= n ; i++){
        if(ch[i] == 0){
            miu[i] = -1;
            qu[++tot] = i;
        }
        for(int j = 1 ; j <= tot && qu[j] * i <= n; j++){
            ch[qu[j] * i] = 1;
            if(i % qu[j] == 0)  break;
            else    miu[i * qu[j]] = miu[i] * -1;
        }
    }
    for(int p = 1 ; p <= n ; p++)
        for(int d = 1 ; d <= n / p ; d++)
            ans += 1LL * p * 2 * miu[d] * (n / (p * d)) * (m / (p * d));
    ans -= 1LL * n * m;
    printf("%lld",ans);
}

但是明显在这样做。。。数据不水的话就可能超时 ， 那么就要用到根号优化：
由于每次都要乘一个向下取整的数 ， 而有一些数向下取整的值是一样的 ， 所以可以将这些数一起运算
代码如下：

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define M 300001
using namespace std;
int n , m;
long long ans = 0;
int miu[M] , ch[M] , qu[M] , tot;
int main(){
    scanf("%d%d",&n,&m);
    if(n > m)swap(n , m);
    miu[1] = 1;
    for(int i = 2 ; i <= n ; i++){
        if(ch[i] == 0){
            miu[i] = -1;
            qu[++tot] = i;
        }
        for(int j = 1 ; j <= tot && qu[j] * i <= n; j++){
            ch[qu[j] * i] = 1;
            if(i % qu[j] == 0)  break;
            else    miu[i * qu[j]] = miu[i] * -1;
        }
    }
    for(int i = 1;  i <= n ; i++)miu[i] += miu[i - 1];
    for(int p = 1 ; p <= n ; p++)
        for(int d = 1 ; d <= n / p ;){
            int nex = min((n / p) / ((n / p) / d) , (m / p) / ((m / p) / d));
            ans += 1LL * p * 2 * (miu[nex] - miu[d - 1]) * (n / (p * d)) * (m / (p * d));
            d = nex + 1;
        } 
    ans -= 1LL * n * m;
    printf("%lld",ans);
}