啥
将相邻点的距离表示为d1,d2…dm,并且和为n
显然当
d
i
+
d
(
i
+
1
)
>
=
(
n
+
1
)
/
2
di+d(i+1)>=(n+1)/2
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);
}
}