1. 程式人生 > >BZOJ3505: [Cqoi2014]數三角形(洛谷P3166)

BZOJ3505: [Cqoi2014]數三角形(洛谷P3166)

組合

直接計算比較煩,考慮全部方案減不合法的方案。不合法的方案就是三個點都在一條直線上。除去橫的和豎的線,其它的枚一枚求一求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; }