bzoj3362&poj1984 Navigation Nightmare

http://www.elijahqi.win/archives/709
农夫约翰有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.

这个题和食物链那个带权并查集有些像 我们在连接的时候可以运用矢量计算来想想

要维护两个点之间的横向距离以及纵向距离

每次压缩距离的时候加上他们父亲的距离

其实关于这个矢量我在写程序的时候相当于强制要求儿子指向父亲

所以d(x|y)[y]-d(x|y)[x] 是最终的结果

#include<cstdio>
#include<algorithm>
#define N 44000
#define M 11000
inline int read(){
    int x=0;char ch=getchar();
    while (ch<'0'||ch>'9') ch=getchar();
    while (ch<='9'&&ch>='0'){x=x*10+ch-'0';ch=getchar();}
    return x;
}
struct node{
    int x,y,z;char d[2];
}data[N];
struct node1{
    int x,y,t,id;
}q[M];
inline bool cmp(node1 a,node1 b){
    return a.t<b.t;
}
int fa[N],dx[N],dy[N];
inline int find(int x){
    if (fa[x]==x) return x;
    int tmp=find(fa[x]);dx[x]+=dx[fa[x]];dy[x]+=dy[fa[x]];
    return fa[x]=tmp;
}
int ans[M],n,m,k;
int main(){
    freopen("bzoj3362.in","r",stdin);
    n=read();m=read();
    for (int i=1;i<=m;++i) data[i].x=read(),data[i].y=read(),data[i].z=read(),scanf("%s",data[i].d);k=read();
    for (int i=1;i<=k;++i) q[i].x=read(),q[i].y=read(),q[i].t=read(),q[i].id=i;
    std::sort(q+1,q+k+1,cmp);
    int now=1;
    for (int i=1;i<=n;++i) fa[i]=i;
    for (int i=1;i<=k;++i){
        while (now<=q[i].t){
            int xx=find(data[now].x),yy=find(data[now].y);
            if (xx==yy) continue;
            fa[xx]=yy;
            dx[xx]=dx[data[now].y]-dx[data[now].x];dy[xx]=dy[data[now].y]-dy[data[now].x];
            if (data[now].d[0]=='E') dx[xx]+=data[now].z;
            if (data[now].d[0]=='W') dx[xx]-=data[now].z;
            if (data[now].d[0]=='N') dy[xx]+=data[now].z;
            if (data[now].d[0]=='S') dy[xx]-=data[now].z;
            ++now;
        }
        int xx=find(q[i].x),yy=find(q[i].y);
        if (xx!=yy) {ans[q[i].id]=-1;continue;}
        ans[q[i].id]+=abs(dy[q[i].y]-dy[q[i].x])+abs(dx[q[i].y]-dx[q[i].x]);
    }
    for(int i=1;i<=k;++i) printf("%d\n",ans[i]);
    return 0;
}