[BZOJ2834]回家的路

Description

Input

Output

Sample Input

2 1
1 2
1 1 2 2

Sample Output

5

 

思路还是很简单的，然而最短路打错各种对拍各种调了一早上

代码：

 #include<iostream>
 #include<cstdio>
 #include<cstring>
 #include<vector>
 #include<queue>
 #include<cstdlib>
 #include<map>
 #include<algorithm>
 #define M 300010
 using namespace std;
 vector<int>S1[M],S2[M];
 map<int,int>f[M];
 int n,m,num,cnt,x1,x2,y1,y2;
 int head[M],dis[M],vis[M],st[M],v[M];
 struct point{int to,next,dis;}e[M<<];
 void add(int from,int to,int dis) {
     e[++num].next=head[from];
     e[num].to=to;
     e[num].dis=dis;
     head[from]=num;
 }
 void SPFA(int s) {
     memset(dis,,sizeof(dis));
     queue<int>Q;dis[s]=;
     Q.push(s);
     while(!Q.empty()) {
         int x=Q.front();Q.pop();vis[x]=false;
         for(int i=head[x];i;i=e[i].next) {
             int to=e[i].to;
             if(dis[x]+e[i].dis<dis[to]) {
                 dis[to]=dis[x]+e[i].dis;
                 if(!vis[to]) vis[to]=true,Q.push(to);
             }
         }
     }
 }
 int get(int x,int y) {
     return *abs(v[st[x]]-v[st[y]]);
 }
 bool cmp(int x,int y) {
     return v[x]<v[y];
 }
 int main() {
     scanf("%d%d",&n,&m);
     for(int i=;i<=m;i++) {
         int x,y;scanf("%d%d",&x,&y);
         if(f[x][y]!=) continue;
         f[x][y]=++cnt;
         S1[x].push_back(cnt);v[cnt]=y;
         S2[y].push_back(++cnt);v[cnt]=x;
         add(cnt-,cnt,),add(cnt,cnt-,);
     }
     scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
     if(!f[x1][y1]) {
         f[x1][y1]=++cnt;
         S1[x1].push_back(cnt);v[cnt]=y1;
         S2[y1].push_back(++cnt);v[cnt]=x1;
         add(cnt-,cnt,),add(cnt,cnt-,);
     }
     else {
         int t=f[x1][y1];
         add(t+,t,),add(t,t+,);
     }
     if(!f[x2][y2]) {
         f[x2][y2]=++cnt;
         S1[x2].push_back(cnt);v[cnt]=y2;
         S2[y2].push_back(++cnt);v[cnt]=x2;
         add(cnt-,cnt,),add(cnt,cnt-,);
     }
     else {
         int t=f[x2][y2];
         add(t+,t,),add(t,t+,);
     }
     for(int i=;i<=(n<<);i++) {
         int tmp=;
         for(int j=;j<S1[i].size();j++) st[++tmp]=S1[i][j];
         if(tmp>) {
             sort(st+,st+tmp+,cmp);
             for(int j=;j<=tmp;j++) {
                 add(st[j-],st[j],get(j-,j));
                 add(st[j],st[j-],get(j,j-));
             }
         }
         tmp=;
         for(int j=;j<S2[i].size();j++) st[++tmp]=S2[i][j];
         if(tmp>) {
             sort(st+,st+tmp+,cmp);
             for(int j=;j<=tmp;j++) {
                 add(st[j-],st[j],get(j-,j));
                 add(st[j],st[j-],get(j,j-));
             }
         }
     }
     int s=f[x1][y1],t=f[x2][y2];
     SPFA(s);
     printf("%d\n",dis[t]>1e8?-:dis[t]);
     return ;
 }

HINT

N<=20000,M<=100000