洛谷 P1967 [NOIP2013 D1T3] 货车运输

题目描述

A 国有 n 座城市，编号从 1 到 n，城市之间有 m 条双向道路。每一条道路对车辆都有重量限制，简称限重。现在有 q 辆货车在运输货物， 司机们想知道每辆车在不超过车辆限重的情况下，最多能运多重的货物。

输入输出格式

输入格式：

输入文件名为 truck.in。

输入文件第一行有两个用一个空格隔开的整数 n，m，表示 A 国有 n 座城市和 m 条道

路。 接下来 m 行每行 3 个整数 x、 y、 z，每两个整数之间用一个空格隔开，表示从 x 号城市到 y 号城市有一条限重为 z 的道路。意：x 不等于 y，两座城市之间可能有多条道路。

接下来一行有一个整数 q，表示有 q 辆货车需要运货。

接下来 q 行，每行两个整数 x、y，之间用一个空格隔开，表示一辆货车需要从 x 城市运输货物到 y 城市，注意：x 不等于 y。

输出格式：

输出文件名为 truck.out。

输出共有 q 行，每行一个整数，表示对于每一辆货车，它的最大载重是多少。如果货

车不能到达目的地，输出-1。

输入输出样例

输入样例#1：

4 3
1 2 4
2 3 3
3 1 1
3
1 3
1 4
1 3

输出样例#1：

3
-1
3

说明

对于 30%的数据，0 < n < 1,000，0 < m < 10,000，0 < q< 1,000； 对于 60%的数据，0 < n < 1,000，0 < m < 50,000，0 < q< 1,000； 对于 100%的数据，0 < n < 10,000，0 < m < 50,000，0 < q< 30,000，0 ≤ z ≤ 100,000。

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

最大生成树+lca~

两点之间可能有多条路径，但只要保留最大的就可以了，所以先用最大生成树（森林？）重建一下图，再跑lca，计算faa数组的时候同时记录一下该点到2^i节点上的最小边权，这样可以很方便地求出总的最小值。每次输入的时候先判断一下两点在不在一棵树上，再输出结果~

（程序好像只能在5.9.2以上运行，低版本会输出error 0，并不明白是为什么……但是在洛谷上是A的~）

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;

int n,m,q,fa[10001],w[50001],v[50001],fi[10001],ne[50001],cnt,tot,deep[10001],faa[10001][20],du[10001][20],xx,yy;
bool b[10001];

struct node{
	int x,y,v;
}ro[50001];

void add(int u,int vv,int z)
{
	w[++cnt]=vv;ne[cnt]=fi[u];fi[u]=cnt;v[cnt]=z;
}

bool cmp(node u,node v)
{
	return u.v>v.v;
}

int findd(int u)
{
	return fa[u]==u ? u:fa[u]=findd(fa[u]);
}

void dfs(int u)
{
	b[u]=1;
	for(int i=1;i<=16;i++)
	{
		if(deep[u]<(1<<i)) break;
		faa[u][i]=faa[faa[u][i-1]][i-1];
		du[u][i]=min(du[u][i-1],du[faa[u][i-1]][i-1]);
	}
	for(int i=fi[u];i;i=ne[i])
	  if(!b[w[i]])
	  {
	  	faa[w[i]][0]=u;du[w[i]][0]=v[i];
	  	deep[w[i]]=deep[u]+1;dfs(w[i]);
	  }
}

int lca(int u,int v)
{
	if(deep[u]<deep[v]) swap(u,v);
	int kkz=deep[u]-deep[v];
	for(int i=0;i<=16;i++)
	  if(kkz&(1<<i)) u=faa[u][i];
	for(int i=16;i>=0;i--)
	  if(faa[u][i]!=faa[v][i]) u=faa[u][i],v=faa[v][i];
	if(u==v) return u;
	return faa[u][0];
}

int ans(int u,int v)
{
	int minn=999999999;
	int kkz=deep[u]-deep[v];
	for(int i=0;i<=16;i++)
	  if(kkz&(1<<i))
	  {
	  	minn=min(minn,du[u][i]);
	  	u=faa[u][i];
	  }
	return minn;
}

int main()
{
	scanf("%d%d",&n,&m);
	memset(du,127/3,sizeof(du));
	for(int i=1;i<=n;i++) fa[i]=i;
	for(int i=1;i<=m;i++) scanf("%d%d%d",&ro[i].x,&ro[i].y,&ro[i].v);
	sort(ro+1,ro+m+1,cmp);
	for(int i=1;i<=m;i++)
	  if(findd(ro[i].x)!=findd(ro[i].y))
	  {
	  	fa[findd(ro[i].x)]=findd(ro[i].y);tot++;
	  	add(ro[i].x,ro[i].y,ro[i].v);
	  	add(ro[i].y,ro[i].x,ro[i].v);
	  	if(tot==n-1) break;
	  }
	for(int i=1;i<=n;i++)
	  if(!b[i]) dfs(i);
	scanf("%d",&q);
	while(q--)
	{
		scanf("%d%d",&xx,&yy);
		if(findd(xx)!=findd(yy)) printf("-1\n");
		else
		{
			int kkz=lca(xx,yy);
			printf("%d\n",min(ans(xx,kkz),ans(yy,kkz))); 
		}
	}
	return 0;
}