排列数
A
n
m
=
n
!
(
n
−
m
)
!
A_n^m = \frac{n!}{(n-m)!}
Anm=(n−m)!n!
组合数
C
n
m
=
n
!
m
!
(
n
−
m
)
!
C_n^m = \frac{n!}{m!(n-m)!}
Cnm=m!(n−m)!n!
性质
C
n
m
=
C
n
n
−
m
C_n^m = C_n^{n-m}
Cnm=Cnn−m
C
n
m
=
C
n
−
1
m
+
C
n
−
1
m
−
1
C_n^m = C_{n-1}^m + C_{n-1}^{m-1}
Cnm=Cn−1m+Cn−1m−1
C
n
0
+
C
n
1
+
C
n
2
+
.
.
.
.
.
.
C
n
n
=
2
n
C_n^0 + C_n^1 + C_n^2 + ......C_n^n = 2^n
Cn0+Cn1+Cn2+......Cnn=2n
例题
添加链接描述
思路
显然的特判 如果
n
+
m
<
k
n+m<k
n+m<k 或者
k
=
1
k=1
k=1 输出肯定是
0
0
0
否则呢? 先按照
a
b
a
b
a
b
ababab
ababab顺序,填充满
k
k
k个,那么容易得到a和b剩余的个数分别为
n
−
(
k
+
1
)
/
2
n-(k+1)/2
n−(k+1)/2、
m
−
k
/
2
m-k/2
m−k/2。
对于
a
a
a来说,也就是在
(
k
+
1
)
/
2
(k+1)/2
(k+1)/2个已经放好的a内,在放剩下的
n
−
(
k
+
1
)
/
2
n-(k+1)/2
n−(k+1)/2,.
那这不就是隔板法吗,所以方案数就是
C
n
−
1
(
k
+
1
)
/
2
−
1
C_{n-1}^{(k+1)/2 - 1}
Cn−1(k+1)/2−1
同理容易得到b的方案数
C
m
−
1
k
/
2
−
1
C_{m-1}^{k/2-1}
Cm−1k/2−1
注意下、前面规定的放置顺序是
a
b
a
b
a
b
ababab
ababab
那么我们也可以按照babababa顺序放置。。方法和上述一样也是分别进行隔板法然后相乘。
code
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll N = 1e5 + 5;
const ll mod = 1e9+7;
ll n, m, k, f[N], ans;
ll qpow(ll a, ll b)
{
ll res = 1 % mod;
for(; b; b >>= 1)
{
if(b & 1)
res = (ll) res * a % mod;
a = (ll) a * a % mod;
}
return res;
}
int main()
{
cin >> n >> m >> k;
if(n + m < k || k == 1)
cout << 0 << endl;
else
{
f[0]=1;
for(int i = 1; i <= N; i++)
f[i] = f[i-1] * i % mod;
//abababa
if(n >= (k/2 + k%2) && m >= (k/2))
ans = f[n-1] * qpow(f[(k+1)/2 - 1] * f[n - (k+1)/2] % mod, mod-2) % mod * (f[m-1] * qpow(f[k/2 - 1] * f[m - (k/2)] % mod, mod-2) % mod);
//bababab
if(m >= (k/2 + k%2) && n >= (k/2))
ans += f[m-1] * qpow(f[(k+1)/2 - 1] * f[m - (k+1)/2] % mod, mod-2) % mod * (f[n-1] * qpow(f[k/2 - 1] * f[n - (k/2)] % mod, mod-2) % mod);
/*if(n>=k/2+k%2 && m>=k/2) ans=f[n-1]*qpow(f[k/2+k%2-1]*f[n-k/2-k%2]%mod,mod-2)%mod*(f[m-1]*qpow(f[k/2-1]*f[m-k/2]%mod,mod-2)%mod);
//bababab...
if(m>=k/2+k%2 && n>=k/2) ans+=f[m-1]*qpow(f[k/2+k%2-1]*f[m-k/2-k%2]%mod,mod-2)%mod*(f[n-1]*qpow(f[k/2-1]*f[n-k/2]%mod,mod-2)%mod);*/
cout << ans%mod << endl;
}
return 0;
}