1 条题解
-
0
刚放这题时有一篇题解(是谁已经记不清了),我就是用的他的做法做出的,但是现在被删了。所以我按自己的理解再发一遍。
给定个操作,每个操作形如 表示将 ,或者是 表示求 的值。
现有 个询问,每个询问给出 ,表示求将数列初始置为 后依次执行内的操作,求所有 操作的答案和。
同阶,,6s。
与官方题解做法不同。 注意区分“查询”和“询问“,查询指2操作。
对操作扫描线,想方设法维护出 时所有 的答案。这样不仅去掉了 的限制,而且只需要求屏蔽掉 的操作后的答案,模型简化。
考虑对序列分块,分别计算以下贡献:
- 修改对散块查询的贡献
- 整块修改对整块查询的贡献
- 散块修改对整块查询的贡献
对于情况1,枚举每个位置的查询,我们可以轻松维护出每个位置最新覆盖的数(整块标记和散块标记),并将这些操作更新进答案中。
具体地,假设当前位操作的时刻是 ,那么如果 那么会对询问造成贡献。注意到这里的 (修改)共有 个,而询问只有 个,所以使用 分块,这部分贡献就被处理完了。
对于情况2,如果当前整块的值整个相同(即上一次是整块修改),实际上就和单个数是一样的,所以同样套用刚刚的数据结构做即可。
对于情况3,实际上我们求的是当前查询整块的值不整个相同时所有查询带来的贡献。注意到部分修改一个块只能是散块修改,并且如果整块修改覆盖了原来的散块修改也要去掉散块修改的贡献。易知这里总的操作数是 (上界 )。
单独考虑所有的块并枚举每个操作。如果当前操作有对当前块的散块修改或者删除散块修改 的修改,就暴力执行。
考虑这里一个散块修改操作的贡献:设其执行时间为 ,那么对于 的询问会造成贡献。贡献就是权值 加入时刻到询问时刻中当前块查询的次数 删除时刻到当前时刻中当前块查询的次数。实际上删除和加入本质相同,只是权值互为相反数,之后不再分类讨论。
对于所谓的加入时刻到询问时刻中当前块查询的次数,可以转化成初始时刻到询问时刻的查询次数-初始时刻到加入时刻的查询次数。
考虑维护两个 分块,分别维护初始时刻到当前所有权值的贡献和以及所有权值乘上当前查询次数的贡献和。
因此总复杂度是。本题可以轻松通过。做完可以去试试加强版,时限少了1s,非常卡常。
情况三比较困难,给一下代码:
struct block2 {//fast query ll s1[N],s2[505],s3[N],s4[505]; V clear() {fo(i,1,m)s1[i]=s3[i]=0;fo(i,0,m>>10)s2[i]=s4[i]=0;} V add(int x,ll y,ll z) { if(!x||!y&&!z) return; int lim=(x>>10)-1,lim2=(x>>10)<<10; fo(i,0,lim) s2[i]+=y,s4[i]+=z; fo(i,lim2,x) s1[i]+=y,s3[i]+=z; } L ques(int x,ll sum) { return sum*(s3[x]+s4[x>>10])-(s1[x]+s2[x>>10]); } }T2; V Modify(int l,int r,int k) { rub[0]=0;int Len=r-l+1,X=bel[l]; fo(i,l,r) { if(s1[i]) { if(!bj[s1[i]]) rub[++rub[0]]=s1[i]; bj[s1[i]]++; } s1[i]=k; } ll t; fo(i,1,rub[0]) {//这里这么写是因为有一个颜色段均摊的复杂度分析,由于总的操作次数是线性的所以颜色段数也是线性 t=1ll*a[rub[i]].k*bj[rub[i]]; T2.add(rub[i],-t*sum,-t);bj[rub[i]]=0; } if(k) T2.add(k,1ll*sum*a[k].k*Len,1ll*a[k].k*Len); if(Vis&&Cov) { T2.add(Cov,1ll*sum*a[Cov].k*(len[X]-Len),1ll*a[Cov].k*(len[X]-Len)); fo(i,le[X],l-1) s1[i]=Cov; fo(i,r+1,ri[X]) s1[i]=Cov; } Vis=0; } V Work() { memset(s1,0,sizeof(s1)); int X,Y,T; fo(x,1,bel[n]) { Vis=1,sum=0,Cov=0; T2.clear(); T=1; fo(i,1,m) { X=bel[a[i].l];Y=bel[a[i].r]; if(a[i].k==-1) { if(X<x&&x<Y) sum++; } else { if(X==Y) { if(X==x) Modify(a[i].l,a[i].r,i); } else { if(X==x) Modify(a[i].l,ri[x],i); if(Y==x) Modify(le[x],a[i].r,i); if(X<x&&x<Y) { if(!Vis) Modify(le[x],ri[x],0),Vis=1; Cov=i; } } } while(q[T].r==i) ans[q[T].id]+=T2.ques(q[T].l,sum),++T; } } }
- 1
信息
- ID
- 7497
- 时间
- 5500ms
- 内存
- 512MiB
- 难度
- 10
- 标签
- 递交数
- 23
- 已通过
- 1
- 上传者