【LOJ10235】【BZOJ4403】序列统计

题目链接 

解法:

          首先,我们设f[i][j]表示满足条件的序列中,长度为i,最后一个数字为L+j-1的序列个数。

          那我们就可以得到一下递推式: f[i][j]=\sum_{k=1}^{j} f[i-1][k](i>1,j>1)(其中对于每一个f[1][n]或者f[n][1]其值都为1),

                                              近一步推到 f[i][j]=f[i-1][j]+f[i][j-1].

           且我们所要求的答案是 \sum_{i=1}^{N}\sum_{j=1}^{R-L+1}f[i][j].

           我们考虑化简答案的式子,由递推式得:f[i][R-L+1]=\sum_{j=1}^{R-L+1}f[i-1][j],

                                                      所以原答案式= \sum_{i=2}^{N+1}f[i][R-L+1].

                                                                            =  \sum_{i=2}^{N+1}f[i][R-L+1].+f[1][R-L+2]-1

            又由 f[i][j]=f[i-1][j]+f[i][j-1]得,原式= f[N+1][R-L+2]-1

            接下来引出一个结论:f[i][j]=C_{i+j-2}^{i-1},

            以下是证明: 在一个正方形中,对于每一条从(1,1)走到(i,j)的最短路,(相邻格子移动那种)都可以把第i行的最后一个经过的点的纵坐标,看作是序列的第i位的值,而这样的最短路方案数正是C_{i+j-2}^{i-1},且符合f[i][j]的定义。

            所以综上,答案为C_{N+R-L+1}^{N}.

(哇。。。。我自己都觉得我的解题思路好烦啊。。。。)

 

Code:

#include<bits/stdc++.h>
#define rep(i,j,k) for(int i=j;i<=k;i++)
#define rep2(i,j,k) for(int i=j;i>=k;i--)
#define mo 1000003
using namespace std;
template<typename T> void read(T &num){
	char c=getchar();num=0;T f=1;
	while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
	while(c>='0'&&c<='9'){num=(num<<3)+(num<<1)+(c^48);c=getchar();}
	num*=f;
}
template<typename T> void qwq(T x){
	if(x>9)qwq(x/10);
	putchar(x%10+'0');
}
template<typename T> void write(T x){
	if(x<0){x=-x;putchar('-');}
	qwq(x);putchar('\n');
}
int powe[1000010];int inv[1000010];
inline int poww(int n,int m){
	int re_value=1;int temp=n;int tmp=m;
	while(tmp){
		if(tmp&1)re_value=(1ll*re_value*temp)%mo;
		temp=(1ll*temp*temp)%mo;
		tmp>>=1;
	}
	return re_value;
}
inline void Pretreatment(){
	powe[0]=powe[1]=inv[0]=inv[1]=1;
	rep(i,2,mo-1){powe[i]=(1ll*powe[i-1]*i)%mo;}
	inv[mo-1]=poww(powe[mo-1],mo-2);
	rep2(i,mo-2,2){inv[i]=(1ll*inv[i+1]*(i+1))%mo;}
	return;
}
inline int C(int n,int m){
	if(n<m)return 0;
	return (1ll*((1ll*powe[n]*inv[m])%mo)*inv[n-m])%mo;
}
inline int Lucas(int n,int m){
	if(!m||n==m)return 1;
	return (1ll*C(n%mo,m%mo)*Lucas(n/mo,m/mo))%mo;
}

int main(){
	int T;read(T);
	Pretreatment();
	while(T--){
		int N,L,R;read(N);read(L);read(R);
		int K=R-L+1;
		write((Lucas(N+K,N)-1+mo)%mo);
	}
	return 0;
}

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值