bzoj 4403: 序列统计 lucas定理+组合数学

题意

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

分析

一开始没想到把不降转为递增来求。
对于一个长度为n的不降序列，若我们把第i个数加上i，那么它就是一个递增序列。也就是说，长度为i且元素在[1,m]之间的不降序列数量就相当于元素在[1+1,m+n]之间的递增序列数量，也就等于$C_{n+m-1}^{m-1}$。
那么我们要求的就是
$\sum\limits_{i=1}^nC_{i+m-1}^{m-1}$
$=(\sum\limits_{i=1}^nC_{i+m-1}^{m-1})+C_m^m-1$
$=(\sum\limits_{i=2}^nC_{i+m-1}^{m-1})+C_{m+1}^m-1$
$=(\sum\limits_{i=3}^nC_{i+m-1}^{m-1})+C_{m+2}^m-1$
…
$=C_{n+m}^m-1$
然后就可以直接套lucas定理了。

代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;

typedef long long LL;

const int MOD=1000003;

int ny[MOD],jc[MOD];

int C(int n,int m)
{
    if (n<m) return 0;
    if (n<MOD&&m<MOD) return (LL)jc[n]*ny[m]%MOD*ny[n-m]%MOD;
    else return (LL)C(n/MOD,m/MOD)*C(n%MOD,m%MOD)%MOD;
}

int main()
{
    jc[0]=ny[0]=jc[1]=ny[1]=1;
    for (int i=2;i<MOD;i++) jc[i]=(LL)jc[i-1]*i%MOD,ny[i]=(LL)(MOD-MOD/i)*ny[MOD%i]%MOD;
    for (int i=2;i<MOD;i++) ny[i]=(LL)ny[i]*ny[i-1]%MOD;
    int T;
    scanf("%d",&T);
    while (T--)
    {
        int n,l,r;
        scanf("%d%d%d",&n,&l,&r);
        int m=r-l+1;
        printf("%d\n",(C(n+m,m)+MOD-1)%MOD);
    }
    return 0;
}