返回顶部

BZOJ 4403序列统计


假设存在一个满足条件的长度为i的不下降序列(显然是一定存在的)那么只需要从中选出i个数即可

(不必在意选出具体数的大小,可以把满足条件的序列写下来,选几个数感受一下)。

但是(nm)里的 m 的是就是 (rl+1) 吗?

乍一看是这样的,但是这样会出现一个问题,单调不下降子序列中的数可以重复,但如果只是从lr去选数,那将会使结果变小。

所以可以先找出(当然是在脑子里找出)满足条件的长度为i的序列,然后和lr这几个数合起来,从这几个数里面选出i个数。即 (rl+1+irl+1)

然后就从1n枚举长度再求和即可

所以答案为

ni=1( r  l + 1 + i i)

但是此时O(n+)的时间复杂度一定是不可行的,而i的值又一直发生改变,无法化简。

所以我们可以把 (rl+1+ii) 写为(rl+1+irl+1),再设x=rl+1那么原式变为

ni=1(x+i i)

也就是 (x+1x)+(x+2x)+(x+3x)+……+(x+nx)

接下来就是化简……

首先看一下(nm) +(nm+1)的结果是什么

(nm)+(nm+1)

=n!m!(nm)!+n!(m+1)!(nm1)!

=n!(m+1)(m+1)!(nm!)+n!(nm)(m+1)!(nm)!

=(m+1+nm)n!(m+1)!(nm)!

=(n+1)!(m+1)!(nm)!

=(n+1m+1)

所以有(nm)+(nm+1)=(n+1m+1)

那么给答案前面加一个(x+1x+1),根据刚刚的公式,(x+1x+1) + (x+1x) = (x+2x+1)

(x+2x) + (x+2x+1) = (x+3x+1)

以此类推,最终答案是(x+nn)1

附上代码
#include<bits/stdc++.h>
using namespace std;
#define int long long
#define what_can_I_say main
#define crash_on_you 0
const int N=1e9+100,p=1000003;
int i,j,n,m,ans,t,l,r,x,fact[p+100],ny[p+100];
int op(int a,int b){
	int ans=1;
	while(b){
		if(b&1)ans=ans*a%p;
		a=a*a%p;
		b>>=1;
	}
	return ans;
}
int c(int n,int m){
	if(n<m)return 0;
	return 1ll*(1ll*fact[n]*op(fact[m],p-2)%p)*op(fact[n-m],p-2)%p;
}
int lucas(int n,int m){
	if(n<=p&&m<=p)return c(n,m);
	return c(n%p,m%p)*lucas(n/p,m/p);
}
signed what_can_I_say(){
	fact[0]=1;
	for(i=1;i<p;i++)fact[i]=fact[i-1]*i%p;
	scanf("%lld",&t);
	while(t--){
		scanf("%lld%lld%lld",&n,&l,&r);
		x=r-l+1;
		printf("%lld\n",(lucas(x+n,n)-1+p)%p);
	}
	return crash_on_you;
}

posted @   无敌の暗黑魔王  阅读(40)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· 分享一个免费、快速、无限量使用的满血 DeepSeek R1 模型,支持深度思考和联网搜索!
· 基于 Docker 搭建 FRP 内网穿透开源项目(很简单哒)
· ollama系列1:轻松3步本地部署deepseek,普通电脑可用
· 按钮权限的设计及实现
· Apache Tomcat RCE漏洞复现(CVE-2025-24813)
点击右上角即可分享
微信分享提示