ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

YbtOJ-序列计数【组合数学,莫队】

2022-02-09 16:35:30  阅读:159  来源: 互联网

标签:leq YbtOJ ll id 计数 binom include 莫队 sum


正题


题目大意

求有多少个长度在\(l,r\)之间,值域是\([1,n]\)的严格上升子序列

\(1\leq T,n\leq 10^5,1\leq l\leq r\leq 10^5\)


解题思路

先转换成两个前缀和的差,那么相当于我们要快速求

\[\sum_{i=0}^m\binom{n}{i} \]

的值。

考虑到我们有组合数恒等式\(\binom n m=\binom{n-1}{m-1}+\binom{n-1}{m}\)。

如果我们知道了\(F(n,m)=\sum_{i=0}^m\binom{n}{i}\),那么有\(F(n+1,m)=2F(n,m)-\binom{n}{m}\)(也就是相当于复制一份左移一位相加)。

然后\(F(n,m)\)都可以\(O(1)\)移动\(n,m\)了,直接上莫队。

时间复杂度:\(O(n\sqrt n)\)


code

#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const ll N=5e5+10,P=998244353,T=500;
struct node{
	ll id,r,x;
}a[N];
ll n,m,fac[N],inv[N],ans[N];
ll C(ll n,ll m)
{return fac[n]*inv[m]%P*inv[n-m]%P;}
bool cmp(node x,node y){
	if(x.x/T==y.x/T)
		return x.r<y.r;
	return x.x/T<y.x/T;
}
signed main()
{
	freopen("sequence.in","r",stdin);
	freopen("sequence.out","w",stdout);
	inv[0]=inv[1]=fac[0]=1;
	for(ll i=2;i<N;i++)inv[i]=P-inv[P%i]*(P/i)%P;
	for(ll i=1;i<N;i++)fac[i]=fac[i-1]*i%P,inv[i]=inv[i-1]*inv[i]%P;
	scanf("%lld",&m);
	for(ll i=1;i<=m;i++){
		ll l,r,x;
		scanf("%lld%lld%lld",&l,&r,&x);
		a[++n]=(node){i,r,x};
		a[++n]=(node){-i,l-1,x};
	}
	sort(a+1,a+1+n,cmp);
	ll x=0,r=0,sum=1;
	for(ll i=1;i<=n;i++){
		while(x<a[i].x)sum=(sum*2ll-C(x,r))%P,x++;
		while(x>a[i].x)x--,sum=(P+1)/2*(sum+C(x,r))%P;
		while(r<a[i].r)r++,(sum+=C(x,r))%=P;
		while(r>a[i].r)(sum-=C(x,r))%=P,r--;
		(ans[abs(a[i].id)]+=sum*(a[i].id/abs(a[i].id))%P)%=P;
	}
	for(ll i=1;i<=m;i++)
		printf("%lld\n",(ans[i]+P)%P);
	return 0;
}

标签:leq,YbtOJ,ll,id,计数,binom,include,莫队,sum
来源: https://www.cnblogs.com/QuantAsk/p/15875561.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有