组合数学 - 序列统计 - AcWing 1312

组合数学 - 序列统计 - AcWing 1312

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

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

输入格式

输入第一行包含一个整数 T，表示数据组数。

第二到第 T+1 行每行包含三个整数 N,L,R。

输出格式

输出包含 T 行，每行有一个数字，表示你所求出的答案对 10⁶+3 取模的结果。

数据范围

$$\begin{gathered} 1\le N,L,R\le 10^9,1\le T\le 100， \\ 输入数据保证L\le R。 \end{gathered}$$

输入样例：

2
1 4 5
2 4 5

输出样例：

2
5

样例解释

对于第一组输入，满足条件的两个序列为 {4},{5}。

---

分析：

$本题目标：求满足$

$$L\le a_1\le a_2\le ...\le a_k\le R，k\in [1,N]①$$

$的序列(a_1,a_2,...,a_k)的数量。$

$等价于：$

$$0\le a_1\le a_2\le ...\le a_k\le R-L，k\in [1,N]②$$

转化方法一(隔板法)：

参考：组合数学(隔板法) + 高精度 - 方程的解 - AcWing 1308

$在《方程的解》这道题中，我们利用隔板法解决了：$

$$x_1+x_2+...+x_k=C，C是常数且x_i>0，$$

$形如这样的问题可以用隔板法解决，计算出不同的$正整数解$的方案总数。$

$本题，我们记：x_1=a_1,x_2=a_2-a_1,...,x_k=a_k-a_{k-1}$

$问题转化为：$

$$x_1+x_2+...+x_k\le C=R-L，x_i\ge 0，$$

$求不同的$非负正整数解$的方案总数。$

$我们先映射y_i=x_i+1，将非负正整数解转化为正整数解，$

$不等式转化为：$

$$y_1+y_2+...+y_k\le C=R-L+k，y_i>0，$$

$求不同的正整数解的方案总数。$

$到这里，我们又能够参照《方程的解》，用隔板法解决这个问题。$

需要注意的是：

$因为这里是小于等于号，故我们不一定要把C=R-L+k个小球全部划分完，$

$这里我们用k个隔板(而非k-1个隔板)来解决这个问题。$

$y_i的大小为第i个隔板与第i-1个隔板之间的小球数量。$

$如上图，y_1+y_2+y_3=4\le C=5，为一组可行解。$

$总结：对于每一个序列长度k，方案总数为在R-L+k个空位中选择k个放隔板，即C_{R-L+k}^k$

$最终答案：$

$$\sum _{k=1}^N{C}_{R-L+k}^k$$

---

转化方法二：

$先考虑如下不等式：$

$$L\le b_1<b_2<...<b_k\le R，k\in [1,N]$$

$上述不等式的解的方案总数，相当于在区间[L,R]之间选择k个不同的数即C_{R-L+1}^k。$

$将①式映射:$

$令b_1=a_1,b_2=a_2+1,b_3=a_3+2,...,b_k=a_k+k-1$

$得到新的不等式$

$$L\le b_1<b_2<...<b_k\le R+k-1，k\in [1,N]③$$

$这相当于我们将数列a的相邻两项之间的差增加的1，因为原来数列a相邻两项的差大于等于0。$

$根据方程③，即在区间[L,R+k-1]之间选择k个数，方案总数为：C_{R-L+k}^k，与方法一得到的结论相同。$

---

最后一步转化：

$因为N\le 10^9，若直接枚举k会TLE。$

$故我们对{\sum }_{k=1}^N{C}_{R-L+k}^k进行转化，记M=R-L$

${\sum }_{k=1}^N{C}_{M+k}^k=C_{M+1}^1+C_{M+2}^2+C_{M+3}^3+...+C_{M+N}^N$

$=C_{M+1}^{M+1-1}+C_{M+2}^{M+2-2}+C_{M+3}^{M+3-3}+...+C_{M+N}^{M+N-N}$

$=C_{M+1}^M+C_{M+2}^M+C_{M+3}^M+...+C_{M+N}^M$

$根据C_n^m=C_{n-1}^{m-1}+C_{n-1}^m$

$原式=C_{M+1}^{M+1}+C_{M+1}^M+C_{M+2}^M+C_{M+3}^M+...+C_{M+N}^M-C_{M+1}^{M+1}$

$=C_{M+2}^{M+1}+C_{M+2}^M+C_{M+3}^M+...+C_{M+N}^M-C_{M+1}^{M+1}$

$=....$

$=C_{M+N+1}^{M+1}-C_{M+1}^{M+1}=C_{M+N+1}^{M+1}-1$

故最终答案为：

$$C_{R-L+N+1}^{R-L+1}-1$$

$模数10^6+3是质数，用Lucas定理求组合数。$

---

总结：$线性方程/不定方程求解的方案总数，都可以通过一系列映射，再用隔板法解决。$

---

代码：

#include<iostream>

#define ll long long

using namespace std;

const int mod=1e6+3;

int T,n,l,r;

int quick_pow(int a,int b,int mod)
{
    int res=1;
    while(b)
    {
        if(b&1) res=(ll)res*a%mod;
        a=(ll)a*a%mod;
        b>>=1;
    }
    return res;
}

int C(int a,int b,int mod)
{
    if(a<b) return 0;
    
    int up=1, down=1;
    for(int i=a,j=1;j<=b;i--,j++)
    {
        up=(ll)up*i%mod;
        down=(ll)down*j%mod;
    }
    return (ll)up*quick_pow(down,mod-2,mod)%mod;
}

int Lucas(int a,int b,int p)
{
    if(a<p&&b<p) return C(a,b,p);
    return (ll)C(a%p,b%p,p)*Lucas(a/p,b/p,p)%p;
}

int main()
{
    cin>>T;
    while(T--)
    {
        cin>>n>>l>>r;
        cout<<(Lucas(r-l+n+1,r-l+1,mod)-1+mod)%mod<<endl;
    }
    
    return 0;
}