【Lucas定理】BZOJ4403[序列统计]题解

最新推荐文章于 2020-02-14 17:23:10 发布

ZigZagK

最新推荐文章于 2020-02-14 17:23:10 发布

阅读量370

点赞数

分类专栏： BZOJ题解 Lucas定理

本文链接：https://blog.csdn.net/zzkksunboy/article/details/78331165

版权

BZOJ题解同时被 2 个专栏收录

117 篇文章 0 订阅

订阅专栏

Lucas定理

3 篇文章 0 订阅

订阅专栏

题目概述

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

解题报告

先考虑长度为 $n$ 方案数。假设现在选好了 $n$ 个数，那么序列是确定的（相同不算），所以直接插板法就可以了，令 $m=R-L+1$ ，则方案数就是 $m+n-1\choose m-1$ 。

那么答案就是 $\sum_{i=1}^{n}{m+i-1\choose m-1}$ 。

添上一个 $m\choose m$ ，由于 ${x\choose y}={x-1\choose y-1}+{x-1\choose y}$ ，所以可以得到：

(m + n - 1 m - 1) + (m + n - 2 m - 1) + \dots + (m m - 1) + (m m) = (m + n - 1 m - 1) + (m + n - 2 m - 1) + \dots + (m + 1 m - 1) + (m + 1 m) = \dots = (m + n m)

${m+n-1\choose m-1}+{m+n-2\choose m-1}+\cdots+{m\choose m-1}+{m\choose m}\\ ={m+n-1\choose m-1}+{m+n-2\choose m-1}+\cdots+{m+1\choose m-1}+{m+1\choose m}\\ =\cdots\\ ={m+n\choose m}$
即答案为

(m+nm)−1 ${m+n\choose m}-1$ 。然后看到了可怕的

109 $10^9$ ……而

106+3 $10^6+3$ 是个素数，直接Lucas定理+乘法逆元。

示例程序

#include<cstdio>
using namespace std;
typedef long long LL;
const int MOD=1000003;

int te,n,L,R,fac[MOD+5],INV[MOD+5];

void Make()
{
    fac[0]=INV[0]=INV[1]=1;for (int i=2;i<MOD;i++) INV[i]=MOD-(LL)(MOD/i)*INV[MOD%i]%MOD;
    for (int i=1;i<MOD;i++) fac[i]=(LL)fac[i-1]*i%MOD,INV[i]=(LL)INV[i-1]*INV[i]%MOD;
}
#define C(x,y) ((x)<(y)?0:((LL)fac[(x)]*INV[(y)]%MOD*INV[(x)-(y)]%MOD))
inline int Locas(int x,int y) {int ans=1;while (y) ans=(LL)ans*C(x%MOD,y%MOD)%MOD,x/=MOD,y/=MOD;return ans;}
int main()
{
    freopen("program.in","r",stdin);
    freopen("program.out","w",stdout);
    for (Make(),scanf("%d",&te);te;te--)
    {
        scanf("%d%d%d",&n,&L,&R);
        printf("%d\n",(Locas(R-L+1+n,n)+MOD-1)%MOD);
    }
    return 0;
}