[BZOJ1251]序列终结者
[BZOJ1251]序列终结者
试题描述
网上有许多题,就是给定一个序列,要你支持几种操作:A、B、C、D。一看另一道题,又是一个序列 要支持几种操作:D、C、B、A。尤其是我们这里的某人,出模拟试题,居然还出了一道这样的,真是没技术含量……这样 我也出一道题,我出这一道的目的是为了让大家以后做这种题目有一个“库”可以依靠,没有什么其他的意思。这道题目 就叫序列终结者吧。 【问题描述】 给定一个长度为N的序列,每个序列的元素是一个整数(废话)。要支持以下三种操作: 1. 将[L,R]这个区间内的所有数加上V。 2. 将[L,R]这个区间翻转,比如1 2 3 4变成4 3 2 1。 3. 求[L,R]这个区间中的最大值。 最开始所有元素都是0。
输入
第一行两个整数N,M。M为操作个数。 以下M行,每行最多四个整数,依次为K,L,R,V。K表示是第几种操作,如果不是第1种操作则K后面只有两个数。
输出
对于每个第3种操作,给出正确的回答。
输入示例
4 4 1 1 3 2 1 2 4 -1 2 1 3 3 2 4
输出示例
2
数据规模及约定
N<=50000,M<=100000。
题解
复习了一下 splay,然而不知为何慢的飞起。。。
#include <iostream> #include <cstdio> #include <cstdlib> #include <cstring> #include <cctype> #include <algorithm> using namespace std; int read() { int x = 0, f = 1; char c = getchar(); while(!isdigit(c)){ if(c == '-') f = -1; c = getchar(); } while(isdigit(c)){ x = x * 10 + c - '0'; c = getchar(); } return x * f; } #define maxn 50010 struct Node { int v, siz, mx, addv; bool rev; Node() {} Node(int _): v(_), addv(0), rev(0) {} } ns[maxn]; int rt, ToT, fa[maxn], ch[2][maxn]; void maintain(int o) { ns[o].siz = 1; ns[o].mx = ns[o].v; for(int i = 0; i < 2; i++) if(ch[i][o]) ns[o].mx = max(ns[o].mx, ns[ch[i][o]].mx + ns[ch[i][o]].addv), ns[o].siz += ns[ch[i][o]].siz; return ; } void build(int& o, int l, int r) { if(l > r){ o = 0; return ; } ns[o = ++ToT] = Node(0); if(l == r) return maintain(o); int mid = l + r >> 1; build(ch[0][o], l, mid - 1); build(ch[1][o], mid + 1, r); if(ch[0][o]) fa[ch[0][o]] = o; if(ch[1][o]) fa[ch[1][o]] = o; return maintain(o); } void pushdown(int o) { ns[o].v += ns[o].addv; for(int i = 0; i < 2; i++) if(ch[i][o]) ns[ch[i][o]].addv += ns[o].addv, ns[ch[i][o]].rev ^= ns[o].rev; ns[o].addv = 0; if(ns[o].rev) swap(ch[0][o], ch[1][o]), ns[o].rev = 0; return maintain(o); } void rotate(int u) { int y = fa[u], z = fa[y], l = 0, r = 1; if(z) ch[ch[1][z]==y][z] = u; if(ch[1][y] == u) swap(l, r); fa[u] = z; fa[y] = u; fa[ch[r][u]] = y; ch[l][y] = ch[r][u]; ch[r][u] = y; maintain(y); maintain(u); return ; } int S[maxn], top; void splay(int u) { int t = u; top = 0; S[++top] = t; while(fa[t]) t = fa[t], S[++top] = t; while(top) pushdown(S[top--]); while(fa[u]) { int y = fa[u], z = fa[y]; if(z) { if(ch[0][y] == u ^ ch[0][z] == y) rotate(u); else rotate(y); } rotate(u); } return ; } int split(int u) { if(!u) return 0; splay(u); int tmp = ch[1][u]; fa[tmp] = 0; ch[1][u] = 0; maintain(u); return tmp; } int merge(int a, int b) { if(!a) return maintain(b), b; if(!b) return maintain(a), a; pushdown(a); while(ch[1][a]) a = ch[1][a], pushdown(a); splay(a); ch[1][a] = b; fa[b] = a; maintain(a); return a; } int qkth(int o, int k) { if(!o) return 0; pushdown(o); int ls = ch[0][o] ? ns[ch[0][o]].siz : 0; if(k == ls + 1) return o; if(k > ls + 1) return qkth(ch[1][o], k - ls - 1); return qkth(ch[0][o], k); } int Find(int k) { int u = 1; while(fa[u]) u = fa[u]; return qkth(u, k); } void Add(int ql, int qr, int v) { int lrt = Find(ql - 1), mrt = Find(qr), rrt; split(lrt); rrt = split(mrt); ns[mrt].addv += v; mrt = merge(mrt, rrt); merge(lrt, mrt); return ; } void Rev(int ql, int qr) { int lrt = Find(ql - 1), mrt = Find(qr), rrt; split(lrt); rrt = split(mrt); ns[mrt].rev ^= 1; mrt = merge(mrt, rrt); merge(lrt, mrt); return ; } int Que(int ql, int qr) { int lrt = Find(ql - 1), mrt = Find(qr), rrt; split(lrt); rrt = split(mrt); int ans = ns[mrt].mx + ns[mrt].addv; mrt = merge(mrt, rrt); merge(lrt, mrt); return ans; } int main() { int n = read(), q = read(); build(rt, 1, n); while(q--) { int tp = read(), l = read(), r = read(); if(tp == 1) { int v = read(); Add(l, r, v); } if(tp == 2) Rev(l, r); if(tp == 3) printf("%d\n", Que(l, r)); } return 0; }