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

题目描述

Dpstr最近迷上了推冰块。冰地是一个n行m列的网格区域，第i行第j列的格子记为(i,j)，也就是左上角为(1,1)，右下角为(n,m)。每个格子可能是冰面、障碍物、减速带三者之一。其中，冰地外围（即第0行、第n+1行、第0列、第m+1列）的所有格子均有障碍物。除此之外，冰地内共有k个障碍物和减速带，其余格子为冰面。
初始时，有一个冰块位于(1,1)处。Dpstr每次可以选择上、下、左、右四个方向之一推动该冰块，推动后该冰块将一直沿此方向移动，直到冰块所在的格子为减速带，或冰块沿运动方向的下一个格子为障碍物时，冰块将停止运动。一旦冰块停在减速带上，该减速带即消失。
Dpstr希望通过尽量少的推动次数使得冰块停在(n,m)处。请计算Dpstr至少要推多少次冰块。

数据范围

对于30%的数据，2≤n≤5，2≤m≤5，0≤k≤5；
对于50%的数据，2≤n≤1,000，2≤m≤1,000，0≤k≤1,000；
对于70%的数据，2≤n≤50,000，2≤m≤50,000，0≤k≤50,000；
对于100%的数据，2≤n≤1,000,000,000，2≤m≤1,000,000,000，0≤k≤50,000，1≤xi≤n，1≤yi≤m，0≤ti≤1。

解法

显然每个格子最多只会走一次。
所以可以使用BFS。
运用二分解决滑行问题。

代码

#include<iostream>
#include<stdio.h>
#include<string.h>
#include<math.h>
#include<algorithm>
#define ll long long
using namespace std;
const char* fin="ice.in";
const char* fout="ice.out";
const ll inf=0x7fffffff;
const ll maxn=50007,maxt=1000007;
const ll fx[4][2]={{0,-1},{-1,0},{1,0},{0,1}};
ll n,m,n1,i,j,k,l,ans=inf;
struct node{
    ll x,y,z;
    void operator =(const node &j){
        x=j.x;
        y=j.y;
        z=j.z;
    }
}f[maxn],g[maxn];
ll b[maxt][2],head,tail;
ll h[maxt],dis[maxt];
ll hash(ll x){
    ll t=x%maxt;
    while (h[t] && h[t]!=x) t=(t+1)%maxt;
    return t;
}
ll pos(ll x,ll y){
    return (x-1)*m+y;
}
void add(ll x,ll y,ll z){
    ll k=pos(x,y),i,j;
    i=hash(k);
    if (!h[i]) {
        h[i]=k;
        dis[i]=z;
        b[++tail][0]=x;
        b[tail][1]=y;
    }
}
void bfs(ll x,ll y){
    ll xx[4],yy[4],i,j,k,dist,now,l,r,mid,o;
    bool bz[4];
    add(x,y,0);
    while (head++<tail){
        j=b[head][0];
        k=b[head][1];
        now=hash(pos(j,k));
        memset(bz,0,sizeof(bz));
        l=1;
        r=n1;
        while (l<r){
            mid=(l+r+1)/2;
            if (f[mid].x<j || f[mid].x==j && f[mid].y<k) l=mid;
            else r=mid-1;
        }
        for (o=l;o<=l+2;o++){
            if (o>n1) break;
            if (f[o].x==j && f[o].y<k) {
                xx[0]=f[o].x;
                yy[0]=f[o].y+(f[o].z==0);
                bz[0]=true;
            }
            if (f[o].x==j && f[o].y>k){
                xx[1]=f[o].x;
                yy[1]=f[o].y-(f[o].z==0);
                bz[1]=true;
                break;
            }
        }
        if (!bz[0]) xx[0]=j,yy[0]=1;
        if (!bz[1]) xx[1]=j,yy[1]=m;
        l=1;
        r=n1;
        while (l<r){
            mid=(l+r+1)/2;
            if (g[mid].y<k || g[mid].y==k && g[mid].x<j) l=mid;
            else r=mid-1;
        }
        for (o=l;o<=l+2;o++){
            if (o>n1) break;
            if (g[o].y==k && g[o].x<j) {
                xx[2]=g[o].x+(g[o].z==0);
                yy[2]=g[o].y;
                bz[2]=true;
            }
            if (g[o].y==k && g[o].x>j){
                xx[3]=g[o].x-(g[o].z==0);
                yy[3]=g[o].y;
                bz[3]=true;
                break;
            }
        }
        if (!bz[2]) xx[2]=1,yy[2]=k;
        if (!bz[3]) xx[3]=n,yy[3]=k;
        for (i=0;i<4;i++){
            add(xx[i],yy[i],dis[now]+1);
            if (xx[i]==n && m==yy[i]) ans=min(ans,dis[now]+1);
        }
    }
}
bool cmp(node a,node b){
    return a.x<b.x || a.x==b.x && a.y<b.y;
}
bool cmp1(node a,node b){
    return a.y<b.y || a.y==b.y && a.x<b.x;
}
int main(){
    freopen(fin,"r",stdin);
    freopen(fout,"w",stdout);
    scanf("%d%d%d",&n,&m,&n1);
    for (i=1;i<=n1;i++) scanf("%d%d%d",&f[i].x,&f[i].y,&f[i].z),g[i]=f[i];
    sort(f+1,f+n1+1,cmp);
    sort(g+1,g+n1+1,cmp1);
    bfs(1,1);
    printf("%lld",ans);
    return 0;
}

启发

对于搜索题，从最简单的搜索入手。
初始搜索最好打bfs，因为最容易优化，而且最短路最快找到。
然后根据客观的时间耗费，优化对应的地方。