LOJ #2610. 「NOIP2013」货车运输

题目描述

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

输入格式

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

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

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

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

输出格式

输出共有 qqq 行，每行一个整数，表示对于每一辆货车，它的最大载重是多少。如果货车不能到达目的地，输出 −1-1−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<10000 < n < 1\,0000<n<1000，0<m<100000 < m < 10\,0000<m<10000，0<q<10000 < q < 1\,0000<q<1000；

对于 60% 的数据，0<n<10000 < n < 1\,0000<n<1000，0<m<500000 < m < 50\,0000<m<50000，0<q<10000 < q < 1\,0000<q<1000；

对于 100% 的数据，0<n<100000 < n < 10\,0000<n<10000，0<m<500000 < m < 50\,0000<m<50000，0<q<300000 < q < 30\,0000<q<30000，0≤z≤1000000 \leq z \leq 100\,0000≤z≤100000。

 

 

kruskal+LCA裸题

Prim改成kruskal 就过了

 

#include<iostream>
#include<cstdio>
#include<algorithm>
using namespace std;
inline int read()
{
	int ret=0;
	char ch=getchar();
	while(ch<'0'||ch>'9') ch=getchar();
	while(ch>='0'&&ch<='9')
		ret=(ret<<1)+(ret<<3)+ch-'0',ch=getchar();
	return ret;
}

int n,m,cnt,cnt1,num;
const int N=1e5+5;
int to[N],nxt[N],w[N],he[N],f[N],gen[N];
int dep[N],lg[N],fi[N][20],fw[N][20];

void dfs(int fa,int u)
{
	dep[u]=(fa==0)?0:dep[fa]+1;
	fi[u][0]=fa; 
	for(int i=1;i<=lg[dep[u]];i++)
		fi[u][i]=fi[fi[u][i-1]][i-1],
		fw[u][i]=min(fw[u][i-1],fw[fi[u][i-1]][i-1]);
	for(int e=he[u];e;e=nxt[e])
	{
		int v=to[e];
		if(fa!=v)
			fw[v][0]=w[e],dfs(u,v);
	}	
}

inline int get(int u,int v)
{
	if(dep[u]>dep[v]) swap(u,v);
	int ret=1e9;
	while(dep[u]<dep[v])
		ret=min(ret,fw[v][lg[dep[u]-dep[v]]]),
		v=fi[v][lg[dep[u]-dep[v]]];	
	for(int i=lg[dep[u]];i>=0;i--)
		if(fi[u][i]!=fi[v][i])
			ret=min(ret,min(fw[v][i],fw[u][i])),
			v=fi[v][i],u=fi[u][i];
	if(u!=v) 
		ret=min(ret,min(fw[v][0],fw[u][0])),
		u=fi[u][0],v=fi[v][0];
	return ret;
}

inline void add(int u,int v,int k)
{
	to[++cnt]=v;
	nxt[cnt]=he[u];
	w[cnt]=k;
	he[u]=cnt;
}

inline int find(int x)
{
	return (x==f[x])?x:f[x]=find(f[x]);
}

struct NA{
	int u,v,w;
}e[N];

bool cmp(NA i,NA j)
{
	return i.w>j.w;
}
int main()
{
	n=read(),m=read();
	for(int i=1;i<=n;i++) f[i]=i;
	for(int i=1;i<=m;i++)
		e[i].u=read(),e[i].v=read(),e[i].w=read();
	sort(e+1,e+m+1,cmp);
	for(int i=1;i<=m;i++)
	{
		int fa1=find(e[i].u),fa2=find(e[i].v);
		if(fa1!=fa2) 
			f[fa1]=fa2,add(e[i].u,e[i].v,e[i].w),add(e[i].v,e[i].u,e[i].w);
	}
	
	for(int i=1;i<=n;i++)
		if(f[i]==i) gen[++num]=i;
	
	lg[0]=-1;
	for(int i=1;i<=n;i++)
		lg[i]=lg[i>>1]+1;
	//	printf("Gj love zyx");
	for(int r=1;r<=num;r++)
		dfs(0,gen[r]);
	//	printf("-1");
	int t=read(); 
	while(t--)
	{
		int u=read(),v=read();
		if(find(u)!=find(v)) 
		{
			printf("-1\n"); continue;
		}
		printf("%d\n",get(u,v));
	}
	return 0;
}