[Lucas定理][组合] Bzoj P4403 序列统计

Description

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

题解

• 我们要求长度为n的元素在[1..r-l+1]中的单调不下降序列的数量
• 考虑将第i个元素加上i-1，那么每个合法序列就会对应一个长度为n，元素为[1..n+m-1]中的单调递增的序列数量
• 显然两者是一一对应的，那么就只用统计后面的贡献
• 那么就是在n+m-1个元素中选n个，就是C(n+m-1,n)
• 所以答案就是C(m,n+m)-1，然后用lucas定理算这个东东就好了

代码

#include <cstdio>
#include <iostream>
#define ll long long
using namespace std;
const ll mo=1e6+3;
ll n,l,r,T,ans,mul[mo],inv[mo];
ll C(ll n,ll m) { return m>n?0:mul[n]*inv[mul[m]]%mo*inv[mul[n-m]]%mo; }
ll lucas(ll n,ll m) { return !m?1:(C(n%mo,m%mo)*lucas(n/mo,m/mo)%mo)%mo; }
int main()
{
    mul[1]=mul[0]=1; for (int i=2;i<=mo;i++) mul[i]=i*mul[i-1]%mo;
    inv[1]=inv[0]=1; for (int i=2;i<=mo;i++) inv[i]=(mo-mo/i)*inv[mo%i]%mo;
    for (scanf("%lld",&T);T;T--) scanf("%lld%lld%lld",&n,&l,&r),ans=((lucas(r-l+1+n,r-l+1)-1)%mo+mo)%mo,printf("%lld
",ans);
}