【bzoj2157】【旅游】

2157: 旅游

Time Limit: 10 Sec Memory Limit: 259 MB

Submit: 419 Solved: 248

[Submit][Status][Discuss]

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

3

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

3

2

1

-1

5

3

思路：跟Tree一样，只是多了区间求和和区间最小值。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=100100;
int n,m,point[N],next[N*10],tot=1,siz[N]={0},fa[N][20]={0},son[N]={0};
int deep[N]={0},num=0,pos[N]={0},belong[N]={0},r[N]={0},map[N]={0},de[N*4]={0};
struct S{
    int st,en,va;
}aa[N*10];
struct C{
    int maxn,minn,sum;
}tr[N*4];
bool use[N];
inline void add(int x,int y,int z)
{
    tot+=1;next[tot]=point[x];point[x]=tot;
    aa[tot].st=x;aa[tot].en=y;aa[tot].va=z;
    tot+=1;next[tot]=point[y];point[y]=tot;
    aa[tot].st=y;aa[tot].en=x;aa[tot].va=z;
}
inline void dfs_1(int x)
{
    int i;
    siz[x]=1;
    use[x]=false;
    for(i=1;i<=14;++i){
        if(deep[x]<(1<<i)) break;
        fa[x][i]=fa[fa[x][i-1]][i-1];
    }
    for(i=point[x];i;i=next[i])
      if(use[aa[i].en]){
        fa[aa[i].en][0]=x;
        deep[aa[i].en]=deep[x]+1;
        dfs_1(aa[i].en);
        siz[x]+=siz[aa[i].en];
      }
}
inline void dfs_2(int x,int y)
{
    int i,k=0;
    num+=1;
    belong[x]=y;
    for(i=point[x];i;i=next[i]){
        if(deep[aa[i].st]>deep[aa[i].en]) map[aa[i].st]=i>>1;
        else map[aa[i].en]=i>>1;
    }
    for(i=point[x];i;i=next[i])
      if(deep[x]<deep[aa[i].en]&&siz[k]<siz[aa[i].en])
        k=aa[i].en,son[x]=i>>1,pos[i>>1]=num;
    if(k==0) {num-=1; return ;}
    dfs_2(k,y);
    for(i=point[x];i;i=next[i])
      if(deep[x]<deep[aa[i].en]&&k!=aa[i].en)
        num+=1,pos[i>>1]=num,dfs_2(aa[i].en,aa[i].en);
}
inline int LCA(int x,int y)
{
    int i;
    if(deep[x]<deep[y])swap(x,y);
    int t=deep[x]-deep[y];
    for(i=0;i<=14;++i)
        if(t&(1<<i))x=fa[x][i];
    for(i=14;i>=0;--i)
        if(fa[x][i]!=fa[y][i])
        {x=fa[x][i];y=fa[y][i];}
    if(x==y)return x;
    else return fa[x][0];
}
#define mid (l+r)/2
#define L k<<1,l,mid
#define R k<<1|1,mid+1,r
inline void solve(int k)
{
    swap(tr[k].maxn,tr[k].minn);
    tr[k].maxn=-tr[k].maxn;
    tr[k].minn=-tr[k].minn;
    tr[k].sum=-tr[k].sum;
}
inline void pushdown(int k)
{
    if(de[k]!=0){
        if((k<<1)<N*4){
            solve(k<<1); de[k<<1]=!de[k<<1];
            solve(k<<1|1); de[k<<1|1]=!de[k<<1|1];
        }
        de[k]=0;
    }
}
inline void insert(int k,int l,int r,int x,int y)
{
    pushdown(k);
    if(l==r&&l==x){
        tr[k].maxn=tr[k].minn=tr[k].sum=y;
        return ;
    }
    if(x<=mid) insert(L,x,y);
    else insert(R,x,y);
    tr[k].maxn=max(tr[k<<1].maxn,tr[k<<1|1].maxn);
    tr[k].minn=min(tr[k<<1].minn,tr[k<<1|1].minn);
    tr[k].sum=tr[k<<1].sum+tr[k<<1|1].sum;
}
inline int qurry(int k,int l,int r,int x,int y,int kind)
{
    int maxn=-210000000,minn=210000000,sum=0;
    pushdown(k);
    if(x<=l&&y>=r){
        if(kind==0) return tr[k].maxn;
        if(kind==1) return tr[k].minn;
        if(kind==2) return tr[k].sum;
    }
    if(x<=mid){
        if(kind==0) maxn=max(maxn,qurry(L,x,y,kind));
        if(kind==1) minn=min(minn,qurry(L,x,y,kind));
        if(kind==2) sum+=qurry(L,x,y,kind);
    }
    if(y>mid){
        if(kind==0) maxn=max(maxn,qurry(R,x,y,kind));
        if(kind==1) minn=min(minn,qurry(R,x,y,kind));
        if(kind==2) sum+=qurry(R,x,y,kind);
    }
    if(kind==0) return maxn;
    if(kind==1) return minn;
    if(kind==2) return sum;
}
inline int work(int x,int y,int kind)
{
    int maxn=-210000000,z,minn=210000000,sum=0;
    if(map[x]==0){if(kind==0) return maxn;if(kind==1) return minn;if(kind==2) return sum;}
    if(belong[x]==x&&x!=y){
        if(kind==0) maxn=max(maxn,qurry(1,1,n-1,pos[map[x]],pos[map[x]],kind));
        if(kind==1) minn=min(minn,qurry(1,1,n-1,pos[map[x]],pos[map[x]],kind));
        if(kind==2) sum+=qurry(1,1,n-1,pos[map[x]],pos[map[x]],kind);
        x=fa[x][0];
    }
    if(map[x]==0){if(kind==0) return maxn;if(kind==1) return minn;if(kind==2) return sum;}
    while(belong[x]!=belong[y]){
        if(kind==0) maxn=max(maxn,qurry(1,1,n-1,pos[son[belong[x]]],pos[map[x]],kind));
        if(kind==1) minn=min(minn,qurry(1,1,n-1,pos[son[belong[x]]],pos[map[x]],kind));
        if(kind==2) sum+=qurry(1,1,n-1,pos[son[belong[x]]],pos[map[x]],kind);
        z=fa[belong[x]][0];
        if(kind==0) maxn=max(maxn,qurry(1,1,n-1,pos[map[belong[x]]],pos[map[belong[x]]],kind));
        if(kind==1) minn=min(minn,qurry(1,1,n-1,pos[map[belong[x]]],pos[map[belong[x]]],kind));
        if(kind==2) sum+=qurry(1,1,n-1,pos[map[belong[x]]],pos[map[belong[x]]],kind);
        x=z;
    }
    if(map[x]==0){if(kind==0) return maxn;if(kind==1) return minn;if(kind==2) return sum;}
    if(kind==0) maxn=max(maxn,qurry(1,1,n-1,pos[son[y]],pos[map[x]],kind));
    if(kind==1) minn=min(minn,qurry(1,1,n-1,pos[son[y]],pos[map[x]],kind));
    if(kind==2) sum+=qurry(1,1,n-1,pos[son[y]],pos[map[x]],kind);
    if(kind==0) return maxn;if(kind==1) return minn;if(kind==2) return sum;
}
inline void change(int k,int l,int r,int x,int y)
{
    pushdown(k);
    if(x<=l&&y>=r){
        solve(k); de[k]=1;
        return ;
    }
    if(x<=mid) change(L,x,y);
    if(y>mid) change(R,x,y);
    tr[k].maxn=max(tr[k<<1].maxn,tr[k<<1|1].maxn);
    tr[k].minn=min(tr[k<<1].minn,tr[k<<1|1].minn);
    tr[k].sum=tr[k<<1].sum+tr[k<<1|1].sum;
}
inline void negate_(int x,int y)
{
    int z;
    if(map[x]==0) return ;
    if(belong[x]==x&&x!=y){
        change(1,1,n-1,pos[map[x]],pos[map[x]]);
        x=fa[x][0];
    }
    if(map[x]==0) return ;
    while(belong[x]!=belong[y]){
        change(1,1,n-1,pos[son[belong[x]]],pos[map[x]]);
        z=fa[belong[x]][0];
        change(1,1,n-1,pos[map[belong[x]]],pos[map[belong[x]]]);
        x=z;
    }
    if(map[x]==0) return ;
    change(1,1,n-1,pos[son[y]],pos[map[x]]);
    return ;
}
int main()
{
    int i,j,x,y,z;
    char ch[10];
    memset(use,1,sizeof(use));
    scanf("%d",&n);
    for(i=1;i<n;++i){
        scanf("%d%d%d",&x,&y,&z);
        x+=1;y+=1;
        add(x,y,z);r[i]=z;
    }
    dfs_1(1);
    dfs_2(1,1);
    for(i=1;i<n;++i) insert(1,1,n-1,pos[i],r[i]);
    scanf("%d",&m);
    while(m--){
        scanf("%*c%s",&ch);
        scanf("%d%d",&x,&y);
        if(ch[0]==‘C‘) insert(1,1,n-1,pos[x],y);
        else{
            x+=1;y+=1;
            int lca=LCA(x,y);
            if(ch[0]==‘N‘) negate_(x,lca),negate_(y,lca);
            if(ch[0]==‘S‘) printf("%d\n",work(x,lca,2)+work(y,lca,2));
            if(ch[0]==‘M‘){
                if(ch[1]==‘A‘) printf("%d\n",max(work(x,lca,0),work(y,lca,0)));
                if(ch[1]==‘I‘) printf("%d\n",min(work(x,lca,1),work(y,lca,1)));
            }
        }
    }
}
时间： 2024-08-05 06:25:47

【bzoj2157】【旅游】的相关文章

bzoj2157旅游

题目: 2157: 旅游 Time Limit: 10 Sec  Memory Limit: 259 MBSubmit: 569  Solved: 310[Submit][Status][Discuss] Description Ray 乐忠于旅游,这次他来到了T 城.T 城是一个水上城市,一共有 N 个景点,有些景点之间会用一座桥连接.为了方便游客到达每个景点但又为了节约成本,T 城的任意两个景点之间有且只有一条路径.换句话说, T 城中只有N − 1 座桥.Ray 发现,有些桥上可以看到美丽

[bzoj2157]旅游 (lct)

这个应该也算裸的模板题吧..主要是边权的问题,对于每条边u->v,我们可以新建一个节点代替他,把边的信息弄到新的点上,就变成u->x->v了... 当然了这样的话要防止u和v这些没用的点影响到实际的结果....这个可以初始化的时候解决 话说如果写链剖的话也可以用这样的姿势,就不用像以前那样特判lca了.. 上传的时候少更新了一个值结果调了半天= = 1 #include<cstdio> 2 #include<iostream> 3 #include<math

bzoj2157: 旅游

又是lct裸题(据说暴力能过?) 写着略蛋疼 1 #include<iostream> 2 #include<cstring> 3 #include<cstdlib> 4 #include<cstdio> 5 #include<algorithm> 6 7 using namespace std; 8 9 const int Maxn=200010; 10 11 int ch[Maxn][2],w[Maxn],sum[Maxn],mn[Maxn]

【树链剖分】【线段树】bzoj2157 旅游

#include<cstdio> #include<algorithm> using namespace std; #define INF 2147483647 #define N 20001 #define lson rt<<1,l,m #define rson rt<<1|1,m+1,r int v[N<<1],first[N],next[N<<1],en,bw[N<<1],dw[N<<1]; void A

BZOJ2157: 旅游 树链剖分 线段树

http://www.lydsy.com/JudgeOnline/problem.php?id=2157 在对树中数据进行改动的时候需要很多pushdown(具体操作见代码),不然会wa,大概原因和线段树区间修改需要很多pushup是一样的. 这个轻重链的方法特别好用,虽然第一次写树链剖分但是容易理解又有优秀复杂度的结构让人情不自禁orz. (后来发现很久以前学lca的时候就学了树链剖分只不过忘了,mdzz) 因为忘了去掉测试代码的freopen,re了4次(虽然有三次就算不re也wa),发现一

【BZOJ2157】旅游 LCT

模板T,SB的DMoon..其实样例也是中国好样例...一开始不会复制,yangyang:找到“sample input”按住shift,按page down.... 1 #include <iostream> 2 #include <cstdio> 3 #define inf 0x7fffffff 4 #define N 20010 5 #define M 20010 6 using namespace std; 7 int n,m; 8 struct node *null; 9

【bzoj2157】旅游

题目描述 Ray 乐忠于旅游,这次他来到了T 城.T 城是一个水上城市,一共有 N 个景点,有些景点之间会用一座桥连接.为了方便游客到达每个景点但又为了节约成本,T 城的任意两个景点之间有且只有一条路径.换句话说, T 城中只有N ? 1 座桥. Ray 发现,有些桥上可以看到美丽的景色,让人心情愉悦,但有些桥狭窄泥泞,令人烦躁.于是,他给每座桥定义一个愉悦度w,也就是说,Ray 经过这座桥会增加w 的愉悦度,这或许是正的也可能是负的.有时,Ray 看待同一座桥的心情也会发生改变. 现在,Ray

【BZOJ2157】旅游 裸树链剖分

#include <stdio.h> int main() { puts("转载请注明出处[vmurder]谢谢"); puts("网址:blog.csdn.net/vmurder/article/details/43965845"); } 重写大发好!!!!! ****什么题解都没有,水题一道, 挂了就去调,调不过就去重写. 代码: #include <cstdio> #include <cstring> #include &l

广州去厦门旅游攻略--（转自 #散文吧网站#）

广州去厦门旅游攻略 发布时间:2016-12-11 17:30 厦门由厦门岛.离岛鼓浪屿.内陆九龙江南岸海沧半岛.集美半岛.翔安区以及同安等组成,陆地面积1699.39Km2,海域面积300多平方公里.厦门的主体--厦门岛南北长13.7公里,东西宽12.5公里,面积约为128.14Km2.是厦门的主要岛屿,也是厦门第一大岛屿.厦门岛是厦门经济特区的发祥地,岛上有厦门的商业和政治中心.各国殖民者最初居住的地方鼓浪屿就在厦门岛西南部.今天小编带给大家的是广州去厦门旅游攻略, 希望对大家有帮助. 厦门