首先，给一个单调不降序列的第i位+i，这样就变成了单调上升序列，设原来数据范围是(l,r)，改过之后变成了(l+1,r+n)
在m个数里选长为n的一个单调上升序列的方案数为( C_m^n )，也就是随便选n个数只能组成惟一的单调上升序列，所以要求的式子就变成了
[sum_{i=1}^{n}C_{r-l+i}^{i}
]
这样看着比较难受，我们把它改成
[sum_{i=1}^{n}C_{r-l+i}^{r-l}
]
我们在开头加一个( C_{r-l+1}^{r-l+1} )，这样根据( C_nm=C_{n-1}m+C_{n-1}^{m-1} )，他就可以和式子的第一项( C_{r-l+1}^{r-l} )合并为( C_{r-l+2}^{r-l+1} )，然后这个又可以可第二个合并，以此类推，最后这个式子就会合并为( C_{r-l+n+1}^{r-l+1} )，然后再减掉( C_{r-l+1}^{r-l+1}=1 )即可
然后用lucas求这个组合数即可
#include<iostream>
#include<cstdio>
using namespace std;
const long long N=1000010,mod=1e6+3;
long long T,n,l,r,fac[N],inv[N];
long long ksm(long long a,long long b)
{
long long r=1;
while(b)
{
if(b&1)
r=r*a%mod;
a=a*a%mod;
b>>=1;
}
return r;
}
long long C(long long n,long long m)
{
if(m>n)
return 0;
return fac[n]*inv[m]%mod*inv[n-m]%mod;
}
long long lucas(long long n,long long m)
{
if(m>n)
return 0;
return n<mod?C(n,m):lucas(n/mod,m/mod)*C(n%mod,m%mod)%mod;
}
int main()
{
scanf("%lld",&T);
fac[0]=1,inv[0]=1;
for(int i=1;i<mod;i++)
fac[i]=fac[i-1]*i%mod;
inv[mod-1]=ksm(fac[mod-1],mod-2);
for(int i=mod-2;i>=1;i--)
inv[i]=inv[i+1]*(i+1)%mod;
while(T--)
{
scanf("%lld%lld%lld",&n,&l,&r);
printf("%lld
",(lucas(r-l+n+1,r-l+1)-1+mod)%mod);
}
return 0;
}