题目链接
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;
}