CF100015C

主要找到环上任意一条边,有比较dis(u,v),dis(u,a)+w+dis(b,v),dis(u,b)+w+dis(a,u)
然后,然后没了
lca求dis(u,v);dis(u,v)=dis[u]+dis[v]-2*dis[lca(u,v)];
我加油慢慢写
咳咳,我可真是写了半天才写过了样例,lca经测试无误。
可是啊

哭了...
还是附上代码

#include<iostream>
#include<cstring>
#include<algorithm>
#include<cstdio>
using namespace std;
struct node{
	int u,v,w,nxt;
}e[1000101];
int fir[1000101],cnt=0,top[1000101];
void add(int x,int y,int z){
	e[++cnt].nxt=fir[x];e[cnt].u=x;e[cnt].v=y;e[cnt].w=z;fir[x]=cnt;
}
int dis[1000101]={0};
int vis[1000101]={0},id[100101],tot=0,p=0,a=0,b=0,ww,siz[1000101],dep[1000101],son[1000101]={0},fa[1000101];
void dfs(int x,int f){
	for(int i=fir[x];i;i=e[i].nxt){
		int vv=e[i].v;//cout<<vv<<' '<<x<<' '<<i<<endl;
		if(vv==f)continue;
		if(vis[vv]==1&&a==0&&b==0){
		p=i;a=e[i].u;b=e[i].v;ww=e[i].w;e[i].w=-1;return;
		}
		if(!vis[vv]){
		vis[vv]=1;	dfs(vv,x);
		}
	}
}
void dfs1(int x,int f){
	siz[x]=1;dep[x]=dep[f]+1;fa[x]=f;
	for(int i=fir[x];i;i=e[i].nxt){
		int vv=e[i].v;//cout<<vv<<' '<<x<<' '<<i<<endl;
		if((a==x&&vv==b)||(vv==a&&x==b))e[i].w=-1;
		if(vv==f||(a==x&&vv==b)||(vv==a&&x==b))continue;
		if(!vis[vv]&&e[i].w!=-1){
		vis[vv]=1;dis[vv]=dis[x]+e[i].w;dfs1(vv,x);
		}
		siz[x]+=siz[vv];
		if(siz[vv]>siz[son[x]])son[x]=vv;
	}
}
void dfs2(int x,int topp){
    top[x]=topp;
    if(son[x])dfs2(son[x],topp);
    for(int i=fir[x];i;i=e[i].nxt){
        int vv=e[i].v;
        if(vv==fa[x]||vv==son[x]||(a==x&&vv==b)||(vv==a&&x==b))continue;
        dfs2(vv,vv);
    }
}
int lca(int x,int y){
	while(top[x]!=top[y]){
		 if(dep[top[x]]>=dep[top[y]])x=fa[top[x]];
            else y=fa[top[y]];
        }
		return dep[x]<dep[y]?x:y;
}
int diss(int x,int y,int zz){
	return dis[x]+dis[y]-2*dis[zz];
}
int n,m;
int main(){
	scanf("%d",&n);
	for(int i=1;i<=n;i++){
		int x,y,z;
		scanf("%d%d%d",&x,&y,&z);
		add(x,y,z);add(y,x,z);
	}
	dfs(0,0);
	//cout<<a<<' '<<b<<endl;
	memset(vis,0,sizeof(vis));siz[0]=0;dep[0]=0;
	dfs1(0,0);
	dfs2(0,0);
	int q;
	cin>>q;
	//for(int i=1;i<=n;i++)cout<<dis[i]<<' ';
	while(q--){
		int x,y;
		scanf("%d%d",&x,&y);
		int zz=lca(x,y),tt=lca(x,a),kk=lca(b,y);
		int ans=min(diss(x,y,zz),diss(x,a,tt)+ww+diss(b,y,kk));
		tt=lca(x,b);kk=lca(a,y);
		ans=min(ans,diss(x,b,tt)+ww+diss(a,y,kk));
		printf("%d\n",ans);
	}
	int lla;
	cin>>lla;
	return 0;
}

我觉得,照理,树剖的lca是跑得草鸡快的,至少比神马倍增跑得快
可是啊...不想说了

posted @ 2018-09-25 16:11  lisuier  阅读(136)  评论(0编辑  收藏  举报