jzoj5983 多边形（组合数学）

啥

将相邻点的距离表示为d1,d2…dm，并且和为n
显然当$di+d(i+1)>=(n+1)/2$时出现一个锐角。
有两结论：k>3时无解，k=3时只有m=3有答案。

考虑求出这样的一组d，那么他对应着n个确定的多边形。 直接求不好求，不如容斥
当k=3时，很容易算出解数就是m==3 ? C(H,2) : 0，可以枚举中间那个。
当k=1时，算至少为两个锐角的有哪些。枚举一个锐角，比如说是di和di+1这个地方。
然后分类讨论，若di>=H，则相当于把n-(H-1)分成m个正整数的方案；
若di=x<=H，则di+1>H-x-1，即相当于把n-x-(H-x-1)分为m-1个正整数的方案。
当k=2时，同理就可以了。要用到组合数的某些求和技巧。
容斥的时候要注意，3个锐角对2个锐角的系数是3,2个锐角对1个锐角的系数是2，三个锐角对一个锐角的系数是3.

#include <cstdio>
#include <iostream>
#include <cstring>
#include <algorithm>
using namespace std;
typedef long long ll;
const ll N = 1e6 + 10, mo = 1000109107;
ll n,m,k,H,ans,T;
ll ksm(ll x,ll y) {
	ll ret = 1; for (; y; y>>=1) {
		if (y&1) ret=ret*x%mo;
		x=x*x%mo;
	}
	return ret;
}

ll S(ll l,ll r) {
	l%=mo,r%=mo; return (l+r)*(r-l+1)%mo*(mo+1>>1)%mo;
}

ll jc[N],njc[N],inv[N];
ll C(ll n,ll m) {
	if (n < m) return 0;
	return jc[n] * njc[m] % mo * njc[n - m] % mo;
}

ll calc1() {
	return (C(n-H,m-2)*(H-1)%mo+C(n-H,m-1))*m%mo;
}

ll calc2() {
	return (C(n-H,m-1)+C(H,m-1))*m%mo;
}

ll calc3() {
	return m==3 ? C(H,2)%mo : 0;
}

int main() {
	freopen("polygon.in","r",stdin);
	// freopen("polygon.out","w",stdout);
	jc[0]=njc[0]=1; for (int i = 1; i < N; i++) jc[i] = jc[i-1] * i % mo;
	njc[N-1]=ksm(jc[N-1],mo-2);
	for (int i = N - 2; ~i; i--) njc[i] = njc[i+1] * (i+1) % mo;
	for (int i = 1; i < N; i++) inv[i] = njc[i] * jc[i-1] % mo;
	for (cin>>T;T;T--) {
		ans=0;
		scanf("%lld %lld %lld",&n,&m,&k), H = (n + 1) / 2;
		// printf("%lld %lld %lld\n",calc1(),calc2(),calc3());
		if (k==3) ans=calc3();
		else if (k==2) ans=calc2()-3*calc3();
		else if (k==1) ans=calc1()-2*(calc2()-3*calc3())-3*calc3();
		if (k==0) ans=C(n-1,m-1)-calc3()-(calc2()-3*calc3())
								-(calc1()-2*(calc2()-3*calc3())-3*calc3());
		ans=ans*inv[m]%mo*n%mo;
		printf("%lld\n",(ans+mo)%mo);
	}
}