LPOJ 欢迎来到实力至上主义的教室 —— 线段树合并

2023-11-02 10:48

本文主要是介绍LPOJ 欢迎来到实力至上主义的教室 —— 线段树合并,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

题目链接:点我啊╭(╯^╰)╮

题目大意:

    中文题

解题思路:

    线段树维护答案,每个叶子节点维护一颗权值线段树
    合并右 k k k 个点就是先 m e r g e merge merge ,然后二分删掉 k k k 左边的点

核心:线段树合并

正常合并:

#include<bits/stdc++.h>
#define rint register int
#define deb(x) cerr<<#x<<" = "<<(x)<<'\n';
using namespace std;
typedef long long ll;
const int maxn = 1e5 + 5;	
vector <int> v;
int n, q, tot, cnt, a[maxn], root[maxn];
int t0[maxn<<2], ls[maxn*40], rs[maxn*40];
struct node1{int op, x, y, l, r, v, k;
} Q[maxn];struct node{int mx, mi, num;node(){mi = 2e9, mx = num = 0;}node operator + (const node &A)const{node ret;ret.num = num + A.num;ret.mi = min(mi, A.mi);ret.mx = max(mx, A.mx);return ret;}
} t[maxn*40];void update(int &rt, int pos, int l, int r){if(l>pos || r<pos) return;if(!rt) rt = ++tot;if(l == r){t[rt].mi = t[rt].mx = l;t[rt].num++;return;}int m = l + r >> 1;update(ls[rt], pos, l, m);update(rs[rt], pos, m+1, r);t[rt] = t[ls[rt]] + t[rs[rt]];
}void update0(int pos, int c, int l, int r, int rt, int op){if(l>pos || r<pos) return;if(l == r){if(op) update(root[l], c, 1, cnt);if(t[root[l]].mi == 2e9) t0[rt] = 0;else t0[rt] = v[t[root[l]].mx] - v[t[root[l]].mi];return;}int m = l + r >> 1;update0(pos,c,l,m,rt<<1,op);update0(pos,c,m+1,r,rt<<1|1,op);t0[rt] = max(t0[rt<<1], t0[rt<<1|1]);
}int query(int L, int R, int l, int r, int rt){if(l>R || r<L) return 0;if(L<=l && r<=R) return max(0, t0[rt]);int m = l + r >> 1, ret = 0;ret = max(ret, query(L, R, l, m, rt<<1));ret = max(ret, query(L, R, m+1, r, rt<<1|1));return ret;
}int merge(int x, int y){if(!x || !y) return x + y;int now = ++tot;ls[now] = merge(ls[x], ls[y]);rs[now] = merge(rs[x], rs[y]);t[now] = t[ls[now]] + t[rs[now]];return now;
}void gao(int rt, int k, int l, int r){if(l == r){t[rt].num = min(t[rt].num, k);return ;}int cntr = t[rs[rt]].num, m = l + r >> 1;if(cntr < k) gao(ls[rt], k-cntr, l, m);else {ls[rt] = 0;gao(rs[rt], k, m+1, r);}t[rt] = t[ls[rt]] + t[rs[rt]];
}int getid(int x){return lower_bound(v.begin(), v.end(), x) - v.begin();
}int main() {scanf("%d", &n);v.push_back(-1);for(int i=1; i<=n; i++){scanf("%d", a+i);v.push_back(a[i]);}scanf("%d", &q);for(int i=1; i<=q; i++){scanf("%d", &Q[i].op);if(Q[i].op == 1) {scanf("%d%d", &Q[i].x, &Q[i].v);v.push_back(Q[i].v);} else if(Q[i].op == 2) scanf("%d%d", &Q[i].l, &Q[i].r);else scanf("%d%d%d", &Q[i].x, &Q[i].y, &Q[i].k);}sort(v.begin(), v.end());v.erase(unique(v.begin(), v.end()), v.end());cnt = v.size();for(int i=1; i<=n; i++) update0(i, getid(a[i]), 1, n, 1, 1);for(int i=1; i<=q; i++){if(Q[i].op == 1) update0(Q[i].x, getid(Q[i].v), 1, n, 1, 1);else if(Q[i].op == 2) printf("%d\n", query(Q[i].l, Q[i].r, 1, n, 1));else {int x = Q[i].x, y = Q[i].y, k = Q[i].k;root[x] = merge(root[x], root[y]);gao(root[x], k, 1, cnt);root[y] = 0;update0(x,0,1,n,1,0);update0(y,0,1,n,1,0);}}
}

暴力合并:

#include<bits/stdc++.h>
#define rint register int
#define deb(x) cerr<<#x<<" = "<<(x)<<'\n';
using namespace std;
typedef long long ll;
const int maxn = 1e5 + 5;	
vector <int> v;
int n, q, tot, cnt, a[maxn], root[maxn];
int t0[maxn<<2], ls[maxn*150], rs[maxn*150];
struct node1{int op, x, y, l, r, v, k;
} Q[maxn];struct node{int mx, mi, num;node(){mi = 2e9, mx = num = 0;}node operator + (const node &A)const{node ret;ret.num = num + A.num;ret.mi = min(mi, A.mi);ret.mx = max(mx, A.mx);return ret;}
} t[maxn*150];void update(int &rt, int pos, int l, int r){if(l>pos || r<pos) return;if(!rt) rt = ++tot;if(l == r){t[rt].mi = t[rt].mx = l;t[rt].num++;return;}int m = l + r >> 1;update(ls[rt], pos, l, m);update(rs[rt], pos, m+1, r);t[rt] = t[ls[rt]] + t[rs[rt]];
}void update0(int pos, int c, int l, int r, int rt, int op){if(l>pos || r<pos) return;if(l == r){if(op) update(root[l], c, 1, cnt);if(t[root[l]].mi == 2e9) t0[rt] = 0;else t0[rt] = v[t[root[l]].mx] - v[t[root[l]].mi];return;}int m = l + r >> 1;update0(pos,c,l,m,rt<<1,op);update0(pos,c,m+1,r,rt<<1|1,op);t0[rt] = max(t0[rt<<1], t0[rt<<1|1]);
}int query(int L, int R, int l, int r, int rt){if(l>R || r<L) return 0;if(L<=l && r<=R) return max(0, t0[rt]);int m = l + r >> 1, ret = 0;ret = max(ret, query(L, R, l, m, rt<<1));ret = max(ret, query(L, R, m+1, r, rt<<1|1));return ret;
}int merge(int x, int y, int k, int l, int r){if((!x && !y) || !k) return 0;if(l == r){int now = ++tot, cnt = t[x].num + t[y].num;t[now].num = min(k, cnt);if(t[now].num) t[now].mi = t[now].mx = l;return now;}int m = l + r >> 1;if(!x){int now = ++tot, cnt = t[rs[y]].num;if(cnt >= k){ls[now] = 0;rs[now] = merge(0, rs[y], k, m+1, r);} else {ls[now] = merge(0, ls[y], k-cnt, l, m);rs[now] = rs[y];}t[now] = t[ls[now]] + t[rs[now]];return now;}	if(!y){int now = ++tot, cnt = t[rs[x]].num;if(cnt >= k){ls[now] = 0;rs[now] = merge(rs[x], 0, k, m+1, r);} else {ls[now] = merge(ls[x], 0, k-cnt, l, m);rs[now] = rs[x];}t[now] = t[ls[now]] + t[rs[now]];return now;}int now = ++tot, cnt = t[rs[x]].num + t[rs[y]].num;if(cnt >= k){ls[now] = 0;rs[now] = merge(rs[x], rs[y], k, m+1, r);} else {ls[now] = merge(ls[x], ls[y], k-cnt, l, m);rs[now] = merge(rs[x], rs[y], cnt, m+1, r);}t[now] = t[ls[now]] + t[rs[now]];return now;
}int getid(int x){return lower_bound(v.begin(), v.end(), x) - v.begin();
}int main() {scanf("%d", &n);v.push_back(-1);for(int i=1; i<=n; i++){scanf("%d", a+i);v.push_back(a[i]);}scanf("%d", &q);for(int i=1; i<=q; i++){scanf("%d", &Q[i].op);if(Q[i].op == 1) {scanf("%d%d", &Q[i].x, &Q[i].v);v.push_back(Q[i].v);} else if(Q[i].op == 2) scanf("%d%d", &Q[i].l, &Q[i].r);else scanf("%d%d%d", &Q[i].x, &Q[i].y, &Q[i].k);}sort(v.begin(), v.end());v.erase(unique(v.begin(), v.end()), v.end());cnt = v.size();for(int i=1; i<=n; i++) update0(i, getid(a[i]), 1, n, 1, 1);for(int i=1; i<=q; i++){if(Q[i].op == 1) update0(Q[i].x, getid(Q[i].v), 1, n, 1, 1);else if(Q[i].op == 2) printf("%d\n", query(Q[i].l, Q[i].r, 1, n, 1));else {int x = Q[i].x, y = Q[i].y, k = Q[i].k;root[x] = merge(root[x], root[y], k, 1, cnt);root[y] = 0;update0(x,0,1,n,1,0);update0(y,0,1,n,1,0);}}
}

这篇关于LPOJ 欢迎来到实力至上主义的教室 —— 线段树合并的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/330077

相关文章

hdu2241(二分+合并数组)

题意:判断是否存在a+b+c = x,a,b,c分别属于集合A,B,C 如果用暴力会超时,所以这里用到了数组合并,将b,c数组合并成d,d数组存的是b,c数组元素的和,然后对d数组进行二分就可以了 代码如下(附注释): #include<iostream>#include<algorithm>#include<cstring>#include<stack>#include<que

poj3468(线段树成段更新模板题)

题意:包括两个操作:1、将[a.b]上的数字加上v;2、查询区间[a,b]上的和 下面的介绍是下解题思路: 首先介绍  lazy-tag思想:用一个变量记录每一个线段树节点的变化值,当这部分线段的一致性被破坏我们就将这个变化值传递给子区间,大大增加了线段树的效率。 比如现在需要对[a,b]区间值进行加c操作,那么就从根节点[1,n]开始调用update函数进行操作,如果刚好执行到一个子节点,

hdu1394(线段树点更新的应用)

题意:求一个序列经过一定的操作得到的序列的最小逆序数 这题会用到逆序数的一个性质,在0到n-1这些数字组成的乱序排列,将第一个数字A移到最后一位,得到的逆序数为res-a+(n-a-1) 知道上面的知识点后,可以用暴力来解 代码如下: #include<iostream>#include<algorithm>#include<cstring>#include<stack>#in

hdu1689(线段树成段更新)

两种操作:1、set区间[a,b]上数字为v;2、查询[ 1 , n ]上的sum 代码如下: #include<iostream>#include<algorithm>#include<cstring>#include<stack>#include<queue>#include<set>#include<map>#include<stdio.h>#include<stdl

hdu 1754 I Hate It(线段树,单点更新,区间最值)

题意是求一个线段中的最大数。 线段树的模板题,试用了一下交大的模板。效率有点略低。 代码: #include <stdio.h>#include <string.h>#define TREE_SIZE (1 << (20))//const int TREE_SIZE = 200000 + 10;int max(int a, int b){return a > b ? a :

hdu 1166 敌兵布阵(树状数组 or 线段树)

题意是求一个线段的和,在线段上可以进行加减的修改。 树状数组的模板题。 代码: #include <stdio.h>#include <string.h>const int maxn = 50000 + 1;int c[maxn];int n;int lowbit(int x){return x & -x;}void add(int x, int num){while

day-51 合并零之间的节点

思路 直接遍历链表即可,遇到val=0跳过,val非零则加在一起,最后返回即可 解题过程 返回链表可以有头结点,方便插入,返回head.next Code /*** Definition for singly-linked list.* public class ListNode {* int val;* ListNode next;* ListNode() {}*

poj 1127 线段相交的判定

题意: 有n根木棍,每根的端点坐标分别是 px, py, qx, qy。 判断每对木棍是否相连,当他们之间有公共点时,就认为他们相连。 并且通过相连的木棍相连的木棍也是相连的。 解析: 线段相交的判定。 首先,模板中的线段相交是不判端点的,所以要加一个端点在直线上的判定; 然后,端点在直线上的判定这个函数是不判定两个端点是同一个端点的情况的,所以要加是否端点相等的判断。 最后

HDU4737线段树

题目大意:给定一系列数,F(i,j)表示对从ai到aj连续求或运算,(i<=j)求F(i,j)<=m的总数。 const int Max_N = 100008 ;int sum[Max_N<<2] , x[Max_N] ;int n , m ;void push_up(int t){sum[t] = sum[t<<1] | sum[t<<1|1] ;}void upd

zoj 1721 判断2条线段(完全)相交

给出起点,终点,与一些障碍线段。 求起点到终点的最短路。 枚举2点的距离,然后最短路。 2点可达条件:没有线段与这2点所构成的线段(完全)相交。 const double eps = 1e-8 ;double add(double x , double y){if(fabs(x+y) < eps*(fabs(x) + fabs(y))) return 0 ;return x + y ;