[GDOI2015]推箱子(状压bfs)

[GDOI2015]推箱子(状压bfs)

题面

题面过长,略

分析

观察到$m \times m =64 \(,那么可以把箱子的01状态压到一个```unsigned long long```里面 然后对于地图上的每一个点\)(x,y)\(,预处理出左上角在\)(x,y)\(,边长为\)m\(的正方形的01状态。如果这个状态和箱子的状态按位与的结果为0,那么就说明箱子可以通过. 然后发现这类似一个分层图上的最短路问题,直接BFS即可.状态\)dist[x][y][k]\(表示箱子左上角在\)(x,y)\(,旋转了\)k \times 90 \degree$.分题目中的三种情况转移即可.

代码

#include<iostream>
#include<cstdio>
#include<cstring> 
#include<queue>
#define maxm 8
#define INF 0x3f3f3f3f
#define maxn 2000
using namespace std;
typedef unsigned long long ull;
int n,m;
int ma[maxn+5][maxn+5],box[maxm+5][maxm+5];
inline void rotate90(){
    static int tmp[maxm+5][maxm+5];
    for(int i=0;i<m;i++){
        for(int j=0;j<m;j++) tmp[m-j-1][i]=box[i][j];
    }
    for(int i=0;i<m;i++){
        for(int j=0;j<m;j++) box[i][j]=tmp[i][j];
    }
} 
ull box_bin(){
    ull ans=0;
    for(int i=0;i<m;i++){
        for(int j=0;j<m;j++) ans=ans<<1|box[i][j];
    }
    return ans;
}
ull ma_bin(int x,int y){
    ull ans=0;
    for(int i=0;i<m;i++){
        for(int j=0;j<m;j++) ans=ans<<1|ma[x+i][y+j]; 
    }
    return ans;
}

ull f[4];
ull g[maxn+5][maxn+5];

const int dirx[4]={1,-1,0,0},diry[4]={0,0,1,-1};
struct node{
    int x;
    int y;//左上角坐标 
    int k;//旋转到哪种状态 
    node(){
        
    }
    node(int _x,int _y,int _k){
        x=_x;
        y=_y;
        k=_k;
    }
};
int dist[maxn+5][maxn+5][4];
void bfs(){
    queue<node>q;
    memset(dist,0x3f,sizeof(dist));
    dist[0][0][0]=0;
    q.push(node(0,0,0));
    while(!q.empty()){
        node p=q.front();
        q.pop();
        int x=p.x,y=p.y,k=p.k;
        int xx,yy,kk; 
        kk=(k+1)&3;//顺时针转
        if(((f[kk]&g[x][y])==0)&&dist[x][y][k]+1<dist[x][y][kk]){
            q.push(node(x,y,kk));
            dist[x][y][kk]=dist[x][y][k]+1;
        }
        kk=(k-1)&3;//逆时针转
        if(((f[kk]&g[x][y])==0)&&dist[x][y][k]+1<dist[x][y][kk]){
            q.push(node(x,y,kk));
            dist[x][y][kk]=dist[x][y][k]+1;
        }
        for(int i=0;i<4;i++){//走 
            xx=x+dirx[i],yy=y+diry[i];
            if(xx>=0&&yy>=0&&xx<=n-m&&yy<=n-m){
                if(((f[k]&g[xx][yy])==0)&&dist[x][y][k]+1<dist[xx][yy][k]){
                    q.push(node(xx,yy,k));
                    dist[xx][yy][k]=dist[x][y][k]+1;
                }
            }
        }
    }
}
int main(){
    scanf("%d %d",&n,&m);
    for(int i=0;i<m;i++){
        for(int j=0;j<m;j++) scanf("%1d",&box[i][j]);
    }
    for(int i=0;i<n;i++){
        for(int j=0;j<n;j++) scanf("%1d",&ma[i][j]);
    }
    for(int i=0;i<4;i++){
        f[i]=box_bin();
        rotate90(); 
    } 
    for(int i=0;i<n-m+1;i++){
        for(int j=0;j<n-m+1;j++) g[i][j]=ma_bin(i,j); 
    } 
    bfs();
    int ans=INF;
    for(int i=0;i<4;i++) ans=min(ans,dist[n-m][n-m][i]);
    if(ans==INF) printf("-1\n");
    else printf("%d\n",ans);
}