BZOJ 2005: [Noi2010]能量采集|容斥原理

做这个题上来就SB..

Ans=∑i=1n∑j=1m2∗gcd(i,j)−1

然后我想成2∗gcd(i,j)的因子个数，各种推式子，各种过不了样例..

gcd(x,y)=i说明，(0,0)到(x,y)能被分成相同的i块，每一块中间都不经过任何整点

龙哥一语道破！！！

求起来类似莫比乌斯反演..很裸的样子

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<queue>
#include<vector>
#include<map>
#include<set>
#include<algorithm>
#include<iostream>
#define N 100010
#define ll long long
using namespace std;
ll ans,n,m,f
;
int main()
{
cin>>n>>m;
if(n>m)swap(n,m);
for(int i=n;i;i--)
{
f[i]=(n/i)*(m/i);
for(int j=i+i;j<=n;j+=i)
f[i]-=f[j];
ans+=f[i]*(2*i-1);
}
cout<<ans;
return 0;
}