2005: [Noi2010]能量采集

#include <iostream>
#include <algorithm>
#include <cstring>
#include <cstdio>
#include <cmath>
using namespace std;
long long f[110000];
int n,m;
long long ans;
int main()
{
scanf("%d%d",&n,&m);
if(n>m) swap(n,m);
for(int i=n;i;i--)
{
f[i]=(long long)(n/i)*(m/i);
for(int j=2;j*i<=n;j++)
{
f[i]-=f[j*i];
}
ans+=f[i]*(2*i-1);
}
printf("%lld\n",ans);
return 0;
}