Snacks HDU 5692 dfs序列+线段树

时间:2021-07-10 14:56:13

Snacks HDU 5692 dfs序列+线段树

题意

百度科技园内有n个零食机,零食机之间通过n−1条路相互连通。每个零食机都有一个值v,表示为小度熊提供零食的价值。

由于零食被频繁的消耗和补充,零食机的价值v会时常发生变化。小度熊只能从编号为0的零食机出发,并且每个零食机至多经过一次。另外,小度熊会对某个零食机的零食有所偏爱,要求路线上必须有那个零食机。

为小度熊规划一个路线,使得路线上的价值总和最大

输入输出:

输入数据第一行是一个整数T(T≤10),表示有T组测试数据。

对于每组数据,包含两个整数n,m(1≤n,m≤100000),表示有n个零食机,m次操作。

接下来n−1行,每行两个整数x和y(0≤x,y<n),表示编号为x的零食机与编号为y的零食机相连。

接下来一行由n个数组成,表示从编号为0到编号为n−1的零食机的初始价值v(|v|<100000)。

接下来m行,有两种操作:0 x y,表示编号为x的零食机的价值变为y;1 x,表示询问从编号为0的零食机出发,必须经过编号为x零食机的路线中,价值总和的最大值。

本题可能栈溢出,辛苦同学们提交语言选择c++,并在代码的第一行加上:

#pragma comment(linker, "/STACK:1024000000,1024000000")

题解思路

dfs序列的题,这是第一道,当然不出意外,我不会做。但是我找到一篇写的非常好的博客,详细解释了这个题。点我链接

思路就是用dfs序列,将树“拍”成一维的一串,这样就可以用线段树或者树状数组进行接下来的操作了。

至于为什么要转化为一维的,因为简单啊,并且用来处理一维数据的方法也很多。

代码实现

两种形似,第一个是数组形式的线段树,有注释。第二种是使用结构体线段树,和第一种基本一样。

#pragma comment(linker, "/STACK:1024000000,1024000000")
#include <cstdio>
#include <cstring>
#include <vector>
#include <iostream> using namespace std; vector <int> a[100001];
int n,m,now = 0,b[100001],xulie[100005],e[100005];
long long v[100005],dis[100005] ={0},maxsum[100005*4],lazy_tag[100005*4] = {0}; void dfs(int x,int pre)//深搜找dfs序
{
b[x] = ++now;//这是这个节点的开始位置
xulie[now] = x;//把这个序号的节点加入到序列中
int len = a[x].size();//这是x节点的出度
for(int i = 0;i <= len-1;i++)//vector是从0开始计数的
{
int y = a[x][i];//x的一个出度。
if (y!=pre)//如果没有往回走
{
dis[y] = dis[x]+v[y];//更新y节点它到0点的价值和
dfs(y,x);//继续dfs
}
}
e[x] = now;//这是节点的结束位置。b[x]..e[x]里面存的都是x的子树节点。
} long long max(long long a,long long b)
{
return a > b?a:b;
} void push_up(int rt)
{
maxsum[rt] = max(maxsum[rt<<1],maxsum[rt<<1|1]);//只要左右区间的最值取max就可以了。
} void build(int l,int r,int rt)
{
lazy_tag[rt] = maxsum[rt] = 0;//初始化懒惰标记和区间最值。
if (l == r)//如果左右指针相同
{
maxsum[rt] = dis[xulie[l]];//这个区间的最大值就是本身
return;//注意不要写成maxsum[rt] = dis[l],l是dfs序,xulie[l]才是这个dfs序所代表的点
}
int m = (l+r)>>1;//取得中点
build(l,m,rt<<1);//左递归建树
build(m+1,r,rt<<1|1);//右递归建树
push_up(rt);//根据儿子更新当前区间的最值。
} void push_down(int rt)
{
if (lazy_tag[rt]!=0) //如果懒惰标记标记的值不为0
{
maxsum[rt<<1]+=lazy_tag[rt];//更新儿子节点们的区间最值(记住是同等更新就好)
maxsum[rt<<1|1] += lazy_tag[rt];
lazy_tag[rt<<1]+=lazy_tag[rt];//更新儿子节点们的懒惰标记。
lazy_tag[rt<<1|1]+=lazy_tag[rt];
lazy_tag[rt]=0;//把当前节点的懒惰标记置为0;
}
} void updata(int l,int r,int num,int begin,int end,int rt) //所要修改的区间为l,r,当前节点rt的区间
{//为begin,end,然后要递增的值为num
if (l <= begin && end <= r) //如果当前节点所在的区间被所要修改的区间覆盖
{
maxsum[rt]+=num;//递增这个区间的最大值(因为是同等递增,所以直接把最大值+递增值就可以);
lazy_tag[rt] +=num;//标一个懒惰标记。
return;
}
push_down(rt);//处理懒惰标记。
int m = (begin+end)>>1;//取中点
if (l <= m)//往左递增值
updata(l,r,num,begin,m,rt<<1);
if (m < r)//往右递增值。
updata(l,r,num,m+1,end,rt<<1|1);
push_up(rt);//根据儿子节点更新当前节点。
} long long query(int l,int r,int begin,int end,int rt)//要找的区间为l,r,当前节点rt所代表的区间为begin,end;
{
if (l <= begin && end <= r)//如果当前节点所在的区间完全在要找的区间内
return maxsum[rt];//就直接返回这个区间的最大值
long long temp_ans= -10000000000;//这是用于取左右区间最大值的较大值。
int m = (begin+end)>>1;
push_down(rt);//先处理懒惰标记!!!!!!!!!!!!
if (l <= m)//如果左区间和所找的区间有交集,就看看左区间是否更大
temp_ans = max(temp_ans,query(l,r,begin,m,rt<<1));
if (m < r)//右区间同理
temp_ans = max(temp_ans,query(l,r,m+1,end,rt<<1|1));
return temp_ans;//返回左右区间的较大值即可。
} int main()
{
//freopen("rush.txt","r",stdin);
int t;
scanf("%d",&t);//输入数据的组数
for (int ii = 1;ii <= t;ii++)//不要写成while (t--)。。。这样写你要怎么输出呢。。。
{
now = 0;//dfs序号从0再重新开始
for (int i = 0;i <= 100000;i++)//初始化某个点的出度信息。
a[i].clear();
printf("Case #%d:\n",ii);//输出这是第几组数据
scanf("%d%d",&n,&m);//输入点和边的信息。
for (int i = 1;i <= n-1;i++)//输入n-1条边(就是树了!)
{
int x,y;
scanf("%d%d",&x,&y);//输入起点和终点
a[x].push_back(y);//建一条双向边。
a[y].push_back(x);
}
for (int i = 0;i <= n-1;i++)//输入点的初始价值。
scanf("%lld",&v[i]);
dis[0] = v[0];//从0到0就是0点的价值。
dfs(0,-1);//进行dfs,获取0到所有点的价值和。
build(1,n,1);//建树。
for (int i = 1;i <= m;i++)//输入m个操作。
{
int op,x,y;
scanf("%d",&op);//读取操作
if (op == 0)//如果要改变某个值
{
scanf("%d%d",&x,&y);
updata(b[x],e[x],y-v[x],1,n,1);//把替换操作改为递增操作。要进行修改的区间就是
v[x] = y;//x的子树所在的区间。
//v[x]下次可能还会用到所以要修改一下。
}
else
{
scanf("%d",&x);//读入x,表示要经过x(然后从x和x的子树节点里面选一个节点到那里)
cout << query(b[x],e[x],1,n,1) << endl;//在x的子树节点所在区间内找最大价值
}
}
}
return 0;
}
#pragma comment(linker, "/STACK:1024000000,1024000000")
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
using namespace std;
typedef long long ll;
const int maxn=1e5+7;
int in[maxn], out[maxn], xulie[maxn];
ll dis[maxn], v[maxn];
vector<int> g[maxn];
struct node{
int l, r;
ll maxx, lazy;
}t[maxn<<2];
int n, m, now;
void init()
{
now=0;
for(int i=0; i<=n; i++)
{
g[i].clear();
}
}
void dfs(int rt, int fa)
{
in[rt]=++now;
xulie[now]=rt;
int len=g[rt].size();
for(int i=0; i<len; i++)
{
int y=g[rt][i];
if(y==fa) continue;
dis[y]=dis[rt]+v[y];
dfs(y, rt);
}
out[rt]=now;
}
void pushup(int rt)
{
t[rt].maxx=max(t[rt<<1].maxx, t[rt<<1|1].maxx);
}
void build(int rt, int l, int r)
{
t[rt].l=l;
t[rt].r=r;
t[rt].lazy=0;
if(l==r)
{
t[rt].maxx=dis[xulie[l]];
return ;
}
int mid=(l+r)>>1;
build(rt<<1, l, mid);
build(rt<<1|1, mid+1, r);
pushup(rt);
}
void down(int rt)
{
if(t[rt].lazy!=0)
{
t[rt<<1].maxx+=t[rt].lazy;
t[rt<<1].lazy+=t[rt].lazy; t[rt<<1|1].maxx+=t[rt].lazy;
t[rt<<1|1].lazy+=t[rt].lazy; t[rt].lazy=0;
}
}
void update(int rt, int l, int r, int y)
{
if(l<=t[rt].l && t[rt].r <=r)
{
t[rt].maxx+=y;
t[rt].lazy+=y;
return ;
}
down(rt);
int mid=(t[rt].l+t[rt].r)>>1;
if(l<=mid) update(rt<<1, l, r, y);
if(r>mid) update(rt<<1|1, l, r, y);
pushup(rt);
}
ll query(int rt, int l, int r)
{
if(l<= t[rt].l && t[rt].r<=r)
{
return t[rt].maxx;
}
ll ans=-1e18;
int mid=(t[rt].l+t[rt].r)>>1;
down(rt);
if(l<=mid) ans=max(ans, query(rt<<1, l, r));
if(r>mid) ans=max(ans, query(rt<<1|1, l, r));
return ans;
}
int main()
{
int t;
scanf("%d", &t);
for(int ca=1; ca<=t; ca++)
{
scanf("%d%d", &n, &m);
init();
printf("Case #%d:\n", ca);
for(int i=1; i<n; i++)
{
int x, y;
scanf("%d%d", &x, &y);
g[x].push_back(y);
g[y].push_back(x);
}
for(int i=0; i<=n-1; i++)
{
scanf("%lld", &v[i]);
}
dis[0]=v[0];
dfs(0, -1);
build(1, 1, n);
for(int i=1; i<=m; i++)
{
int op, x, y;
scanf("%d", &op);
if(op==0)
{
scanf("%d%d", &x, &y);
update(1,in[x], out[x], y-v[x]);
v[x]=y;
}
else
{
scanf("%d", &x);
printf("%lld\n", query(1, in[x], out[x]));
}
}
}
return 0;
}