題目鏈接：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; }

　　

轉載于:https://www.cnblogs.com/JoeFan/p/4336035.html

總結

以上是生活随笔為你收集整理的[BZOJ 1500] [NOI2005] 维修数列的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。