题目链接:BZOJ - 1500
题目分析
我要先说一下,这道题我写了一晚上,然后Debug了一整个白天..........再一次被自己的蒟蒻程度震惊= =
这道题是传说中的Splay维护数列的Boss题目。
前面的几个操作和询问看起来比较正常,就是最后一个维护最大区间和比较复杂。
其实这个也并不是十分复杂,只是要多维护一点东西,事实证明,我代码里的错误都不是与这个询问有关的。
维护每个节点的 Lmx[x] ,即这个节点的子树代表的区间的从左端开始的最大权值和。Rmx[x],同理。Mx[x],这个区间的最大权值和。
然而我的最大权值和都是可以为空的,也就是说如果只能取负数我就什么都不取,就是0。
题目规定的最大权值区间是不能为空的,所以我多维护了一个 Max[x],询问的时候判断如果 Max[x] < 0 ,就输出 Max[x],否则输出我维护的 Mx[x]。
............................然后就是我Debug了一天的东西.........................
经过一天悲剧的探索与尝试,我发现我的错误是出现在了Reverse操作有关的东西。
PushDown写得有非常严重的错误,下面来梳理一下正确的写法,以后要固定一下写法,不能一个代码一个写法...
首先,翻转一个区间的时候,假如这个区间子树的根节点是 x,就Reverse(x)。
之后,PushDown(x) 的时候,判断 Rev[x] 是否为 1,如果是1,就Reverse(Son[x][0]); Reverse(Son[x][1]);
注意,交换左右儿子的操作在Reverse(x)里完成。Rev[x]为1表示的是x的两个儿子还应被 Reverse,但是x的两个儿子的顺序已经是对的。是x的孙子的顺序需要被交换。
代码:
void Reverse(int x)
{
Rev[x] ^= 1;
swap(Son[x][0], Son[x][1]);
} void PushDown(int x)
{
if (Rev[x] == 0) return;
if (Son[x][0]) Reverse(Son[x][0]);
if (Son[x][1]) Reverse(Son[x][1]);
Rev[x] = 0;
}
代码
#include <iostream>
#include <cstdlib>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm> using namespace std; const int MaxN = 1000000 + 5, INF = 999999999; inline void Read(int &Num)
{
char c = getchar();
bool Neg = false;
while (c < '0' || c > '9')
{
if (c == '-') Neg = true;
c = getchar();
}
Num = c - '0'; c = getchar();
while (c >= '0' && c <= '9')
{
Num = Num * 10 + c - '0';
c = getchar();
}
if (Neg) Num = -Num;
} int n, m, Root, Index, Tot1, Tot2, Len;
int V[MaxN], A[MaxN], Father[MaxN], Son[MaxN][2], Sum[MaxN], Rev[MaxN], Stack[MaxN];
int Max[MaxN], Size[MaxN], Lmx[MaxN], Rmx[MaxN], Mx[MaxN], D[MaxN]; bool Rep[MaxN]; char Str[15]; inline int gmax(int a, int b) {return a > b ? a : b;}
inline int gmin(int a, int b) {return a < b ? a : b;} inline void Update(int x)
{
Size[x] = Size[Son[x][0]] + Size[Son[x][1]] + 1;
Sum[x] = Sum[Son[x][0]] + Sum[Son[x][1]] + A[x];
Max[x] = gmax(A[x], gmax(Max[Son[x][0]], Max[Son[x][1]]));
Lmx[x] = gmax(Lmx[Son[x][0]], Sum[Son[x][0]] + A[x] + Lmx[Son[x][1]]);
Rmx[x] = gmax(Rmx[Son[x][1]], Sum[Son[x][1]] + A[x] + Rmx[Son[x][0]]);
Mx[x] = gmax(Rmx[Son[x][0]] + A[x] + Lmx[Son[x][1]], gmax(Rmx[x], Lmx[x]));
Mx[x] = gmax(Mx[x], gmax(Mx[Son[x][0]], Mx[Son[x][1]]));
} inline void Replace(int x, int Num)
{
Rep[x] = true;
D[x] = Num;
A[x] = Num;
Sum[x] = Num * Size[x];
Max[x] = Num;
Lmx[x] = Rmx[x] = Mx[x] = gmax(0, Sum[x]);
} inline void Reverse(int x)
{
Rev[x] ^= 1;
swap(Son[x][0], Son[x][1]);
swap(Lmx[x], Rmx[x]);
} inline void PushDown(int x)
{
if (Rep[x])
{
if (Son[x][0]) Replace(Son[x][0], D[x]);
if (Son[x][1]) Replace(Son[x][1], D[x]);
Rep[x] = false;
}
if (Rev[x])
{
if (Son[x][0]) Reverse(Son[x][0]);
if (Son[x][1]) Reverse(Son[x][1]);
Rev[x] = 0;
}
} int NewNode(int a)
{
int x;
if (Tot2 > 0) x = Stack[Tot2--];
else x = ++Tot1;
Size[x] = 1;
Son[x][0] = Son[x][1] = 0;
Father[x] = 0;
A[x] = Sum[x] = Max[x] = a;
Lmx[x] = Rmx[x] = Mx[x] = gmax(a, 0);
Rev[x] = 0; Rep[x] = false;
return x;
} int Build(int s, int t)
{
int x, m = (s + t) >> 1;
x = NewNode(V[m]);
if (s < m)
{
Son[x][0] = Build(s, m - 1);
Father[Son[x][0]] = x;
}
if (t > m)
{
Son[x][1] = Build(m + 1, t);
Father[Son[x][1]] = x;
}
Update(x);
return x;
} void Rotate(int x, int f)
{
int y = Father[x];
if (y == 0) return;
Son[y][f ^ 1] = Son[x][f];
if (Son[x][f]) Father[Son[x][f]] = y;
Father[x] = Father[y];
if (Father[y])
{
if (y == Son[Father[y]][0]) Son[Father[y]][0] = x;
else Son[Father[y]][1] = x;
}
Son[x][f] = y;
Father[y] = x;
Update(y);
Update(x);
} void Splay(int x, int d)
{
if (x == d) return;
int y;
while (Father[x] != d)
{
y = Father[x];
if (Father[y] == d)
{
if (x == Son[y][0]) Rotate(x, 1);
else Rotate(x, 0);
break;
}
if (y == Son[Father[y]][0])
{
if (x == Son[y][0])
{
Rotate(y, 1);
Rotate(x, 1);
}
else
{
Rotate(x, 0);
Rotate(x, 1);
}
}
else
{
if (x == Son[y][1])
{
Rotate(y, 0);
Rotate(x, 0);
}
else
{
Rotate(x, 1);
Rotate(x, 0);
}
}
}
if (Father[x] == 0) Root = x;
} int Find(int Num)
{
int x = Root, k = Num;
PushDown(x);
while (Size[Son[x][0]] + 1 != k)
{
if (Size[Son[x][0]] + 1 > k)
{
x = Son[x][0];
}
else
{
k -= Size[Son[x][0]] + 1;
x = Son[x][1];
}
PushDown(x);
}
return x;
} void Delete(int x)
{
if (x == 0) return;
if (Son[x][0]) Delete(Son[x][0]);
if (Son[x][1]) Delete(Son[x][1]);
Stack[++Tot2] = x;
} int main()
{
scanf("%d%d", &n, &m);
Tot1 = Tot2 = 0;
Len = n;
for (int i = 1; i <= n; ++i) Read(V[i]);
Max[0] = -INF;
int Temp;
Root = Build(1, n);
Splay(Find(1), 0);
Temp = NewNode(-INF);
Son[Root][0] = Temp;
Father[Temp] = Root;
Update(Root);
Splay(Find(n + 1), 0);
Temp = NewNode(-INF);
Son[Root][1] = Temp;
Father[Temp] = Root;
Update(Root);
int Pos, Tot, New, Num, x, y;
for (int i = 1; i <= m; ++i)
{
scanf("%s", Str);
if (strcmp(Str, "INSERT") == 0)
{
Read(Pos); Read(Tot);
Len += Tot;
for (int j = 1; j <= Tot; ++j) scanf("%d", &V[j]);
New = Build(1, Tot);
x = Find(Pos + 1);
y = Find(Pos + 2);
Splay(x, 0);
Splay(y, x);
Son[y][0] = New;
Father[New] = y;
Update(y);
Update(x);
}
else if (strcmp(Str, "DELETE") == 0)
{
Read(Pos); Read(Tot);
Len -= Tot;
x = Find(Pos);
y = Find(Pos + Tot + 1);
Splay(x, 0);
Splay(y, x);
Delete(Son[y][0]);
Son[y][0] = 0;
Update(y);
Update(x);
}
else if (strcmp(Str, "MAKE-SAME") == 0)
{
Read(Pos); Read(Tot); Read(Num);
x = Find(Pos);
y = Find(Pos + Tot + 1);
Splay(x, 0);
Splay(y, x);
Replace(Son[y][0], Num);
Update(y);
Update(x);
}
else if (strcmp(Str, "REVERSE") == 0)
{
Read(Pos); Read(Tot);
x = Find(Pos);
y = Find(Pos + Tot + 1);
Splay(x, 0);
Splay(y, x);
Reverse(Son[y][0]);
Update(y);
Update(x);
}
else if (strcmp(Str, "GET-SUM") == 0)
{
Read(Pos); Read(Tot);
x = Find(Pos);
y = Find(Pos + Tot + 1);
Splay(x, 0);
Splay(y, x);
printf("%d\n", Sum[Son[y][0]]);
}
else if (strcmp(Str, "MAX-SUM") == 0)
{
x = Find(1);
y = Find(Len + 2);
Splay(x, 0);
Splay(y, x);
if (Max[Son[y][0]] <= 0) printf("%d\n", Max[Son[y][0]]);
else printf("%d\n", Mx[Son[y][0]]);
}
}
return 0;
}
[BZOJ 1500] [NOI2005] 维修数列的更多相关文章
-
bzoj 1500: [NOI2005]维修数列 splay
1500: [NOI2005]维修数列 Time Limit: 10 Sec Memory Limit: 64 MBSubmit: 6556 Solved: 1963[Submit][Status ...
-
BZOJ 1500: [NOI2005]维修数列 (splay tree)
1500: [NOI2005]维修数列 Time Limit: 10 Sec Memory Limit: 64 MBSubmit: 4229 Solved: 1283[Submit][Status ...
-
BZOJ 1500 [NOI2005]维修数列 FHQ Treap
终于A了这题...这题还是很好...但是我太菜...重构了三遍qwq FHQ Treap大法好!qwq...~~ Ins:直接拿输入造一棵树,把原来的树split成[1,pos],[pos+1,n], ...
-
【BZOJ】1500: [NOI2005]维修数列
[算法]splay [题解]数据结构 感谢Occult的模板>_<:HYSBZ 1500 维修数列 #include<cstdio> #include<cctype> ...
-
【BZOJ】1500: [NOI2005]维修数列(splay+变态题)
http://www.lydsy.com/JudgeOnline/problem.php?id=1500 模板不打熟你确定考场上调试得出来? 首先有非常多的坑点...我遇到的第一个就是,如何pushu ...
-
1500: [NOI2005]维修数列
Description Input 输入的第1 行包含两个数N 和M(M ≤20 000),N 表示初始时数列中数的个数,M表示要进行的操作数目.第2行包含N个数字,描述初始时的数列.以下M行,每行一 ...
-
bzoj千题计划221:bzoj1500: [NOI2005]维修数列(fhq treap)
http://www.lydsy.com/JudgeOnline/problem.php?id=1500 1.覆盖标记用INF表示无覆盖标记,要求可能用0覆盖 2.代表空节点的0号节点和首尾的两个虚拟 ...
-
[NOI2005] 维修数列
1500: [NOI2005]维修数列 Time Limit: 10 Sec Memory Limit: 64 MBSubmit: 8397 Solved: 2530 Description In ...
-
[BZOJ1500][NOI2005]维修数列---解题报告
Portal Gun:[BZOJ1500][NOI2005]维修数列 有一段时间没写博客了,最近在刚数据结构......各种板子背得简直要起飞,题目也是一大堆做不完,这里就挑一道平衡树的题来写写好了 ...
随机推荐
-
ASP.Net中无刷新执行Session身份验证
在写一个客户的B/S结构应用程序时,突然发现一个技巧,不知道是否是MS的一个BUG,给相关的有研究的朋友原先考虑写一个检查Session的类,Session失效后,必须转向登陆页面,可每一个调用该类的 ...
-
(整理)streamWriter、streamReader与FileStream
今天偶然使用VS代码分析,发现CA2000警告,然后其中一条为streamWriter.streamReader与FileStream相关内容,特查询并记录一下. 引文地址:http://bbs.cs ...
-
js统计字符串,并且判断出现次数最多的
<!DOCTYPE html> <html lang="en"> <head> <meta charset="UTF-8&quo ...
-
再起航,我的学习笔记之JavaScript设计模式15(组合模式)
组合模式 组合模式(Composite): 又称部分-整体模式,将对象组合成树形结构以表示"部分整体"的层次结构.组合模式使得用户对单个对象和组合对象的使用具有一致性. 如果有一个 ...
-
Python股票分析系列——基础股票数据操作(二).p4
该系列视频已经搬运至bilibili: 点击查看 欢迎来到Python for Finance教程系列的第4部分.在本教程中,我们将基于Adj Close列创建烛台/ OHLC图,这将允许我介绍重新采 ...
-
python 的基础 学习 第五天 基础数据类型的操作方法
1,列表的基本操作方法 1,列表是python中的基础数据类型之一,其他语言中也有类似于列表的数据类型,比如js中叫数组,他是以[ ]括起来,每个元素以逗号隔开,而且他里面可以存放各种数据类型比如: ...
-
k8s搭建问题(1)--OOMKilled
kubectl describe pods ****** --namespace=****** 现象 Host Port: /TCP State: Waiting Reason: CrashLoopB ...
-
微信小程序--地图上添加图片
如何在微信小程序地图添加上,添加图片? 在微信小程序中,地图的层级最高,所以我们没有办法,通过定位,在地图上添加图片等信息; 处理办法: 添加控件:controls; 其中有个属性position,进 ...
-
SQL 时间戳转换为日期
, '1970-01-01 00:00:00') 其中Timestamp为10位的时间戳,+8*3600是获取中国北京时间(东八区)
-
Ubuntu16安装QQ
安装教程: 一:安装依赖库 在终端输入sudo apt-get install libgtk2.0-0:i386 另外,如果是64位系统还要安装ia32-libs 这里我们选择安装lib32ncurs ...