2020-2021:teams:legal_string:jxm2001:莫队算法_2

差别

这里会显示出您选择的修订版和当前版本之间的差别。

到此差别页面的链接

--- 2020-2021:teams:legal_string:jxm2001:莫队算法_2 [2021/08/01 19:47]
jxm2001 [算法实现]
+++ 2020-2021:teams:legal_string:jxm2001:莫队算法_2 [2021/08/01 19:50] (当前版本)
jxm2001 [算法例题]
@@ 行 213: / 行 213: @@
 另外注意当跑到另一个块时需要消除之前的右端点影响，可以考虑暴力清空(因为只有 $O(S)$ 个块)，也可以回滚上一个块的右端点。
-==== 算法例题 ====
+==== 例题1 ====
 [[https://loj.ac/p/2874|Loj 2874]]
@@ 行 291: / 行 291: @@
 </hidden>
+==== 例题2 ====
+[[https://ac.nowcoder.com/acm/contest/11256/I|2021牛客暑期多校训练营5 I]]
+=== 题意 ===
+给定一个长度为 $n$ 的序列 $A$，接下来有 $q$ 个询问。每次询问给定 $l,r,k$，求
+$$
+\sum_{i=0}^{k-1}f(l-i,r+i)(n+1)^i
+$$
+其中 $f$ 定义如下
+$$
+f(l,r)=\max(\{k|\exists x\forall i(0\le i\le k-1\to \exists j(l\le j\le r,a_j=x+i))\})
+$$
+=== 题解 ===
+考虑回滚莫队处理询问，难点在于怎么维护最大的 $k$。
+实际上，这等价于维护数轴上的最长连续线段。可以令 $\text{vl}(x)$ 表示以数 $x$ 为左端点的线段在数轴上的右端点。
+$\text{vr}(x)$ 表示以数 $x$ 为右端点的线段在数轴上的左端点。于是只需要保证每条线段的两个端点的 $\text{vl},\text{vr}$ 正确性即可，不难 $O(1)$ 维护。
+时间复杂度 $O\left(n\sqrt q+\sum k\right)$。
+<hidden 查看代码>
+<code cpp>
+const int MAXN=1e5+5,mod=998244353;
+int blk_sz,a[MAXN];
+struct query{
+	int l,r,k,idx;
+	bool operator < (const query &b)const{
+		if(l/blk_sz!=b.l/blk_sz)return l<b.l;
+		return r<b.r;
+	}
+}opt[MAXN];
+struct node{
+	int p1,v1,p2,v2,ans;
+}st[MAXN];
+int curlen,vis[MAXN],vl[MAXN],vr[MAXN],tp;
+void add(int v){
+	if(vis[v]++)return;
+	int l=vis[v-1]?vl[v-1]:v;
+	int r=vis[v+1]?vr[v+1]:v;
+	st[++tp]=node{l,vr[l],r,vl[r],curlen};
+	vr[l]=r;
+	vl[r]=l;
+	curlen=max(curlen,r-l+1);
+}
+void del(int v){
+	if(--vis[v])return;
+	vr[st[tp].p1]=st[tp].v1;
+	vl[st[tp].p2]=st[tp].v2;
+	curlen=st[tp].ans;
+	tp--;
+}
+int solve(int l,int r,int k,int n){
+	int ans=curlen,base=1;
+	_for(i,1,k){
+		add(a[l-i]);
+		add(a[r+i]);
+		base=1LL*base*(n+1)%mod;
+		ans=(ans+1LL*curlen*base)%mod;
+	}
+	for(int i=k-1;i;i--){
+		del(a[r+i]);
+		del(a[l-i]);
+	}
+	return ans;
+}
+int ans[MAXN];
+int main()
+{
+	int n=read_int(),q=read_int();
+	_rep(i,1,n)
+	a[i]=read_int();
+	_for(i,0,q){
+		int l=read_int(),r=read_int(),k=read_int();
+		opt[i]=query{l,r,k,i};
+	}
+	blk_sz=n/sqrt(q)+1;
+	sort(opt,opt+q);
+	_for(i,0,q){
+		if(opt[i].l/blk_sz==opt[i].r/blk_sz){
+			_rep(j,opt[i].l,opt[i].r)
+			add(a[j]);
+			ans[opt[i].idx]=solve(opt[i].l,opt[i].r,opt[i].k,n);
+			for(int j=opt[i].r;j>=opt[i].l;j--)
+			del(a[j]);
+		}
+	}
+	int lef=1,rig=0,lst=-1;
+	_for(i,0,q){
+		if(opt[i].l/blk_sz!=opt[i].r/blk_sz){
+			if(opt[i].l/blk_sz!=lst){
+				while(lef<=rig)del(a[rig--]);
+				lst=opt[i].l/blk_sz;
+				rig=min(lst*blk_sz+blk_sz-1,n);
+				lef=rig+1;
+			}
+			while(rig<opt[i].r)add(a[++rig]);
+			int tlef=lef;
+			while(tlef>opt[i].l)add(a[--tlef]);
+			ans[opt[i].idx]=solve(opt[i].l,opt[i].r,opt[i].k,n);
+			while(tlef<lef)del(a[tlef++]);
+		}
+	}
+	_for(i,0,q)enter(ans[i]);
+	return 0;
+}
+</code>
+</hidden>
 ===== $\text{bitset}$ 与莫队 =====

2020-2021/teams/legal_string/jxm2001/莫队算法_2.1627818460.txt.gz · 最后更改: 2021/08/01 19:47 由 jxm2001