bzoj3505（排列组合）

一直对排列组合有点头疼，总是有情况没有考虑到

考虑补集，先把答案赋为c(n*m)(3),然后减去三点共线的情况。

首先三点在一条水平或竖直的直线上非常好处理。直接减去c(n)(3)*m+c(m)(3)*n即可。

然后考虑不斜着的情况。

我们枚举一下边上两个点的横纵坐标之差(i,j)。那么中间的点可选的位置就是gcd(i,j)-1;然后再乘上这种直线的条数即可。注意要先处理组合数，要不然会爆的很惨。

#include<cstdio>
#include<algorithm>
#include<cstring>
#define ll long long 
using namespace std;

ll c[1000009][4],n,m;
ll gcd(ll a,ll b)
{
	return b? gcd(b,a%b):a;
}
int main()
{
	scanf("%lld%lld",&n,&m);n++;m++;
	c[0][0]=1;
	for (int i=1;i<=n*m;i++)
	{
		c[i][0]=1;
		for (int j=1;j<=3;j++)c[i][j]=c[i-1][j-1]+c[i-1][j];
	}
	ll ans=0;
	ans=c[n*m][3];
	ans-=n*c[m][3]+m*c[n][3];
	for (int i=1;i<n;i++)
	for (int j=1;j<m;j++)
	{
		ll t=gcd(i,j);
		if (t>1) ans-=(t-1)*(n-i)*(m-j)*2;//斜着的情况！，先是点的个数和这样线段的条数，和两种斜法
	}
	printf("%lld",ans);
	return 0;
}