#代码是最为耐心、最能忍耐和最令人愉快的伙伴,在任何艰难困苦的时刻,它都不会抛弃你。#
今天每日一句好有趣啊。。感觉WA的代码只会一点点掏空我。。。
不说了,今天写一个LCA(最近公共祖先)就是关于一个跳跳跳的东西???
洛谷题面_(:зゝ∠)_ noip2013day1第三题 货车运输???(老司机【滑稽】??)
A 国有 n 座城市,编号从 1 到 n,城市之间有 m 条双向道路。每一条道路对车辆都有重量限制,简称限重。现在有 q 辆货车在运输货物, 老司机们想知道每辆车在不超过车辆限重的情况下,最多能运多重的货物。
【输入】
第一行有两个用一个空格隔开的整数n,m,表示A国有n座城市和m条道路。接下来m行每行3个整数x、y、z,每两个整数之间用一个空格隔开,表示从x号城市到y号城市有一条限重为z的道路。注意:x不等于y,两座城市之间可能有多条道路。接下来一行有一个整数q,表示有q辆货车需要运货。接下来q行,每行两个整数x、y,之间用一个空格隔开,表示一辆货车需要从x城市运输货物到y城市,注意:x不等于y。
【输出】
共有q行,每行一个整数,表示对于每一辆货车,它的最大载重是多少。如果货车不能到达目的地,输出-1。
所以是一道奇怪的题,首先是双向道路也就是无向图,addedge的时候要写两个。可以看到如果货车不能到达目的地的话就输出-1,所以可以骗分并查集来判断连通性。
所以LCA是啥?
LCA就是(Least Common Ancestors)最近公共祖先
比如像
比如图中21和15的LCA就是1
计算LCA分为5部分
1、把低的那个(就是深度较大)的那个跳到另一个点的深度(如果这一步完了重合也就是高的是LCA)
2、这一步21跳到13以后,两个分别跳到6和7也就是LCA下面一个的点(这里是两步)
3、最后跳一步合并维护(所以这里也是两步)
如果货车想运输尽可能多的东西,就要走限重尽可能大的路,所以说构成的树是最大生成树,先从大到小排序,然后预处理。然后用并查集来维护是否联通,如果不连通直接输出-1,然后用树上倍增求出答案,详见代码。
#include<iostream>#include<cstdio>
#include<algorithm>
#include<climits>
using namespace std;
int n,m;
struct road
{
int u,v,w;
}c[50010];
bool cmp(road x,road y)//存图
{
return x.w>y.w;
}
int p[10010];
int find(int x)//并查集
{
if(p[x]==0) return x;
return p[x]=find(p[x]);
}
struct node //内存池
{
int v,w;
int next;
}pool[30010];
int h[10010],cnt;
void addedge(int u,int v,int w) //建树
{
cnt++; pool[cnt].v=v; pool[cnt].w=w; pool[cnt].next=h[u]; h[u]=cnt;
cnt++; pool[cnt].v=u; pool[cnt].w=w; pool[cnt].next=h[v]; h[v]=cnt;
}
int dep[10010],flag[10010],fa[10010][15],mx[10010][15];
void dfs(int u)// 预处理
{
flag[u]=1;
int v;
for(int i=h[u];i;i=pool[i].next)
if(flag[v=pool[i].v]==0)
{
dep[v]=dep[u]+1;
fa[v][0]=u;
mx[v][0]=pool[i].w;
for(int j=1;fa[v][j]=fa[fa [v][j-1] ][j-1];j++)
mx[v][j]=min(mx[v][j-1],mx[ fa[v][j-1] ][j-1]);
dfs(v);
}
}
int lca(int u,int v)
{
if(dep[u]<dep[v]) swap(u,v);
int ans=INT_MAX;
for(int i=14;i>=0;i--)
if(dep[fa[u][i]]>=dep[v])
{
ans=min(ans,mx[u][i]);
u=fa[u][i];
}
if(u==v) return ans;
for(int i=14;i>=0;i--)
if(fa[u][i]!=fa[v][i])
{
ans=min(ans,mx[u][i]);
ans=min(ans,mx[v][i]);
u=fa[u][i];
v=fa[v][i];
}
return min(ans,min(mx[u][0],mx[v][0]));
}
int q;
int main()
{
scanf("%d%d",&n,&m);
for(int i=1;i<=m;i++)
scanf("%d%d%d",&c[i].u,&c[i].v,&c[i].w);
sort(c+1,c+m+1,cmp);
int u,v,w;
for(int i=1;i<=m;i++)
if(find(u=c[i].u)!=find(v=c[i].v))
{
addedge(u,v,c[i].w);
p[find(u)]=find(v);
}
for(int i=1;i<=n;i++)
if(flag[i]==0)
{
dep[i]=1;
dfs(i);
}
scanf("%d",&q);
for(int i=1;i<=q;i++)
{
scanf("%d%d",&u,&v);
if(find(u)!=find(v)) printf("-1\n");
else printf("%d\n",lca(u,v));
}
return 0;
}