bzoj 3362: [Usaco2004 Feb]Navigation Nightmare 导航噩梦（加权并查集）

3362: [Usaco2004 Feb]Navigation Nightmare 导航噩梦

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 121  Solved: 63
[ Submit][ Status][ Discuss]

Description

    农夫约翰有N(2≤N≤40000)个农场，标号1到N,M(2≤M≤40000)条的不同的垂直或水

平的道路连结着农场，道路的长度不超过1000.这些农场的分布就像下面的地图一样，

图中农场用F1..F7表示, 每个农场最多能在东西南北四个方向连结4个不同的农场．此外，农场只处在道路的两端．道路不会交叉且每对农场间有且仅有一条路径．邻居鲍伯要约翰来导航，但约翰丢了农场的地图，他只得从电脑的备份中修复了．每一条道路的信息如下：

从农场23往南经距离10到达农场17

从农场1往东经距离7到达农场17

    当约翰重新获得这些数据时，他有时被的鲍伯的问题打断：“农场1到农场23的曼哈顿距离是多少？”所谓在(XI，Yi)和(X2，y2)之间的“曼哈顿距离”，就是lxl - X21+lyl - y21．如果已经有足够的信息，约翰就会回答这样的问题（在上例中答案是17），否则他会诚恳地抱歉并回答-1.

Input

    第1行：两个分开的整数N和M.

    第2到M+1行：每行包括4个分开的内容，F1，F2，三，D分别描述两个农场的编号，道路的长度，F1到F2的方向N，E，S，w．

    第M+2行：一个整数，K(1≤K≤10000)，表示问题个数．

    第M+3到M+K+2行：每行表示一个问题，由3部分组成：Fi，F2，，．其中Fi和F2表示两个被问及的农场．而/(1≤J≤M)表示问题提出的时刻．J为1时，表示得知信息1但未得知信息2时．

Output

    第1到K行：每行一个整数，回答问题．表示两个农场间的曼哈顿距离．不得而知则输出-1.

Sample Input

7 6
1 6 13 E
6 3 9 E
3 5 7 S
4 1 3 N
2 4 20 W
4 7 2 S
3
1 6 1
1 4 3
2 6 6

Sample Output

13
-1
10

HINT

   在时刻1，约翰知道1到6的距离为13;在时刻3，1到4的距离仍然不知道；在时刻6，位置6向

北3个距离，向西7个距离于位置2，所以距离为10.

Source

Green

[ Submit][ Status][ Discuss]

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#define N 100003
using namespace std;
int fa[N],disx[N],disy[N],ans[N];
int n,m,q;
struct data{
	int t,x,y,l,opt,mark;
}a[N],e[N];
int cmp(data a,data b)
{
	return a.t<b.t;
}
int find(int x)
{
	if (fa[x]==x) return x;
	int t=find(fa[x]);
	disx[x]+=disx[fa[x]];
	disy[x]+=disy[fa[x]];
	fa[x]=t;
	return fa[x];
}
int main()
{
	freopen("a.in","r",stdin);
	//freopen("my.out","w",stdout);
	scanf("%d%d",&n,&m);
	for (int i=1;i<=m;i++) 
	{
	 char s[10];
	 scanf("%d%d%d",&a[i].x,&a[i].y,&a[i].l,s);
	 char c=getchar();
	 while (c<'A'||c>'Z') c=getchar();
	 switch(c){
	 	case 'E':a[i].mark=1;break;
	 	case 'W':a[i].mark=2;break;
	 	case 'S':a[i].mark=3;break;
	 	case 'N':a[i].mark=4;break;
	 }
	 a[i].t=i;
	 a[i].opt=0;
    }
    scanf("%d",&q);
    for (int i=1;i<=q;i++) {
    	scanf("%d%d%d",&e[i].x,&e[i].y,&e[i].t);
    	e[i].l=i;
	}
	for (int i=1;i<=n;i++) fa[i]=i;
	sort(e+1,e+q+1,cmp); int now=1;
	for (int i=1;i<=q;i++) {
		//cout<<e[i].x<<" "<<e[i].y<<endl;
		while (now<=m&&now<=e[i].t) {
			int r1=find(a[now].x); int r2=find(a[now].y);
		//	cout<<a[now].x<<" "<<a[now].y<<"!"<<endl;
			if (r1!=r2) {	
				fa[r2]=r1; int ans=1;
				if (a[now].mark==2||a[now].mark==3) ans=-1;
				if (a[now].mark<=2) disx[r2]=disx[a[now].x]-disx[a[now].y]-ans*a[now].l,disy[r2]=disy[a[now].x]-disy[a[now].y];
				else disy[r2]=disy[a[now].x]-disy[a[now].y]-ans*a[now].l,disx[r2]=disx[a[now].x]-disx[a[now].y];
		    }
			now++; 
		}
		int r1=find(e[i].x); int r2=find(e[i].y);
		if (r1!=r2) ans[e[i].l]=-1;
		else {
			int x=disx[e[i].x]-disx[e[i].y];
			int y=disy[e[i].x]-disy[e[i].y];
			ans[e[i].l]=abs(x)+abs(y);
		}
	}
	for (int i=1;i<=q;i++) printf("%d\n",ans[i]);
}