hiho1804 - 整数分解、组合数、乘法逆元

题目链接

题目叙述很啰嗦，可以简化为：n个球[1-1e5]，放到m个不同的桶里，一共多少种不同的放法。【桶里可以不放】

------------------------------------------------------------------------------------------------

解C（n+m-1, m-1）

由于m,n可能很大，所以需要用逆元。扩展欧几里得。

#include <set>
#include <map>
#include <stack>
#include <queue>
#include <cmath>
#include <vector>
#include <string>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <algorithm>
typedef long long LL;
using namespace std;
const int N = ;
LL MOD = 1e9 + ;
LL F[N], Inv[N];

LL exgcd(LL a, LL b, LL& x, LL& y){
    if (b == ){
        x = , y = ; return a;
    }
    else{
        auto gcd = exgcd(b, a%b, y, x);
        y -= (a / b) * x;
        return gcd;
    }
}

int main(){
    LL n, k; cin >> n >> k;
    int left = n - k, last = MOD, seg = ;
    n += MOD;
    for (int i = ; i < left; ++i){
        int val; scanf("%d", &val);
        if (val>last){
            last = val;
        }
        else{
            ++seg;
            n -= last;
            last = val;
        }
    } n -= last;

    if (n < ){
        puts("No");
    }
    else{
        F[] = ; for (int i = ; i < N; ++i){
            F[i] = (i * F[i - ]) % MOD;
        }
        LL nouse; exgcd(F[N - ], MOD, Inv[N-], nouse);
        while (Inv[N - ] < ) Inv[N - ] += MOD;

        for (int i = N - ; i; --i){
            Inv[i] = (Inv[i + ] * (i + )) % MOD;
        }

        std::cout << (F[n + seg - ] * Inv[n]) % MOD*Inv[seg - ] % MOD << std::endl;
        //C(n+seg-1,seg-1);
    }
}