【NOIP2016提高A组集训第7场11.4】推冰块

Description

Dpstr最近迷上了推冰块。冰地是一个n行m列的网格区域，第i行第j列的格子记为(i,j)，也就是左上角为(1,1)，右下角为(n,m)。每个格子可能是冰面、障碍物、减速带三者之一。其中，冰地外围（即第0行、第n+1行、第0列、第m+1列）的所有格子均有障碍物。除此之外，冰地内共有k个障碍物和减速带，其余格子为冰面。
初始时，有一个冰块位于(1,1)处。Dpstr每次可以选择上、下、左、右四个方向之一推动该冰块，推动后该冰块将一直沿此方向移动，直到冰块所在的格子为减速带，或冰块沿运动方向的下一个格子为障碍物时，冰块将停止运动。一旦冰块停在减速带上，该减速带即消失。
如下图，当n = 5，m = 5时，若冰块位于(3,1)，且(3,4)处有一个障碍物，那么向右推动冰块（图1中A处），冰块将停在(3,3)（图1中B处）；而如果(3,4)处不是障碍物而是减速带，那么冰块将停在(3,4)上（图2中B处），之后(3,4)处的减速带就会消失。

Solution

一看就是暴搜，一个减速带只能经过一次，这个很显然。
一开始想spfa，结果一直没有打出来。
才发现spfa是多么的麻烦。
直接bfs，然后队首的点向4周二分的找出所有点就好了。

Code

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define fo(i,a,b) for(i=a;i<=b;i++)
#define fod(i,a,b) for(i=a;i>=b;i--)
using namespace std;
typedef long long ll;
const int maxn=60007,mo=6000007,mxx=1007;
int i,j,k,l,r,t,n,m,ans,u,v,z,mid;
int num,head,tail,x,y,er;
struct node{int x,y,t;
}a[maxn],b[maxn],c[maxn];
int fang[4][2]={1,0,0,1,-1,0,0,-1};
int fan[5]={0,3,4,1,2};
bool cmp(node x,node y){return x.x<y.x||x.x==y.x&&x.y<y.y;}
bool cmp1(node x,node y){return x.y<y.y||x.y==y.y&&x.x<y.x;}
int data[maxn*100][3];
ll h[mo];
int d[4],e[4],f[4];
bool bz[4];
ll suan(int x,int y){ll z=(x-1)*m+y;return z;
}
int hash(int x,int y){ll z=suan(x,y);int o=z%mo;while(h[o]!=z&&h[o]!=0)o=(o+1)%mo;return o;
}
void jia(int a,int b,int c){if(!a)return;if(a==n&&b==m){ans=min(ans,c+1);return;}er=hash(a,b);if(h[er])return;h[er]=suan(a,b);data[++tail][0]=a,data[tail][1]=b,data[tail][2]=c+1;
}
int main(){freopen("ice.in","r",stdin);freopen("ice.out","w",stdout);//  freopen("fan.in","r",stdin);scanf("%d%d%d",&n,&m,&k);fo(i,1,k){scanf("%d%d%d",&a[i].x,&a[i].y,&a[i].t);b[i]=a[i],c[i]=a[i];}sort(b+1,b+1+k,cmp);sort(c+1,c+1+k,cmp1);ans=0x7fffffff;jia(1,1,-1);while(head<tail){x=data[++head][0],y=data[head][1];memset(d,0,sizeof(d));memset(e,0,sizeof(e));l=1,r=k;while(l<r){mid=(l+r+1)/2;if(b[mid].x>x||b[mid].x==x&&b[mid].y>y)r=mid-1;else l=mid;}fo(i,max(1,l-2),l+2){if(b[i].x==x&&b[i].y<y){d[3]=b[i].x,e[3]=b[i].y+(b[i].t==0);  }if(b[i].x==x&&b[i].y>y){d[1]=b[i].x,e[1]=b[i].y-(b[i].t==0);break;}if(i==k)break;}if(!d[1])d[1]=x,e[1]=m;if(!d[3])d[3]=x,e[3]=1;l=1,r=k;while(l<r){mid=(l+r+1)/2;if(c[mid].y>y||c[mid].y==y&&c[mid].x>x)r=mid-1;else l=mid;}fo(i,max(1,l-2),l+2){if(c[i].y==y&&c[i].x<x){d[2]=c[i].x+(c[i].t==0),e[2]=c[i].y;  }if(c[i].y==y&&c[i].x>x){d[0]=c[i].x-(c[i].t==0),e[0]=c[i].y;break;}if(i==k)break;}if(!d[0])d[0]=n,e[0]=y;if(!d[2])d[2]=1,e[2]=y;fo(i,0,3)jia(d[i],e[i],data[head][2]);}printf("%d\n",ans);
}