Bzoj P4403 序列统计___组合计数+lucas定理

最新推荐文章于 2020-07-16 22:40:13 发布

disPlayLzy_

最新推荐文章于 2020-07-16 22:40:13 发布

阅读量322

点赞数

分类专栏： C++ 排列组合 lucas定理

本文链接：https://blog.csdn.net/Gx_Man_VIP/article/details/99639815

版权

C++ 同时被 3 个专栏收录

368 篇文章 2 订阅

订阅专栏

排列组合

11 篇文章 0 订阅

订阅专栏

lucas定理

2 篇文章 0 订阅

订阅专栏

题目大意：

给定三个正整数N、L和R，统计长度在1到N之间，元素大小都在L到R之间的单调不降序列的数量。输出答案对10^6+3取模的结果。
$1≤N,L,R≤10^9，1≤T≤100，L≤R$

分析：

元素在 $[L, R]$ 选其实相当于在 $[1, R - L + 1]$ 选
问题转换成，
将 $i$ 个( $i \in [1, n]$ )球放到 $R - L + 1$ 个盒子内，盒子可以为空的方案数
$ans=\sum_{i=1}^{n}C_{i+(R-L)}^{R-L}$
也就是 $ans=C_{n+(R-L)+1}^{R-L+1}-1$
求C(n,m).但n,m过大时，可用lucas定理

代码：

#include <iostream>
#include <cstdio>
#include <cmath>
#include <queue>
#include <vector>
#include <cstring>
#include <algorithm>
 
#define rep(i, st, ed) for (int i = st; i <= ed; i++)
#define rwp(i, ed, st) for (int i = ed; i >= st; i--)
#define lson(x) x * 2
#define rson(x) x * 2 + 1
 
#define N 1000005
  
using namespace std;
  
typedef long long ll;
 
const int mo = 1000003;
 
int fac[mo+5], inv[mo+5], ans, T;
 
void Pre_Work() {
    fac[0] = 1; rep(i, 1, mo) fac[i] = (ll)fac[i - 1] * i % mo;
    inv[0] = inv[1] = 1; rep(i, 2, mo) inv[i] = (ll)(mo - mo / i) * inv[mo % i] % mo;
    inv[0] = 1; rep(i, 1, mo) inv[i] = (ll)inv[i - 1] * inv[i] % mo;
}
 
int C(int n, int m) {
    if (n < m) return 0;
    if (n > mo || m > mo) return 1ll * C(n / mo, m / mo) * C(n % mo, m % mo) % mo;
    return 1ll * fac[n] * inv[n - m] % mo * inv[m] % mo;
}
 
int main() {
    Pre_Work();
    scanf("%d", &T);
    int n, L, R;
    while (T--) {
        scanf("%d %d %d", &n, &L, &R);
        R = R - L + 1;
        ans = (C(n + R, R) - 1) % mo;
        ans = (ans % mo + mo) % mo;
        printf("%d\n", ans);
    }
    return 0;
}

HOME Back