洛谷P1680 奇怪的分组（组合数学 + Lucas定理）

2020.6.15
这道题是随机跳到的，一看思路非常显然，不就是Piot天天讲的“star and bar problem”么，高中的时候老谭其实也讲过，叫隔板法，印象深刻（毕竟那次期末炸了）。分组非空，且对特定组的放置有最低要求，可以转为如下问题，在n个坑里先放上必须要放的数目。总数减去必须要放的就是自由分组的部分，然后用隔板法公式C(n - 1, m - 1)即可求出正确答案。

剩下用然后写组合数取模实现花了半小时，抄模板都没抄对。。。。我太菜+1. 用rep的m范围应该是m-1而不是m，wa了9个点，排完之后就ac了。卢卡斯和普通运算都翘了翘

代码：

#include <bits/stdc++.h>
using namespace std;
#define limit (1000 + 5)//防止溢出
#define INF 0x3f3f3f3f
#define inf 0x3f3f3f3f3f
#define lowbit(i) i&(-i)//一步两步
#define EPS 1e-6
#define FASTIO  ios::sync_with_stdio(false);cin.tie(0);
#define ff(a) printf("%d\n",a );
#define pi(a,b) pair<a,b>
#define rep(i, a, b) for(int i = a; i <= b ; ++i)
#define per(i, a, b) for(int i = b ; i >= a ; --i)
#define mint(a,b,c) min(min(a,b), c)
#define MOD 998244353
#define FOPEN freopen("C:\\Users\\tiany\\CLionProjects\\acm_01\\data.txt", "rt", stdin)
typedef long long ll;
typedef unsigned long long ull;
ll read(){
    ll sign = 1, x = 0;char s = getchar();
    while(s > '9' || s < '0' ){if(s == '-')sign = -1;s = getchar();}
    while(s >= '0' && s <= '9'){x = x * 10 + s - '0';s = getchar();}
    return x * sign;
}//快读
void write(ll x){
    if(x < 0) putchar('-'),x = -x;
    if(x / 10) write(x / 10);
    putchar(x % 10 + '0');
}
const int mod = 1000000007;
ll quickPow(ll base, ll expo){
    ll ans = 1;
    while (expo){
        if(expo & 1)(ans *= base) %= mod;
        expo >>= 1;
        base = base * base;
        base %= mod;
    }
    return ans % mod;
}
ll C(ll n, ll m){
    if(n < m)return 0;
    ll x = 1, y = 1;
    if(m > n - m)m = n - m;
    rep(i ,0 , m - 1){
        x = x * (n - i) % mod;
        y = y *  (i + 1) % mod;
    }
    return x * quickPow(y, mod - 2) % mod;
}
ll lucas(ll n, ll m){
    return !m || n == m ? 1 : C(n % mod, m % mod) * lucas(n / mod, m / mod) % mod;
}
ll n,k;
int main() {
#ifdef LOCAL
    FOPEN;
#endif
    n = read(), k = read();
    rep(i ,1 ,k){
        n -= read();
    }
    write(C(n - 1 , k - 1));
    return 0;
}