• [20191004机房测试] ZGY的早餐


    ZGY 每天早上要从宿舍走路到机房,顺便从学校小卖部购买早饭,当然机智的 ZGY 一定会走最短路
    学校的路可以看成一无向联通张图,图上有 n 个点,m 条边,每一个点都有一个唯一的编号 1~n
    每一条边有一个边权,表示两个点之间的距离,ZGY 的宿舍在 S 点,机房在 T点,而小卖部在 H 点
    现在 ZGY 想知道从宿舍经过小卖部到达机房的最短距离
    不过因为在这个世界上有 Q个 ZGY,所以你必须回答 Q 个问题
    

    很棒的数据分治题
    读入里面说了会读入测试点编号……
    其实是很明显的暗示了……

    一半的数据保证是连通图,1 ≤ n ≤ 200,1 ≤ Q ≤ (10^5)
    这么小的n,直接Floyd跑一遍就好了,相当于Floyd的板子题

    另一半的数据保证图连通,无环,1 ≤ n ≤ (10^5),1 ≤ Q ≤ (10^5)
    联通图无环……那不就是一棵树吗!
    dfs预处理出到根节点的距离和LCA,每次询问输出就好

    代码:

    #include<bits/stdc++.h>
    #define ll long long
    #define inf 1000000000000000000
    using namespace std;
    
    int t,n,m,u,v,w;
    int q,S,H,T;
    
    struct Edge
    {
    	int next,to,dis;
    }edge[1000005];
    int cnt=0,head[100005];
    
    inline void add_edge(int from,int to,int dis)
    {
    	edge[++cnt].next=head[from];
    	edge[cnt].to=to;
    	edge[cnt].dis=dis;
    	head[from]=cnt;
    }
    
    template<class T>inline void read(T &res)
    {
    	char c;T flag=1;
    	while((c=getchar())<'0'||c>'9')if(c=='-')flag=-1;res=c-'0';
    	while((c=getchar())>='0'&&c<='9')res=res*10+c-'0';res*=flag;
    }
    
    ll dis[205][205];
    void Floyd()
    {
    	for(register int k=1;k<=n;++k)
    		for(register int i=1;i<=n;++i)
    			for(register int j=1;j<=n;++j)
    				dis[i][j]=min(dis[i][j],dis[i][k]+dis[k][j]);
    }
    
    ll dep[100005];
    ll ldep[100005];
    int anc[400005][25];
    void dfs(int u,int fa)
    {
    	anc[u][0]=fa;
    	for(int i=1;i<=20;i++) anc[u][i]=anc[anc[u][i-1]][i-1];
    	for(register int i=head[u];i;i=edge[i].next)
    	{
    		int v=edge[i].to;
    		if(v==fa) continue;
    		if(dep[v]>dep[u]+edge[i].dis)
    		{
    			dep[v]=dep[u]+edge[i].dis;
    			ldep[v]=ldep[u]+1;
    			dfs(v,u);
    		}
    	}
    }
    
    int LCA(int x,int y)
    {
    	if(ldep[x]<ldep[y]) swap(x,y);
    	int deltaDep=ldep[x]-ldep[y];
    	for(int i=0;deltaDep;deltaDep>>=1,i++)
    		if(deltaDep&1) x=anc[x][i];
    	if(x==y) return x;
    	for(int i=20;anc[x][0]!=anc[y][0];i--)
    	{
    		if(anc[x][i]!=anc[y][i])
    		{
    			x=anc[x][i];
    			y=anc[y][i];
    		}
    	}
    	return anc[x][0];
    }
    
    ll Dis(int u,int v) {return dep[u]+dep[v]-2LL*dep[LCA(u,v)];}
    
    int main()
    {
    	freopen("mindis.in","r",stdin);
    	freopen("mindis.out","w",stdout);
    	read(t);read(n);read(m);
    	if(1<=t&&t<=5)
    	{
    		for(register int i=1;i<=n;++i)
    			for(register int j=1;j<=n;++j)
    				if(i!=j) dis[i][j]=inf;
    		for(register int i=1;i<=m;++i)
    		{
    			read(u);read(v);read(w);
    			dis[u][v]=dis[v][u]=w;
    		}
    		Floyd();
    		read(q);
    		while(q--)
    		{
    			read(S);read(H);read(T);
    			printf("%lld
    ",dis[S][H]+dis[H][T]);
    		}
    	}
    	else if(6<=t&&t<=10)
    	{
    		for(register int i=1;i<=m;++i)
    		{
    			read(u);read(v);read(w);
    			add_edge(u,v,w);
    			add_edge(v,u,w);
    		}
    		for(register int i=1;i<=n;++i) dep[i]=ldep[i]=inf;
    		dep[1]=0;
    		ldep[1]=1;
    		dfs(1,0);
    		read(q);
    		while(q--)
    		{
    			read(S);read(H);read(T);
    			printf("%lld
    ",Dis(S,H)+Dis(H,T));
    		}
    	}
    	return 0;
    }
    
  • 相关阅读:
    MIPS笔记
    花生壳动态域名解析工具原理
    SEE MIPS RUN 第六章 内存管理与TLB
    C/C++动态内存创建与内存管理
    ngclass 用法
    看了一个烟花的html作品 引用:http://www.w3cfuns.com/blog54440495404365.html
    ngclip angualr 的copy功能
    学习技术的方法
    restful restAPI 的定义方式
    我与计算机
  • 原文地址:https://www.cnblogs.com/tqr06/p/11622898.html
Copyright © 2020-2023  润新知