bzoj 2005: [Noi2010]能量采集

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