【BZOJ2005】【codevs1937】能量采集，数论练习之二维公约数求和

传送门1
传送门2

其实phi可以前缀和处理，(n/d)*(m/d)只有(n^0.5+m^0.5)个取值，然而我这里是O(n)的偏暴力流……
代码：

#include<cstdio>
#include<cstring>
#include<iostream>
#include<cmath>
using namespace std;
int phi[100010],prime[100010];
bool pd[100010];
long long n,m;
long long ans;
main()
{
    scanf("%d%d",&n,&m);
    if (n>m) swap(n,m);
    for (int i=1;i<=n;i++)
    {
        if (!pd[i]) phi[i]=i-1,prime[++prime[0]]=i;
        for (int j=1;j<=prime[0];j++)
        {
            if (prime[j]*i>n) break;
            pd[i*prime[j]]=1;
            if (i%prime[j]==0)
            phi[i*prime[j]]=phi[i]*prime[j];
            else phi[i*prime[j]]=phi[i]*phi[prime[j]];
        }
    }//欧拉筛求φ
    ans=n*m;
    for (int i=1;i<=n;i++)
    ans+=2*(n/i)*(m/i)*phi[i];
    printf("%lld",ans);
}