TLE求神犇帮忙
查看原帖
TLE求神犇帮忙
545918
hgzxgzx楼主2022/7/23 10:01
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<vector>
#include<set>
#include<queue>
#include<map>
#define File(x) freopen(""#x".in","r",stdin)
#define rep(i,l,r) for(int i=l; i<r; i++)
#define drep(i,r,l) for(int i=r; i>l; i--)
#define pii pair<int,int>
#define mdp make_pair
using namespace std;
typedef long long LL;
const int N=2e5+5e4+10;
const int W=0x7f7f7f7f;
vector<pii > e1[N],e2[N];// vector存图 
int fa[N][22],dep[N];
int times,dfn[N],s[N],id[N];
int sta[N],top;
bool mark[N];
LL dp[N];
template <class T>
void R(T &a){
	a=0;
	int f=1;char c=getchar();
	while(c<'0'||c>'9') {if(c=='-') f=-1;c=getchar();}
	while(c>='0'&&c<='9') {a=(a<<1)+(a<<3)+c-'0';c=getchar();}
	a*=f;
}
void dfs(int u,int ba){
	dfn[u]=++times;
	for(int i=0; i<e1[u].size(); i++){
		int v=e1[u][i].first;
		if(v!=ba){
			dep[v]=dep[u]+1;
			fa[v][0]=u;
			dfs(v,u);
		}
	}
}
 
inline bool comp(const int &a,const int &b){
	return dfn[a]<dfn[b];
}
int Lca(int x,int y){
	if(dep[x]<dep[y]) swap(x,y);
	drep(i,21,-1) if(dep[fa[x][i]]>=dep[y]) x=fa[x][i];
	if(x==y) return x;
	drep(i,21,-1) if(fa[x][i]!=fa[y][i]) x=fa[x][i],y=fa[y][i];
	return fa[x][0];
}
int obt_w(int u,int v){// 因为虚树中一些点之间的边在原图不存在,所以贪心地从 u 到 v 的一条路径上选最小边权 
	if(dep[u]<dep[v]) swap(u,v);
	int ans=W;
	while(u!=v){
		ans=min(ans,e1[fa[u][0]][id[u]].second);
		u=fa[u][0];
	}
	return ans;
}
void Build_tree(int k){
	sort(s+1,s+k+1,comp);
	sta[top=1]=1;
	rep(i,1,k+1){
		/*
		建树思路:
		如果单调栈栈顶和本节点在一条链上,直接加入
		否则判断单调栈栈顶和本节点的 lca 是否已入栈,不断从栈顶弹出元素直到次大节点的深度小于等于 lca 的深度,如果不是 lca 就加入 lca 
		*/ 
		int lca=Lca(sta[top],s[i]);
		if(lca==sta[top]) sta[++top]=s[i];
		else{
			while(dep[sta[top-1]]>dep[lca]){
				e2[sta[top-1]].push_back(mdp(sta[top],obt_w(sta[top],sta[top-1])));
				top--;
			}
			if(sta[top-1]==lca){
				e2[lca].push_back(mdp(sta[top],obt_w(sta[top],lca)));
				top--;
			}
			else{
				e2[lca].push_back(mdp(sta[top],obt_w(sta[top],lca)));
				sta[top]=lca;
			}
			sta[++top]=s[i];
		}
	}
	while(top>1){
		e2[sta[top-1]].push_back(mdp(sta[top],obt_w(sta[top],sta[top-1])));
		top--;
	}
}
void run_dp(int u){
	for(int i=0; i<e2[u].size(); i++){
		int v=e2[u][i].first;
		run_dp(v);
		if(!mark[v]) dp[u]+=min(dp[v],(LL)e2[u][i].second);
		else dp[u]+=e2[u][i].second;
	}
}
void clear(int num){
	rep(i,1,n+1) e2[i].clear();// 删除一颗树(只想出来这么做) 
	memset(mark,0,sizeof(mark));// 记录本次任务的关键点 
}
void print(int n){
	rep(i,1,n+1){
		rep(j,0,e2[i].size()) printf("%d%c",e2[i][j].first,j==e2[i].size()-1?'\n':' ');
	}	
}
int main(){
	//File(ts);
	int n;
	R(n);
	int u,v,w;
	rep(i,1,n){
		R(u),R(v),R(w);
		e1[u].push_back(mdp(v,w));
		e1[v].push_back(mdp(u,w));
	}
	dep[1]=1;
	dfs(1,0);// 求出dfs序 
	rep(i,1,n){
		int cnt=0;
		rep(j,0,e1[i].size()) if(i==fa[e1[i][j].first][0]) id[e1[i][j].first]=j;// 记录每个节点在其父亲节点上的位置 
	}
	rep(i,1,21) rep(j,1,n+1) fa[j][i]=fa[fa[j][i-1]][i-1];
	int m;
	R(m);
	int num,p;
	rep(i,1,m+1){
		memset(dp,0,sizeof(dp));
		R(num);
		rep(j,1,num+1) R(s[j]),mark[s[j]]=true;
		Build_tree(num); 
		run_dp(1);
		clear(n);
		printf("%lld\n",dp[1]);
	}
	return 0;
}
2022/7/23 10:01
加载中...