【BZOJ1967】[AHOI2005]穿越磁场（最短路）
题面
题解
一个显然的思路是这样的，我们的正方形的边长把整个平面割成了若干块，显然每个联通块都可以看着做一个点，那么接下来只需要把所有的有相邻边的联通块全部连上一条长度为(1)的边表示如果要从这个联通块到达另外一个联通块，需要穿过这一条边，那么最终的答案就变成了起点和终点所在联通块的最短路。因为找联通块不好搞，其实离散之后再当成网格图做就好了。
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<queue>
using namespace std;
#define ll long long
#define MAX 250
inline int read()
{
int x=0;bool t=false;char ch=getchar();
while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
if(ch=='-')t=true,ch=getchar();
while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
return t?-x:x;
}
struct Matrix{int x1,x2,y1,y2;}p[MAX<<1];
int n,Sx[MAX<<3],Sy[MAX<<3],tx,ty;
int Bx,By,Ex,Ey;
bool visx[MAX<<1][MAX<<1],visy[MAX<<1][MAX<<1];
struct Line{int v,next,w;}e[MAX*MAX*10];
int h[MAX*MAX<<2],cnt=1;
inline void Add(int u,int v,int w){e[cnt]=(Line){v,h[u],w};h[u]=cnt++;}
int d[4][2]={1,0,-1,0,0,1,0,-1};
int id(int x,int y){return x*ty+y-ty;}
int dis[MAX*MAX<<1];bool vis[MAX*MAX<<1];
void SPFA()
{
queue<int> Q;Q.push(id(Bx,By));
memset(dis,63,sizeof(dis));dis[id(Bx,By)]=0;
while(!Q.empty())
{
int u=Q.front();Q.pop();
for(int i=h[u];i;i=e[i].next)
{
int v=e[i].v;
if(dis[v]>dis[u]+e[i].w)
{
dis[v]=dis[u]+e[i].w;
if(!vis[v])vis[v]=true,Q.push(v);
}
}
vis[u]=false;
}
}
int main()
{
n=read();
for(int i=1;i<=n;++i)
{
int x=read(),y=read(),c=read();
p[i]=(Matrix){x,x+c,y+c,y};
}
Bx=read();By=read();Ex=read();Ey=read();
Sx[++tx]=Bx;Sx[++tx]=Ex;Sy[++ty]=By;Sy[++ty]=Ey;
for(int i=1;i<=n;++i)Sx[++tx]=p[i].x1,Sx[++tx]=p[i].x2;
for(int i=1;i<=n;++i)Sy[++ty]=p[i].y1,Sy[++ty]=p[i].y2;
sort(&Sx[1],&Sx[tx+1]);sort(&Sy[1],&Sy[ty+1]);
tx=unique(&Sx[1],&Sx[tx+1])-Sx-1;ty=unique(&Sy[1],&Sy[ty+1])-Sy-1;
for(int i=2;i<=tx;++i)if(Sx[i]-Sx[i-1]>1)Sx[++tx]=Sx[i]-1;
for(int i=2;i<=ty;++i)if(Sy[i]-Sy[i-1]>1)Sy[++ty]=Sy[i]-1;
Sx[++tx]=-10;Sy[++ty]=-10;Sx[++tx]=1e9;Sy[++ty]=1e9;
sort(&Sx[1],&Sx[tx+1]);sort(&Sy[1],&Sy[ty+1]);
tx=unique(&Sx[1],&Sx[tx+1])-Sx-1;ty=unique(&Sy[1],&Sy[ty+1])-Sy-1;
for(int i=1;i<=n;++i)
{
p[i].x1=lower_bound(&Sx[1],&Sx[tx+1],p[i].x1)-Sx;
p[i].x2=lower_bound(&Sx[1],&Sx[tx+1],p[i].x2)-Sx;
p[i].y1=lower_bound(&Sy[1],&Sy[ty+1],p[i].y1)-Sy;
p[i].y2=lower_bound(&Sy[1],&Sy[ty+1],p[i].y2)-Sy;
for(int j=p[i].x1;j<=p[i].x2;++j)visx[j][p[i].y1-1]=visx[j][p[i].y2-1]=true;
for(int j=p[i].y2;j<=p[i].y1;++j)visy[p[i].x1][j]=visy[p[i].x2][j]=true;
}
Bx=lower_bound(&Sx[1],&Sx[tx+1],Bx)-Sx;By=lower_bound(&Sy[1],&Sy[ty+1],By)-Sy;
Ex=lower_bound(&Sx[1],&Sx[tx+1],Ex)-Sx;Ey=lower_bound(&Sy[1],&Sy[ty+1],Ey)-Sy;
for(int i=1;i<=tx;++i)
for(int j=1;j<=ty;++j)
for(int k=0;k<4;++k)
{
int x=i+d[k][0],y=j+d[k][1];
if(x<1||y<1||x>tx||y>ty)continue;
if(k==0)Add(id(i,j),id(x,y),visy[x][y]);
if(k==1)Add(id(i,j),id(x,y),visy[i][j]);
if(k==2)Add(id(i,j),id(x,y),visx[i][j]);
if(k==3)Add(id(i,j),id(x,y),visx[x][y]);
}
SPFA();
printf("%d
",dis[id(Ex,Ey)]);
return 0;
}