【BZOJ2157】旅游 树链剖分+线段树
【BZOJ2157】旅游
Description
Ray 乐忠于旅游,这次他来到了T 城。T 城是一个水上城市,一共有 N 个景点,有些景点之间会用一座桥连接。为了方便游客到达每个景点但又为了节约成本,T 城的任意两个景点之间有且只有一条路径。换句话说, T 城中只有N − 1 座桥。Ray 发现,有些桥上可以看到美丽的景色,让人心情愉悦,但有些桥狭窄泥泞,令人烦躁。于是,他给每座桥定义一个愉悦度w,也就是说,Ray 经过这座桥会增加w 的愉悦度,这或许是正的也可能是负的。有时,Ray 看待同一座桥的心情也会发生改变。现在,Ray 想让你帮他计算从u 景点到v 景点能获得的总愉悦度。有时,他还想知道某段路上最美丽的桥所提供的最大愉悦度,或是某段路上最糟糕的一座桥提供的最低愉悦度。
Input
输入的第一行包含一个整数N,表示T 城中的景点个数。景点编号为 0...N − 1。接下来N − 1 行,每行三个整数u、v 和w,表示有一条u 到v,使 Ray 愉悦度增加w 的桥。桥的编号为1...N − 1。|w| <= 1000。输入的第N + 1 行包含一个整数M,表示Ray 的操作数目。接下来有M 行,每行描述了一个操作,操作有如下五种形式: C i w,表示Ray 对于经过第i 座桥的愉悦度变成了w。 N u v,表示Ray 对于经过景点u 到v 的路径上的每一座桥的愉悦度都变成原来的相反数。 SUM u v,表示询问从景点u 到v 所获得的总愉悦度。 MAX u v,表示询问从景点u 到v 的路径上的所有桥中某一座桥所提供的最大愉悦度。 MIN u v,表示询问从景点u 到v 的路径上的所有桥中某一座桥所提供的最小愉悦度。测试数据保证,任意时刻,Ray 对于经过每一座桥的愉悦度的绝对值小于等于1000。
Output
对于每一个询问(操作S、MAX 和MIN),输出答案。
Sample Input
0 1 1
1 2 2
8
SUM 0 2
MAX 0 2
N 0 1
SUM 0 2
MIN 0 2
C 1 3
SUM 0 2
MAX 0 2
Sample Output
2
1
-1
5
3
HINT
一共有10 个数据,对于第i (1 <= i <= 10) 个数据, N = M = i * 2000。
题解:又是树剖模板题,已经不知道该注意什么了,仍然会因为没写pushdown而狂WA不止~
搬运个数据生成器吧~ //from GXZlegend
- #include <cstdio>
- #include <cstring>
- #include <cstdlib>
- #include <ctime>
- #include <iostream>
- #include <string>
- #include <algorithm>
- using namespace std;
- int n , m;
- void build_tree() {
- for( int i = 1 ; i < n ; ++i )
- printf( "%d %d %d\n" , i , rand() % i , rand() % 11 * ( rand() & 1 ? -1 : 1 ) );
- }
- string s[ 5 ] = { "C" , "N" , "SUM" , "MAX" , "MIN" };
- void build_query() {
- while( m-- ) {
- int op = rand() % 5;
- cout << s[ op ] << ' ';
- if( ! op ) {
- printf( "%d %d" , rand() % ( n - 1 ) + 1 , rand() % 11 * ( rand() & 1 ? -1 : 1 ) );
- } else {
- int u = rand() % n , v = rand() % n;
- while( v == u ) v =rand() % n;
- printf( "%d %d" , u , v );
- }
- putchar( '\n' );
- }
- }
- int main() {
- srand( time( NULL ) );
- n =1000 , m = 2000;
- cout << n << ' ' << "\n";
- build_tree();
- cout << m << "\n";
- build_query();
- return 0;
- }
- #include <cstdio>
- #include <cstring>
- #include <iostream>
- #include <queue>
- #define lson x<<1
- #define rson x<<1|1
- using namespace std;
- const int maxn=200010;
- int head[maxn],next[maxn<<1],to[maxn<<1],val[maxn<<1],v[maxn<<1],p[maxn];
- int dep[maxn],siz[maxn],son[maxn],top[maxn],fa[maxn],bel[maxn];
- int s[maxn],rs[maxn],sm[maxn],sn[maxn];
- int n,m,cnt,tot;
- char str[5];
- int rd()
- {
- int ret=0,flag=1; char gc=getchar();
- while(gc<'0'||gc>'9'){if(gc=='-') flag=-flag; gc=getchar();}
- while(gc>='0'&&gc<='9') ret=ret*10+gc-'0',gc=getchar();
- return ret*flag;
- }
- void add(int a,int b,int c)
- {
- to[cnt]=b,val[cnt]=c,next[cnt]=head[a],head[a]=cnt++;
- }
- void pushup(int x)
- {
- s[x]=s[lson]+s[rson],sm[x]=max(sm[lson],sm[rson]),sn[x]=min(sn[lson],sn[rson]);
- }
- void pushdown(int x)
- {
- if(rs[x])
- {
- sm[lson]=-sm[lson],sm[rson]=-sm[rson],sn[lson]=-sn[lson],sn[rson]=-sn[rson];
- swap(sm[lson],sn[lson]),swap(sm[rson],sn[rson]);
- s[lson]=-s[lson],s[rson]=-s[rson],rs[lson]^=1,rs[rson]^=1,rs[x]=0;
- }
- }
- void dfs1(int x)
- {
- siz[x]=1;
- for(int i=head[x];i!=-1;i=next[i])
- {
- if(to[i]!=fa[x])
- {
- fa[to[i]]=x,dep[to[i]]=dep[x]+1,v[to[i]]=val[i],bel[i>>1]=to[i];
- dfs1(to[i]);
- siz[x]+=siz[to[i]];
- if(siz[to[i]]>siz[son[x]]) son[x]=to[i];
- }
- }
- }
- void updata(int l,int r,int x,int a,int b)
- {
- if(l==r)
- {
- s[x]=sm[x]=sn[x]=b;
- return ;
- }
- pushdown(x);
- int mid=l+r>>1;
- if(a<=mid) updata(l,mid,lson,a,b);
- else updata(mid+1,r,rson,a,b);
- pushup(x);
- }
- void dfs2(int x,int tp)
- {
- top[x]=tp,p[x]=++tot;
- updata(1,n,1,p[x],v[x]);
- if(son[x]) dfs2(son[x],tp);
- for(int i=head[x];i!=-1;i=next[i]) if(to[i]!=fa[x]&&to[i]!=son[x]) dfs2(to[i],to[i]);
- }
- void uprs(int l,int r,int x,int a,int b)
- {
- if(a<=l&&r<=b)
- {
- sm[x]=-sm[x],sn[x]=-sn[x],swap(sm[x],sn[x]),s[x]=-s[x],rs[x]^=1;
- return ;
- }
- pushdown(x);
- int mid=l+r>>1;
- if(a<=mid) uprs(l,mid,lson,a,b);
- if(b>mid) uprs(mid+1,r,rson,a,b);
- pushup(x);
- }
- int qs(int l,int r,int x,int a,int b)
- {
- if(a<=l&&r<=b) return s[x];
- pushdown(x);
- int mid=l+r>>1;
- if(b<=mid) return qs(l,mid,lson,a,b);
- if(a>mid) return qs(mid+1,r,rson,a,b);
- return qs(l,mid,lson,a,b)+qs(mid+1,r,rson,a,b);
- }
- int qm(int l,int r,int x,int a,int b)
- {
- if(a<=l&&r<=b) return sm[x];
- pushdown(x);
- int mid=l+r>>1;
- if(b<=mid) return qm(l,mid,lson,a,b);
- if(a>mid) return qm(mid+1,r,rson,a,b);
- return max(qm(l,mid,lson,a,b),qm(mid+1,r,rson,a,b));
- }
- int qn(int l,int r,int x,int a,int b)
- {
- if(a<=l&&r<=b) return sn[x];
- pushdown(x);
- int mid=l+r>>1;
- if(b<=mid) return qn(l,mid,lson,a,b);
- if(a>mid) return qn(mid+1,r,rson,a,b);
- return min(qn(l,mid,lson,a,b),qn(mid+1,r,rson,a,b));
- }
- void N(int x,int y)
- {
- while(top[x]!=top[y])
- {
- if(dep[top[x]]<dep[top[y]]) swap(x,y);
- uprs(1,n,1,p[top[x]],p[x]),x=fa[top[x]];
- }
- if(x==y) return ;
- if(dep[x]>dep[y]) swap(x,y);
- uprs(1,n,1,p[x]+1,p[y]);
- }
- int S(int x,int y)
- {
- int ret=0;
- while(top[x]!=top[y])
- {
- if(dep[top[x]]<dep[top[y]]) swap(x,y);
- ret+=qs(1,n,1,p[top[x]],p[x]),x=fa[top[x]];
- }
- if(x==y) return ret;
- if(dep[x]>dep[y]) swap(x,y);
- return ret+qs(1,n,1,p[x]+1,p[y]);
- }
- int MAX(int x,int y)
- {
- int ret=-(1<<30);
- while(top[x]!=top[y])
- {
- if(dep[top[x]]<dep[top[y]]) swap(x,y);
- ret=max(ret,qm(1,n,1,p[top[x]],p[x])),x=fa[top[x]];
- }
- if(x==y) return ret;
- if(dep[x]>dep[y]) swap(x,y);
- return max(ret,qm(1,n,1,p[x]+1,p[y]));
- }
- int MIN(int x,int y)
- {
- int ret=1<<30;
- while(top[x]!=top[y])
- {
- if(dep[top[x]]<dep[top[y]]) swap(x,y);
- ret=min(ret,qn(1,n,1,p[top[x]],p[x])),x=fa[top[x]];
- }
- if(x==y) return ret;
- if(dep[x]>dep[y]) swap(x,y);
- return min(ret,qn(1,n,1,p[x]+1,p[y]));
- }
- int main()
- {
- memset(head,-1,sizeof(head));
- memset(sm,0x80,sizeof(sm));
- memset(sn,0x3f,sizeof(sn));
- n=rd();
- int i,a,b,c;
- for(i=1;i<n;i++)
- {
- a=rd()+1,b=rd()+1,c=rd();
- add(a,b,c),add(b,a,c);
- }
- dep[1]=1;
- dfs1(1),dfs2(1,1);
- m=rd();
- for(i=1;i<=m;i++)
- {
- scanf("%s",str);
- a=rd(),b=rd();
- if(str[0]=='C') updata(1,n,1,p[bel[a-1]],b);
- if(str[0]=='N') N(a+1,b+1);
- if(str[0]=='S') printf("%d\n",S(a+1,b+1));
- if(str[0]=='M'&&str[1]=='A') printf("%d\n",MAX(a+1,b+1));
- if(str[0]=='M'&&str[1]=='I') printf("%d\n",MIN(a+1,b+1));
- }
- return 0;
- }
【BZOJ2157】旅游 树链剖分+线段树的更多相关文章
- BZOJ2157旅游——树链剖分+线段树
题目描述 Ray 乐忠于旅游,这次他来到了T 城.T 城是一个水上城市,一共有 N 个景点,有些景点之间会用一座桥连接.为了方便游客到达每个景点但又为了节约成本,T 城的任意两个景点之间有且只有一条路 ...
- bzoj 2157: 旅游【树链剖分+线段树】
裸的树链剖分+线段树 但是要注意一个地方--我WA了好几次才发现取完相反数之后max值和min值是要交换的-- #include<iostream> #include<cstdio& ...
- 【BZOJ-2325】道馆之战 树链剖分 + 线段树
2325: [ZJOI2011]道馆之战 Time Limit: 40 Sec Memory Limit: 256 MBSubmit: 1153 Solved: 421[Submit][Statu ...
- 【BZOJ2243】[SDOI2011]染色 树链剖分+线段树
[BZOJ2243][SDOI2011]染色 Description 给定一棵有n个节点的无根树和m个操作,操作有2类: 1.将节点a到节点b路径上所有点都染成颜色c: 2.询问节点a到节点b路径上的 ...
- BZOJ2243 (树链剖分+线段树)
Problem 染色(BZOJ2243) 题目大意 给定一颗树,每个节点上有一种颜色. 要求支持两种操作: 操作1:将a->b上所有点染成一种颜色. 操作2:询问a->b上的颜色段数量. ...
- POJ3237 (树链剖分+线段树)
Problem Tree (POJ3237) 题目大意 给定一颗树,有边权. 要求支持三种操作: 操作一:更改某条边的权值. 操作二:将某条路径上的边权取反. 操作三:询问某条路径上的最大权值. 解题 ...
- bzoj4034 (树链剖分+线段树)
Problem T2 (bzoj4034 HAOI2015) 题目大意 给定一颗树,1为根节点,要求支持三种操作. 操作 1 :把某个节点 x 的点权增加 a . 操作 2 :把某个节点 x 为根的子 ...
- HDU4897 (树链剖分+线段树)
Problem Little Devil I (HDU4897) 题目大意 给定一棵树,每条边的颜色为黑或白,起始时均为白. 支持3种操作: 操作1:将a->b的路径中的所有边的颜色翻转. 操作 ...
- Aizu 2450 Do use segment tree 树链剖分+线段树
Do use segment tree Time Limit: 1 Sec Memory Limit: 256 MB 题目连接 http://www.bnuoj.com/v3/problem_show ...
- 【POJ3237】Tree(树链剖分+线段树)
Description You are given a tree with N nodes. The tree’s nodes are numbered 1 through N and its edg ...
随机推荐
- CDH安装kafka
摘要:前言其实cloudera已经做了这个事了,只是把kafka的包和cdh的parcel包分离了,只要我们把分离开的kafka的服务描述jar包和服务parcel包下载了,就可以实现完美集成了.具体 ...
- DOS建立带点的文件夹
md 表示make directory 创建文件夹 rd 表示remove directory 删除文件夹 我们用一个实例来说明一切吧. 1.按 键盘Windows徽标键+R 来打开 运行 对话框 2 ...
- PHP程序员的10个有用的技巧和教程
PHP被定义为一个通用的服务器端脚本语言,它基本上是专为Web开发的重要目的.借助PHP可以创建动态和惊人的网页效果,他是被嵌入在一个HTML源文件的服务器端的脚本语言之一.因此,它成为最流行的,也是 ...
- js学习笔记19----getElementsByClassName函数封装
js里面的getElementsByClassName()方法可通过某个class名获取到元素,在标准浏览器下可使用,在非标准浏览器下不可用.为了能够让这个方法兼容所有的浏览器,可以封装成如下函数: ...
- SQL计算表的列数
select count(syscolumns.name) from syscolumns , sysobjects where syscolumns.id = sysobjects.id and s ...
- Python3.4下使用sqlalchemy
一. 1.用sudo apt-get install python3-numpy之后,会默认把numpy安装到 /usr/lib/python3/dist-packages目录下,而且版本比较低. ...
- html中可以自定义属性,,,妈的竟然才知道..
html中可以自定义属性,,,妈的竟然才知道.. <input userinfo="没见过帅哥呀" />
- pclzip 压缩文件与解压
类PclZip.class.php下载:PclZip.rar<?php header("Content-type: text/html; charset=utf-8"); f ...
- 字符串池化 python
前言 在 Python 中经常通过内存池化技术来提高其性能,那么问题来了,在什么情况下会池化呢? 让我们通过几个例子进行一下理解一下. 预备知识 在查看例子之前,首先要提 python 中的一个函数 ...
- How many virtual users do I need? 计算需要的vuser
基本公式:--------English:Total Transations = TPS * Vuser==>Vuser = Total Transations / TPS=========== ...