组合数学&&乘法逆元

排列数
A n m = n ! ( n − m ) ! A_n^m = \frac{n!}{(n-m)!} Anm=(nm)!n!
组合数
C n m = n ! m ! ( n − m ) ! C_n^m = \frac{n!}{m!(n-m)!} Cnm=m!(nm)!n!
性质
C n m = C n n − m C_n^m = C_n^{n-m} Cnm=Cnnm
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=Cn1m+Cn1m1
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 mk/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} Cn1(k+1)/21
同理容易得到b的方案数 C m − 1 k / 2 − 1 C_{m-1}^{k/2-1} Cm1k/21
注意下、前面规定的放置顺序是 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;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值