BZOJ3505: [Cqoi2014]數三角形(洛谷P3166)
阿新 • • 發佈:2018-12-14
組合
直接計算比較煩,考慮全部方案減不合法的方案。不合法的方案就是三個點都在一條直線上。除去橫的和豎的線,其它的枚一枚求一求gcd乘一乘就好了。
程式碼:
#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 1005
using namespace std;
typedef long long LL;
LL n,m,ans;
#define C(x) (x*(x-1)*(x-2)/6)
LL gcd(LL a,LL b){ return !b?a:gcd(b,a%b); }
int main (){
scanf("%lld%lld",&n,&m),n++,m++;
ans=C(n*m)-m*C(n)-n*C(m);
for (LL i=1;i<n;i++)
for (LL j=1,r;j<m;j++)
if ((r=gcd(i,j)+1)>2) ans-=(r-2)*(n-i)*(m-j)*2;
return printf("%lld\n",ans),0;
}