点此看题面

大致题意: 有单点修改数字和区间着色两种修改操作,询问你某段区间内包含所有颜色且数字和最小的子区间的数字和,或某段区间内没有重复颜色且数字和最大的子区间的数字和。数据随机。

\(ODT\)维护颜色

看到区间着色且题目中强调数据随机,容易想到使用\(ODT\)去求解。

于是,我们就可以考虑用\(ODT\)来对颜色进行维护。

线段树维护数字和

但是考虑到要求区间数字和,\(ODT\)就很难搞了。

考虑到我们其实可以对于每个询问,每次找到一个合法区间再询问数字和更新答案。

也就是说,可以把维护颜色和维护数字和分开。

那维护数字和自然可以用线段树喽。

其实本来用树状数组似乎更优,但这题由于写法问题可能还会需要求区间\(Min\)或\(Max\),树状数组搞起来就很吃力了。

处理询问

以上大致介绍了,可以用线段树维护数字和,\(ODT\)维护颜色。

接下来,再具体介绍一下如何去处理询问。

考虑到这两种询问实际上都具有单调性,因此可以直接用双指针搞。

我们枚举右端点,然后移动左端点。

  • 对于第一种询问,我们需要在保证颜色数量等于\(c\)的情况下才能移动左端点,然后在移动的同时更新答案(不然统计答案的过程会略显麻烦)。

    而\(c=1\)的情况可能要特判,直接输出这段区间的最小值。

  • 对于第二种询问,我们只需在出现不合法情况,即出现重复颜色时移动左端点。考虑我们此时刚把右端点的颜色加入,若出现重复颜色,必然是由右端点导致的,因此不断移动左端点直至右端点所属颜色出现次数为\(1\)即可。

    但注意,在处理第二种询问时,一旦出现右端点\(Size>1\)的情况,我们在统计完其答案后,需将左端点移动到右端点的位置上,不然就会出现重复颜色。

    同理,在统计第二种询问的答案时,要求左端点的右边界右端点的左边界之间的区间和,不然同样会出现重复颜色。

    而且,在第二种情况时容易漏考虑单个颜色的贡献,因此要初始化\(res\)为整段区间的\(Max\)。

大致就是这些吧。

代码

#include<bits/stdc++.h>
#define Tp template<typename Ty>
#define Ts template<typename Ty,typename... Ar>
#define Reg register
#define RI Reg int
#define Con const
#define CI Con int&
#define I inline
#define W while
#define N 100000
#define INF 1e9
#define Gmax(x,y) (x<(y)&&(x=(y)))
#define Gmin(x,y) (x>(y)&&(x=(y)))
using namespace std;
int n,c,a[N+5];
class FastIO
{
private:
#define FS 100000
#define tc() (A==B&&(B=(A=FI)+fread(FI,1,FS,stdin),A==B)?EOF:*A++)
#define pc(c) (C^FS?FO[C++]=c:(fwrite(FO,1,C,stdout),FO[(C=0)++]=c))
#define tn (x<<3)+(x<<1)
#define D isdigit(c=tc())
int T,C;char c,*A,*B,FI[FS],FO[FS],S[FS];
public:
I FastIO() {A=B=FI;}
Tp I void read(Ty& x) {x=0;W(!D);W(x=tn+(c&15),D);}
Tp I void write(Ty x) {x<0&&(pc('-'),x=-x);W(S[++T]=x%10+48,x/=10);W(T) pc(S[T--]);}
Ts I void read(Ty& x,Ar&... y) {read(x),read(y...);}
Tp I void writeln(Con Ty& x) {write(x),pc('\n');}
I void clear() {fwrite(FO,1,C,stdout),C=0;}
}F;
class SegmentTree//线段树
{
private:
#define STO l,hl,rt<<1
#define ORZ hl+1,r,rt<<1|1
#define PU(x) (O[x]=O[x<<1]+O[x<<1|1])
int n,v[N+5];
struct Interval//维护区间信息
{
int S,Mx,Mn;I Interval(CI s=0,CI mx=-INF,CI mn=INF):S(s),Mx(mx),Mn(mn){}
I Interval operator + (Con Interval& t) Con {return Interval(S+t.S,max(Mx,t.Mx),min(Mn,t.Mn));}
}O[N<<2];
I void Build(CI l,CI r,CI rt)//建树
{
if(!(l^r)) return (void)(O[rt]=Interval(v[l],v[l],v[l]));
RI hl=l+r>>1;Build(STO),Build(ORZ),PU(rt);
}
I int qsu(CI l,CI r,CI rt,CI ql,CI qr)//询问区间和
{
if(ql<=l&&r<=qr) return O[rt].S;RI hl=l+r>>1;
return (ql<=hl?qsu(STO,ql,qr):0)+(qr>hl?qsu(ORZ,ql,qr):0);
}
I int qmx(CI l,CI r,CI rt,CI ql,CI qr)//询问区间最大值
{
if(ql<=l&&r<=qr) return O[rt].Mx;RI hl=l+r>>1,t,res=-INF;
return ql<=hl&&(t=qmx(STO,ql,qr),Gmax(res,t)),qr>hl&&(t=qmx(ORZ,ql,qr),Gmax(res,t)),res;
}
I int qmn(CI l,CI r,CI rt,CI ql,CI qr)//询问区间最小值
{
if(ql<=l&&r<=qr) return O[rt].Mn;RI hl=l+r>>1,t,res=INF;
return ql<=hl&&(t=qmn(STO,ql,qr),Gmin(res,t)),qr>hl&&(t=qmn(ORZ,ql,qr),Gmin(res,t)),res;
}
public:
I void Init(CI x,int* s) {for(RI i=1;i<=x;++i) v[i]=s[i];Build(1,n=x,1);}
I void Update(CI x,CI y)//单点修改(这里使用非递归写法)
{
RI l=1,r=n,rt=1,hl;W(l^r) hl=l+r>>1,x<=hl?(r=hl,rt<<=1):(l=hl+1,(rt<<=1)|=1);
O[rt]=Interval(y,y,y);W(rt>>=1) PU(rt);
}
I int QSum(CI x,CI y) {return qsu(1,n,1,x,y);}
I int QMax(CI x,CI y) {return qmx(1,n,1,x,y);}
I int QMin(CI x,CI y) {return qmn(1,n,1,x,y);}
}T;
class ODT
{
private:
#define IT set<Il>::iterator
#define ins insert
#define era erase
#define fir first
#define LB lower_bound
#define add(x) (!cnt[x]++&&++tot)
#define del(x) (!--cnt[x]&&--tot)
int cnt[N+5];
struct Il//维护区间信息
{
int l,r,v;I Il(CI x=0,CI y=0,CI p=0):l(x),r(y),v(p){}
I bool operator < (Con Il& t) Con {return l<t.l;}
};set<Il> S;
I IT Sp(CI x)//分裂操作
{
IT t;if((t=S.LB(Il(x)))!=S.end()&&!(t->l^x)) return t;
RI l=(--t)->l,r=t->r,v=t->v;S.era(t),S.ins(Il(l,x-1,v));
return S.ins(Il(x,r,v)).fir;
}
public:
I void Init(CI x,int* s)//初始化节点信息
{
for(RI i=(s[0]=s[x+1]=-1,1),t=0;i<=x+2;++i)
s[i]^s[i-1]&&(S.insert(Il(t,i-1,s[i-1])),t=i);
}
I void Assign(CI x,CI y,CI v)//推平操作
{
IT tr=Sp(y+1),tl=Sp(x);S.era(tl,tr),S.ins(Il(x,y,v));
}
I int Q1(CI x,CI y)//处理第一种询问
{
if(!(c^1)) return T.QMin(x,y);//特判c=1
memset(cnt,0,sizeof(cnt));IT tr=Sp(y+1),tl=Sp(x),ti=tl;RI t,res=INF,tot=0;//初始化
W(ti!=tr) {add(ti->v);W(!(tot^c)) t=T.QSum(tl->r,ti->l),Gmin(res,t),del((tl++)->v);++ti;}//移动右端点,在保证颜色数量等于c的情况下才能移动左端点,然后在移动的同时更新答案
return res==INF?-1:res;//判断是否无解,无解返回-1
}
I int Q2(CI x,CI y)//处理第二种询问
{
memset(cnt,0,sizeof(cnt));IT tr=Sp(y+1),tl=Sp(x),ti=tl;RI t,res=T.QMax(x,y);//初始化
W(ti!=tr)
{
++cnt[ti->v];W(ti!=tl&&cnt[ti->v]>1) --cnt[(tl++)->v];//不断移动左端点直至右端点所属颜色出现次数为1
ti!=tl&&(t=T.QSum(tl->r,ti->l),Gmax(res,t));//更新答案
if(ti->l^ti->r) W(ti!=tl) --cnt[(tl++)->v];++ti;//当出现右端点Size>1的情况,将左端点移动到右端点的位置上
}return res;//返回答案
}
}O;
int main()
{
RI Qtot,i,op,x,y,z;for(F.read(n,Qtot,c),i=1;i<=n;++i) F.read(a[i]);T.Init(n,a);//读入数据
for(i=1;i<=n;++i) F.read(a[i]);O.Init(n,a);W(Qtot--)
{
switch(F.read(op,x,y),op)//处理操作
{
case 1:T.Update(x,y);break;case 2:F.read(z),O.Assign(x,y,z);break;
case 3:F.writeln(O.Q1(x,y));break;case 4:F.writeln(O.Q2(x,y));break;
}
}return F.clear(),0;
}

【洛谷5251】[LnOI2019] 第二代图灵机(线段树+ODT)的更多相关文章

  1. 【题解】Luogu P5251 [LnOI2019]第二代图灵机

    原题传送门 前置芝士:珂朵莉树 珂朵莉树的主要功能是区间赋值 这道题还算明显(操作2) 一开始看见这题觉得很毒瘤,但仔细想想发现颜色和数字之间没有什么关系 我们一共要维护三个东西: 1.区间和:树状数 ...

  2. 洛谷P3928 Sequence2(dp,线段树)

    题目链接: 洛谷 题目大意在描述底下有.此处不赘述. 明显是个类似于LIS的dp. 令 $dp[i][j]$ 表示: $j=1$ 时表示已经处理了 $i$ 个数,上一个选的数来自序列 $A[0]$ 的 ...

  3. 洛谷P3434 [POI2006]KRA-The Disks(线段树)

    洛谷题目传送门 \(O(n)\)的正解算法对我这个小蒟蒻真的还有点思维难度.洛谷题解里都讲得很好. 考试的时候一看到300000就直接去想各种带log的做法了,反正不怕T...... 我永远只会有最直 ...

  4. 洛谷.3733.[HAOI2017]八纵八横(线性基 线段树分治 bitset)

    LOJ 洛谷 最基本的思路同BZOJ2115 Xor,将图中所有环的异或和插入线性基,求一下线性基中数的异或最大值. 用bitset优化一下,暴力的复杂度是\(O(\frac{qmL^2}{w})\) ...

  5. 【题解】洛谷P4145 花神游历各国(线段树)

    洛谷P4145:https://www.luogu.org/problemnew/show/P4145 思路 这道题的重点在于sqrt(1)=1 一个限制条件 与正常线段树不同的是区间修改为开方 那么 ...

  6. 洛谷P4588 [TJOI2018]数学计算 【线段树】

    题目链接 洛谷P4588 题解 用线段树维护即可 #include<algorithm> #include<iostream> #include<cstring> ...

  7. 【洛谷3822】[NOI2017] 整数(线段树压位)

    题目: 洛谷 3822 分析: 直接按题意模拟,完了. 将每次加 / 减拆成不超过 \(32\) 个对单独一位的加 / 减. 考虑给一个二进制位(下称「当前位」)加 \(1\) 时,如果这一位本来就是 ...

  8. 洛谷P4556 雨天的尾巴(线段树合并)

    洛谷P4556 雨天的尾巴 题目链接 题解: 因为一个点可能存放多种物品,直接开二维数组进行统计时间.空间复杂度都不能承受.因为每一个点所拥有的物品只与其子树中的点有关,所以可以考虑对每一个点来建立一 ...

  9. 洛谷P1712 [NOI2016]区间 尺取法+线段树+离散化

    洛谷P1712 [NOI2016]区间 noi2016第一题(大概是签到题吧,可我还是不会) 链接在这里 题面可以看链接: 先看题意 这么大的l,r,先来个离散化 很容易,我们可以想到一个结论 假设一 ...

随机推荐

  1. Linux用户登录信息

    1.用户登录日志信息 /var/run/utmp:记录当前正在登录系统的用户信息,默认由who和w记录当前登录用户的信息,uptime记录系统启动时间: /var/log/wtmp:记录当前正在登录和 ...

  2. vue.js请求数据(axios)

    使用npm安装axios npm install axios --save 在main.js中引入axios import axios from "axios"; 注册axios到 ...

  3. Win10家庭版组策略gpedit.msc的问题

    大家都认为,Windows家庭版中并不包含组策略,其实不然,它是有相关文件的,只是不让你使用而已.那么我们让系统允许你使用就好了.首先你需要在桌面上新建一个txt文本文档.然后将以下代码复制到这个新建 ...

  4. pl/sql 如何配置连接远程一个或多个数据库

    参考链接 https://blog.csdn.net/yy_love_my/article/details/45720277

  5. 2019.03.22 读书笔记 Linq中的IEnumerable与IQueryable

    Linq主要分为 Linq to object .Linq to XML .Linq to sql(现在都用EF了). 针对object和sql,微软分别设计了两套接口IEnumerable.IQue ...

  6. RTT之AT命令组件

    包含客户端和服务器:用于GPRS和3G的通讯命令格式.四种基本功能 测试功能:AT+<x>=? 用于查询命令参数格式及取值范围: 查询功能:AT+<x>? 用于返回命令参数当前 ...

  7. Choose and divide(唯一分解定理)

    首先说一下什么是唯一分解定理 唯一分解定理:任何一个大于1的自然数N,如果N不是质数,那么N可以分解成有限个素数的乘积:例:N=(p1^a1)*(p2^a2)*(p3^a3)......其中p1< ...

  8. Perf -- Linux下的系统性能调优工具,第 2 部分 特殊用法以及内核调优示例

    https://www.ibm.com/developerworks/cn/linux/l-cn-perf2/ 本文内容简介 之前介绍了 perf 最常见的一些用法,关注于 Linux 系统上应用程序 ...

  9. [OpenStack] [Liberty] Neutron单网卡桥接模式访问外网

    环境配置: * Exsi一台 * Exsi创建的单网卡虚拟机一台 * Ubuntu 14LTS 64位操作系统 * OpenStack Liberty版本 * 使用Neutron网络而非Nova网络 ...

  10. Neutron命令测试1

    Refer: http://wenku.baidu.com/link?url=DtrbhO0A393hg8kOWKX0XYuZtSC8Iu0occn8NF1pYcUwNzlaSq5qXCQoNEBDM ...