hdu_3804_Query on a tree(树链剖分)

题目连接:http://acm.hdu.edu.cn/showproblem.php?pid=3804

题意:给你一棵树,然后给出树上边的价值,然后给出x,y,问从1到x的边上不超过y的最大值为多少

题解:这题如果直接写裸线段树来在线更新会T飞,正解是把边的值存下来,把询问也存下来,然后都按照价值从小到大排序,然后从最小的开始找,找之前把边的值插进线段树,这样就成了查询区间最大值了,很巧妙的方法。还有就是这题的内存卡的很紧,用前向星爆内存了,我没想通,用vector就过了。很尴尬

#include<cstdio>
#include<vector>
#include<algorithm>
#define F(i,a,b) for(int i=a;i<=b;++i)
#define root 1,n,1
#define ls l,m,rt<<1
#define rs m+1,r,rt<<1|1
using namespace std;
const int N=1e5+5;
vector<int>g[N];
struct vs{
	int val,u,v;
	bool operator<(const vs &b)const{return val<b.val;}
}va[N];
struct ser{
	int x,k,id;
	bool operator<(const ser &b)const{return k<b.k;}
}ask[N];
int t,n,m,x,y,c,tot,ans[N],sum[N*4],sz[N],fa[N],dep[N],hs[N],top[N],tid[N],idx;
//线段树部分
void update(int x,int k,int l,int r,int rt){
	if(l==r){sum[rt]=k;return;}
	int m=(l+r)>>1;
	if(x<=m)update(x,k,ls);
	else update(x,k,rs);
	sum[rt]=max(sum[rt<<1],sum[rt<<1|1]);
}

int query(int L,int R,int l,int r,int rt){
	if(L<=l&&r<=R)return sum[rt];
	int m=(l+r)>>1,maxs=0;
	if(L<=m)maxs=max(query(L,R,ls),maxs);
	if(m<R)maxs=max(query(L,R,rs),maxs);
	return maxs;
}
//树链部分
void dfs1(int u,int pre){
	sz[u]=1,fa[u]=pre,dep[u]=dep[pre]+1,hs[u]=0;
	for(int i=0;i<g[u].size();i++){
		int vv=g[u][i];
		if(vv!=pre){
			dfs1(vv,u);
			if(sz[vv]>sz[hs[u]])hs[u]=vv;
			sz[u]+=sz[vv];
		}
	}
}

void dfs2(int u,int tp){
	tid[u]=++idx,top[u]=tp;
	if(hs[u])dfs2(hs[u],tp);
	for(int i=0;i<g[u].size();i++){
		int vv=g[u][i];
		if(vv!=fa[u]&&vv!=hs[u])dfs2(vv,vv);
	}
}

int findmax(int x){
	int fx=top[x],an=0;
	while(fx!=1)an=max(query(tid[fx],tid[x],root),an),x=fa[fx],fx=top[x];
	if(x==1)return an;
	return max(query(2,tid[x],root),an);
}

int main(){
	scanf("%d",&t);
	while(t--){
		scanf("%d",&n);
		F(i,1,n)g[i].clear();
		dep[0]=0,sz[0]=0;
		F(i,0,4*N-1)sum[i]=0;
		F(i,1,n-1)scanf("%d%d%d",&va[i].u,&va[i].v,&va[i].val),g[va[i].u].push_back(va[i].v),g[va[i].v].push_back(va[i].u);
		dfs1(1,0),idx=0,dfs2(1,1);
		scanf("%d",&m),tot=1;
		F(i,1,m)scanf("%d%d",&ask[i].x,&ask[i].k),ask[i].id=i;
		sort(va+1,va+n),sort(ask+1,ask+1+m);
		F(i,1,m){
			while(tot<n&&va[tot].val<=ask[i].k){
				if(dep[va[tot].u]>dep[va[tot].v])c=va[tot].u;else c=va[tot].v;
				update(tid[c],va[tot].val,root),tot++;
			}
			ans[ask[i].id]=findmax(ask[i].x);
		}
		F(i,1,m)printf("%d\n",ans[i]==0?-1:ans[i]);
	}
	return 0;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 2
    评论
评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值