伸展数最基本操作的模板,区间求和,区间更新。为了方便理解,特定附上一自己搞的搓图

这是样例中的数据输入后建成的树,其中的1,2是加入的边界顶点,数字代表节点编号,我们如果要对一段区间[l, r]进行操作,只需要把第l-1位的数旋转到0节点下面,把r+1位的数旋转到当前的root下面,就如上图所示,那么椭圆里表示的就是区间[l, r]。

附上注释代码。指针版本的比静态数组的快1s多。。

/* **********************************************
Author : JayYe
Created Time: 2013-8-16 11:14:36
File Name : zzz.cpp
*********************************************** */ #include <stdio.h>
#include <string.h>
#include <algorithm>
using namespace std; #define LL __int64
#define keytree (ch[ ch[root][1] ][0])
#define lson ch[x][0]
#define rson ch[x][1] const int maxn = 111111; struct Splaytree { int pre[maxn], ch[maxn][2], sz[maxn];
int root, top;
// 旋转操作, c = 0代表左旋, c = 1代表右旋
void Rotate(int x, int c) {
int y = pre[x];
push_down(y), push_down(x);
ch[y][!c] = ch[x][c];
pre[ch[x][c]] = y;
if(pre[y]) ch[pre[y]][ ch[pre[y]][1] == y ] = x;
pre[x] = pre[y];
ch[x][c] = y;
pre[y] = x;
push_up(y);
}
// Splay 操作, 把 x 节点转到go的下面
void Splay(int x, int go) {
while(pre[x] != go) {
if(pre[pre[x]] == go) {
Rotate(x, ch[pre[x]][0] == x);
}
else {
int y = pre[x] , z = pre[y];
int f = (ch[z][1] == y);
if(ch[y][f] == x)
Rotate(y, !f), Rotate(x, !f); // 一字型旋转
else
Rotate(x, f), Rotate(x, !f); // 之字型旋转
}
}
push_up(x);
if(go == 0) root = x;
}
// 将第k位的数旋转到go的下面
void RotateTo(int k, int go) {
int x = root;
push_down(root);
while(sz[lson] != k) {
if(k < sz[lson]) {
x = lson;
}
else {
k -= sz[lson] + 1;
x = rson;
}
push_down(x);
}
Splay(x, go);
} void debug() {printf("%d\n",root);Treaval(root);}
void Treaval(int x) {
if(x) {
Treaval(ch[x][0]);
printf("结点%2d:左儿子 %2d 右儿子 %2d 父结点 %2d size = %2d ,val = %2d\n",x,ch[x][0],ch[x][1],pre[x],sz[x],val[x]);
Treaval(ch[x][1]);
}
} int val[maxn], add[maxn], a[maxn];
LL sum[maxn];
// 把儿子节点的信息更新上来
void push_up(int x) {
sz[x] = sz[lson] + sz[rson] + 1;
sum[x] = val[x] + add[x] + sum[lson] + sum[rson];
}
// 标记下传
void push_down(int x) {
if(add[x]) {
val[x] += add[x];
add[lson] += add[x];
add[rson] += add[x];
sum[lson] += (LL)add[x]*sz[lson];
sum[rson] += (LL)add[x]*sz[rson];
add[x] = 0;
}
} void newnode(int &x, int c) {
x = ++top;
lson = rson = 0;
sz[x] = 1;
val[x] = sum[x] = c;
add[x] = 0;
} void build(int &x, int l, int r, int f) {
if(l > r) return ;
int mid = (l+r)/2;
newnode(x, a[mid]);
build(lson, l, mid-1, x);
build(rson, mid+1, r, x);
pre[x] = f;
push_up(x);
} void init(int n) {
ch[0][0] = ch[0][1] = pre[0] = 0;
top = root = 0;
newnode(root, -1);
newnode(ch[root][1], -1);
// 为了方便处理边界,加两个边界顶点
pre[top] = root;
sz[root] = 2; for(int i = 0;i < n; i++) scanf("%d", &a[i]);
build(keytree, 0, n-1, ch[root][1]);
push_up(ch[root][1]);
push_up(root);
} void update(int l ,int r, int c) {
RotateTo(l-1, 0);
RotateTo(r+1, root);
add[keytree] += c;
sum[keytree] += (LL)c*sz[keytree];
} LL query(int l, int r) {
RotateTo(l-1, 0);
RotateTo(r+1, root);
return sum[keytree];
} }spt; int main() {
int n, m, l, r, c;
char s[2];
scanf("%d%d", &n, &m);
spt.init(n);
while(m--) {
scanf("%s%d%d", s, &l, &r);
if(s[0] == 'Q')
printf("%I64d\n", spt.query(l, r));
else {
scanf("%d", &c);
spt.update(l, r, c);
}
}
return 0;
}
/* **********************************************
Author : JayYe
Created Time: 2013-8-16 15:19:38
File Name : zzz.cpp
*********************************************** */ #include <stdio.h>
#include <string.h>
#include <algorithm>
using namespace std; #define LL __int64
#define keytree (root->ch[1]->ch[0])
#define lson x->ch[0]
#define rson x->ch[1] const int maxn = 111111; struct NODE {
struct NODE *ch[2], *pre;
int add, val, sz, id;
LL sum;
void push_down() {
if(add) {
val += add;
if(ch[0]) {
ch[0]->add += add;
ch[0]->sum += (LL)add*ch[0]->sz;
}
if(ch[1]) {
ch[1]->add += add;
ch[1]->sum += (LL)add*ch[1]->sz;
}
add = 0;
}
} void push_up() {
sz = ch[0]->sz + ch[1]->sz + 1;
sum = val + add + ch[0]->sum + ch[1]->sum;
} }node[maxn], *null = &node[0], *root; struct Splaytree { int top;
// 旋转操作, c = 0代表左旋, c = 1代表右旋
void Rotate(NODE *x, int c) {
NODE *y = x->pre;
y->push_down(), x->push_down();
y->ch[!c] = x->ch[c];
x->ch[c]->pre = y;
x->pre = y->pre;
if(y->pre != NULL) y->pre->ch[ y->pre->ch[1] == y] = x;
x->ch[c] = y;
y->pre = x;
y->push_up();
}
// Splay 操作, 把 x 节点转到go的下面
void Splay(NODE *x, NODE *go) {
while(x->pre != go) {
if(x->pre->pre == go) {
Rotate(x, x->pre->ch[0] == x);
}
else {
NODE *y = x->pre, *z = y->pre;
int f = (z->ch[1] == y);
if(y->ch[f] == x)
Rotate(y, !f) , Rotate(x, !f); // 一字型旋转
else
Rotate(x, f) , Rotate(x, !f); // 之字型旋转
}
}
x->push_up();
if(go == null) root = x;
}
// 将第k位的数旋转到go的下面
void RotateTo(int k, NODE *go) {
NODE *x = root;
x->push_down();
while(lson->sz != k) {
if(lson->sz > k) {
x = lson;
}
else {
k -= lson->sz + 1;
x = rson;
}
x->push_down();
}
Splay(x, go);
} void debug(NODE *x) {
if(x != null) {
printf("节点: %2d 左儿子: %2d 右儿子: %2d size = %2d val = %2d\n",
x->id, x->ch[0]->id, x->ch[1]->id, x->sz, x->val);
debug(x->ch[0]);
debug(x->ch[1]);
}
} int a[maxn]; NODE *newnode(NODE* f, int c) {
NODE *x = &node[++top];
x->id = top;
x->val = x->sum = c;
x->ch[0] = x->ch[1] = null;
x->sz = 1;
x->add = 0;
x->pre = f;
return x;
} void build(NODE* &x, int l, int r, NODE *f) {
if(l > r) return ;
int mid = (l+r)/2;
x = newnode(f, a[mid]);
build(lson, l, mid-1, x);
build(rson, mid+1, r, x);
x->push_up();
} void init(int n) {
null->id = 0;
null->ch[0] = null->ch[1] = NULL;
null->sz = null->add = null->sum = null->val = 0;
// null->pre = NULL;
top = 0;
root = newnode(null, -1);
root->ch[1] = newnode(root, -1); for(int i = 0;i < n; i++) scanf("%d", &a[i]);
build(keytree, 0, n-1, root->ch[1]);
root->ch[1]->push_up(); root->push_up();
} void update() {
int l, r, c;
scanf("%d%d%d", &l, &r, &c);
RotateTo(l-1, null);
RotateTo(r+1, root);
keytree->add += c;
keytree->sum += (LL)c*keytree->sz;
} void query() {
int l, r;
scanf("%d%d", &l, &r);
RotateTo(l-1, null);
RotateTo(r+1, root);
printf("%I64d\n", keytree->sum);
}
}spt; int main() {
int n, m;
char s[2];
scanf("%d%d", &n, &m);
spt.init(n);
while(m--) {
scanf("%s", s);
if(s[0] == 'Q')
spt.query();
else
spt.update();
}
return 0;
}

POJ 3468 A Simple Problem with Integers (伸展树区间更新求和操作 , 模板)的更多相关文章

  1. [POJ] 3468 A Simple Problem with Integers [线段树区间更新求和]

    A Simple Problem with Integers   Description You have N integers, A1, A2, ... , AN. You need to deal ...

  2. poj 3468 A Simple Problem with Integers (线段树区间更新求和lazy思想)

    A Simple Problem with Integers Time Limit: 5000MS   Memory Limit: 131072K Total Submissions: 75541   ...

  3. poj 3468 A Simple Problem with Integers 线段树区间更新

    id=3468">点击打开链接题目链接 A Simple Problem with Integers Time Limit: 5000MS   Memory Limit: 131072 ...

  4. POJ 3468 A Simple Problem with Integers(线段树,区间更新,区间求和)

    A Simple Problem with Integers Time Limit: 5000MS   Memory Limit: 131072K Total Submissions: 67511   ...

  5. (简单) POJ 3468 A Simple Problem with Integers , 线段树+区间更新。

    Description You have N integers, A1, A2, ... , AN. You need to deal with two kinds of operations. On ...

  6. POJ 3468 A Simple Problem with Integers(线段树区间更新)

    题目地址:POJ 3468 打了个篮球回来果然神经有点冲动. . 无脑的狂交了8次WA..竟然是更新的时候把r-l写成了l-r... 这题就是区间更新裸题. 区间更新就是加一个lazy标记,延迟标记, ...

  7. POJ 3468 A Simple Problem with Integers(线段树区间更新,模板题,求区间和)

    #include <iostream> #include <stdio.h> #include <string.h> #define lson rt<< ...

  8. POJ 3468 A Simple Problem with Integers 线段树 区间更新

    #include<iostream> #include<string> #include<algorithm> #include<cstdlib> #i ...

  9. poj 3468 A Simple Problem with Integers 线段树区间加,区间查询和

    A Simple Problem with Integers Time Limit: 1 Sec  Memory Limit: 256 MB 题目连接 http://poj.org/problem?i ...

随机推荐

  1. Netsharp快速入门(之15) 销售管理(报表B 销售季度表)

    作者:秋时 杨昶   转载须说明出处 4.6.2  销售季度表(交叉表) 1.1.1.1 交叉表带数据源和不带数据源区别 带数据源的可以自定义数据源,可以从实体,也可以从Sql脚本得到数据源,并能自定 ...

  2. 设计模式之状态模式(State)

    状态模式原理:随着状态的变化,对象的行为也发生变化 代码如下: #include <iostream> #include <string> #include <list& ...

  3. BZOJ 3143 HNOI2013 游走 高斯消元 期望

    这道题是我第一次使用高斯消元解决期望类的问题,首发A了,感觉爽爽的.... 不过笔者在做完后发现了一些问题,在原文的后面进行了说明. 中文题目,就不翻大意了,直接给原题: 一个无向连通图,顶点从1编号 ...

  4. if in hlsl

    seems that in HLSL_4, we can use if https://msdn.microsoft.com/en-us/library/bb313972(v=xnagamestudi ...

  5. AngularJs学习笔记--Guide教程系列文章索引

    在很久很久以前,一位前辈向我推荐AngularJs.但当时我没有好好学习,仅仅是讲文档浏览了一次.后来觉醒了……于是下定决心好好理解这系列的文档,并意译出来(英文水平不足……不能说是翻译,有些实在是看 ...

  6. python 二分法查找

    这个也是之前写的程序,现在把它贴上来 #!/usr/bin/python import os os.system('clear') def binsearch(seq,x,low,high): mid ...

  7. Unity3D判断鼠标向右或向左滑动,响应不同的事件

    private var first = Vector2.zero; private var second = Vector2.zero; function Update () { } function ...

  8. NWR协议

    NWR是一种在分布式存储系统中用于控制一致性级别的一种策略.在Amazon的Dynamo云存储系统中,就应用NWR来控制一致性. 让我们先来看看这三个字母的含义:N:在分布式存储系统中,有多少份备份数 ...

  9. 使用 polyfills 的简易方法

    本文作者为 Andrew Betts 与 Robert Nyman.Andrew 是金融时报(Financial Times)实验室主任,该实验室旨在金融时报开发并推广实践性的 Web 技术.Robe ...

  10. 深入JS第一天:原型和它的小伙伴们(一)

    我在这里不说定义,找点问题,再解决问题. 一.原型 Q1:这样做输出的结果是什么? jQuery= String; jQuery.prototype.say = function () { alert ...