最短路+离散化——BZOJ1967/Luogu2537 [AHOI2005]CROSS 穿越磁场

题面：BZOJ1967 Luogu2537
考虑到坐标范围非常大，但是矩阵数目很少
所以考虑离散化（离散化这个过程真的恶心）
离散化的时候要注意把矩形的边左右或者上下都给留下来，这样才能给机器人跑路
离散化之后呢跑一遍spfa，注意边界处理
就是说不能连续在一条边界上面走
至于两个矩形的边交叉处不必担心，因为往交叉处走的话必定不是最优的
最后输出终点的最短路径即可

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<queue>
#include<cstdlib>
#include<ctime>
#include<climits>
#include<string>
using namespace std;
const int dx[5]={0,1,-1,0,0};
const int dy[5]={0,0,0,1,-1};
struct ppap{int x,y;};
queue<ppap>q;
int sx,sy,zx,zy;
int n,x[1001],y[1001],c[1001],nn,mm;
int aa[1000005]={0},bb[1000005]={0};
int disx[10001],disy[10001],dist[301][301];
bool vis[301][301];int b[301][301];
int ans=0;
inline int read(){
    int k=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){k=k*10+ch-'0';ch=getchar();}
    return k*f;
}
inline int spfa(int sx,int sy){
    for(int i=1;i<=nn;i++)
        for(int j=1;j<=mm;j++)dist[i][j]=1e9;
    memset(vis,0,sizeof vis);
    dist[sx][sy]=0;ppap rp;rp.x=sx;rp.y=sy;q.push(rp);
    while(!q.empty()){
        ppap now=q.front();q.pop();
        int xx=now.x,yy=now.y;
        for(int i=1;i<=4;i++){
            int p,nx=xx+dx[i],ny=yy+dy[i];
            if(!b[nx][ny])p=0;
            else if(b[nx][ny]==b[xx][yy])continue;
            else p=1;
            if(dist[nx][ny]>dist[xx][yy]+p){
                dist[nx][ny]=dist[xx][yy]+p;
                if(!vis[nx][ny]){
                    ppap pr;pr.x=nx;pr.y=ny;q.push(pr);
                    vis[nx][ny]=1;
                }
            }
        }
        vis[xx][yy]=0;
    }
}
int main()
{
    n=read();
    for(int i=1;i<=n;i++){
        x[i]=read();y[i]=read();c[i]=read();
        aa[x[i]+1]=aa[x[i]]=aa[x[i]-1]=1;
        bb[y[i]]=bb[y[i]-1]=bb[y[i]+1]=1;
        aa[x[i]+c[i]+1]=aa[x[i]+c[i]]=aa[x[i]+c[i]-1]=1;
        bb[y[i]+c[i]+1]=bb[y[i]+c[i]]=bb[y[i]+c[i]-1]=1;
    }
    sx=read();sy=read();zx=read();zy=read();
    aa[sx]=aa[zx]=1;bb[sy]=bb[zy]=1;
    int p=0,q=0;
    for(int i=0;i<=1000001;i++)if(aa[i]){
        q++;aa[i]=q;if(q>1)disx[q]=i-p;
        p=i;
    }
    nn=q;p=q=0;
    for(int i=0;i<=1000001;i++)if(bb[i]){
        q++;bb[i]=q;if(q>1)disy[q]=i-p;
        p=i;
    }
    mm=q;
    sx=aa[sx];sy=bb[sy];zx=aa[zx];zy=bb[zy];
    for(int i=1;i<=n;i++){
        int xx=x[i],dx=x[i]+c[i],xy=y[i],dy=y[i]+c[i];
        xx=aa[xx];xy=bb[xy];dx=aa[dx];dy=bb[dy];
        for(int j=xx;j<=dx;j++)b[j][xy]=b[j][dy]=i;
        for(int j=xy;j<=dy;j++)b[xx][j]=b[dx][j]=i;
    }//离散化
    spfa(sx,sy);
    printf("%d",dist[zx][zy]);
    return 0;
}