bzoj2005 能量采集

题目链接：http://www.lydsy.com/JudgeOnline/problem.php?id=2005

题目大意：在左边(0,0)处有一个发射器，在横坐标1~x，纵坐标1~y这xy点有接收器，发射器要对每一个点发射信息，发射到(x0,y0)的能量损失为(0,0)(x0,y0)这条线段上的点数(不包括端点)*2+1,问能量损失之和

题目分析：一个点(x0,y0)和点(0,0)中间隔的点的数量为gcd(x0,y0)，之后用容斥原理求解即可

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<iostream>
#include<algorithm>

using namespace std;
typedef long long LL;
LL n,m,ans=0;
LL f[1100000];
int main()
{
	cin>>n>>m;
	LL d=min(n,m);
	for (LL i=d;i>=1;i--)
	{
		f[i]=(m/i)*(n/i);
		LL dd=d/i;
		for (LL j=2;j<=dd;j++)
		    f[i]-=f[i*j];
		ans+=f[i]*(2*i-1);
	}
	cout<<ans<<endl;
	return 0;
}