[Lucas 原理+逆元]BZOJ 4403——序列统计

题目描述

给定三个正整数N、L和R，统计长度在1到N之间，元素大小都在L到R之间的单调不降序列的数量。

输出答案对10^6+3取模的结果。

解题思路

显然L,R的大小对答案没有影响，但是相对大小有影响，设$m=R-L+1$

用隔板法可推出长度为n的方案数$(^{n+m-1}_{m-1})$

于是总方案就是$\sum_{i=1}^{n}(^{i+m-1}_{m-1})$

通过$(^{x}_{y})=(^{x-1}_{y-1})+(^{x-1}_{y})$，我们在总方案后面加上$(^{m}_{m})$

于是就可以得到总方案为$(^{n+m}_{m})-1$

然后直接套用Lucas定理就可以了。

#include<cstdio>
#define LL long long
using namespace std;
const int maxn=1000010,tt=1000003;
LL ji[maxn];
int T,L,R,n;
void work(){ji[0]=1;for (int i=1;i<maxn;i++) ji[i]=(ji[i-1]*i)%tt;}
LL qsm(LL w,int b){
    LL num=1;
    while(b>0){
        if (b%2==1) num=num*w%tt;
        w=w*w%tt;
        b>>=1;
    }
    return num;
}
LL c(int x,int y){if (x>y) return 0;return ji[y]*qsm(ji[x]*ji[y-x]%tt,tt-2)%tt;}
LL Q(int x,int y){LL num=1;while(y) num=num*c(x%tt,y%tt)%tt,x/=tt,y/=tt;return num;}
int main(){
    freopen("exam.in","r",stdin);
    freopen("exam.out","w",stdout);
    scanf("%d",&T);
    work();
    while(T--){
        scanf("%d%d%d",&n,&L,&R);
        printf("%lld\n",(Q(R-L+1,R-L+1+n)-1+tt)%tt);
    }
    return 0;
}