• 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
    ",ans);
    	}
    	int lla;
    	cin>>lla;
    	return 0;
    }
    

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

  • 相关阅读:
    第8章 降维
    第7章 集成学习和随机森林
    JS利用async、await处理少见的登录业务逻辑
    SQL SERVER 实现多行转多列
    Mysql函数----控制流函数介绍
    继承----静态代码快、构造方法、代码块、普通方法的执行顺序
    RBAC----基于角色的访问权限控制
    秋招-----思特沃克视频面试总结
    tomcat启动失败的三种方法
    索引之----mysql联合索引
  • 原文地址:https://www.cnblogs.com/lisuier/p/9700073.html
Copyright © 2020-2023  润新知