P1447 [NOI2010]能量采集（暴力+容斥原理）

题目链接
80分n^2暴力做法

#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio>
using namespace std;
int n,m,ans=0;
int gcd(int a,int b)
{
	return b==0?a:gcd(b,a%b);
}
int main()
{
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++)
	for(int j=1;j<=m;j++)
	{
//		if(gcd(i,j)!=1) ans+=(gcd(i,j)-1)*2+1;
//		if(gcd(i,j)==1) ans++;
		ans+=gcd(i,j)*2-1;
	}
	printf("%d\n",ans);
}

AC（容斥原理）
f[i]记录存在点对(x,y)因数为i的个数
显然f[i]=(n/i)*(m/i);
再用容斥减去它的倍数的情况

#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio>
using namespace std;
int n,m;
long long f[100005],ans=0;
int main()
{
	cin>>n>>m;
	int k=min(n,m);
	for(int i=k;i>=1;i--)
	{
		f[i]=(long long)(n/i)*(m/i);
		for(int j=2;j*i<=k;j++)
		{
			f[i]-=f[i*j]; 
		}
		ans+=f[i]*(2*i-1);
	}
	cout<<ans<<endl;
}