能量采集energy {数学}

【题目描述】
在一个宽阔的平原上，Samha建立了一个巨大的能量收集阵。在平原上建立一个坐标系，收集点分布在n行m列，坐标从(1,1)到(n,m)，一共n * m个。在(0,0)上有整个收集阵的能量核心，每一个收集点都能为能量核心提供2k+1单位的能量，k是收集点和能量核心连线上收集点的个数（不包含自身）。例如点(2,4)和(0,0)的连线上有一个点(1,2)，那么它提供的能量是3。点(3,5)提供的能量是1。 现在你要求出，所有n*m个收集点一共能提供多少单位的能量。

【Sample Input】（一行，两个数n,m）
5 4
{样例解释}
用*表示核心的位置，各个位置上的数字为该点提供的能量。左上为(0,0)。
*
1 1 1 1 1
1 3 1 3 1
1 1 5 1 1
1 3 1 7 1

【Sample Output】（一行，能量总数）
36

【数据范围】
对于10%的数据：1 ≤ n, m ≤ 10；
对于50%的数据：1 ≤ n, m ≤ 100；
对于80%的数据：1 ≤ n, m ≤ 1000
对于90%的数据：1 ≤ n, m ≤ 10,000；
对于100%的数据：1 ≤ n, m ≤ 100,000。

---

【题解】
本题首先要看某一个点和 (0,0) 连线上有多少个整点。
若有点 (x,y)，且 gcd(x,y)=1，即 x 与 y 互质，则这个点不会被挡住，而其倍数点都会被这个点挡住，这些点就会在一条线上。
所以我们用f[i]表示和 (0,0) 连线上有i个点的点数。满足 i|x && i|y的点对数是 [x/i]*[y/i]，但注意到其中会有重复统计到 f[2 *i], f[3 *i]… 这一类的点。所以从大到小统计，再将重复统计的点扣去即可。
//记得开long long

---

#include <cstdio>
#include <iostream>
#define LL long long
LL n,m,ans,f[1000005];
int main()
{
    scanf("%lld%lld\n",&n,&m);
    if (n<m) std::swap(n,m);
    for (int i=m;i;--i)
    {
        f[i]=(n/i)*(m/i);
        for (int j=2;i*j<=m;++j) f[i]-=f[i*j];
        ans+=f[i]*i;
    }
    printf("%lld\n",ans*2-n*m);
    return 0;
}