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

来源:互联网 发布:广联达计价软件多少钱 编辑:程序博客网 时间:2024/05/16 06:25

题目描述

给定三个正整数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<cstdio>#define LL long longusing 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;}
原创粉丝点击