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

题目描述

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

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

解题思路

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

用隔板法可推出长度为n的方案数 (n+m1m1)

于是总方案就是 ni=1(i+m1m1)

通过 (xy)=(x1y1)+(x1y) ,我们在总方案后面加上 (mm)

于是就可以得到总方案为 (n+mm)1

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

#include
#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;i1]*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;
}

你可能感兴趣的:(BZOJ,Lucas定理,CHNJZ的OI学习总结)