1. 程式人生 > >NOIP2013 華容道 (棋盤建圖+spfa最短路)

NOIP2013 華容道 (棋盤建圖+spfa最短路)

sizeof bool 定位 wid 開始 build ron 內容 com

  1 #include <cstdio>
  2 #include <algorithm>
  3 #include <cstring>
  4 #include <queue>
  5 #define inf 0x3f3f3f3f
  6 #define N 35
  7 #define maxn 5000
  8 #define mod 1000000007
  9 #define ll long long 
 10 using namespace std;
 11 
 12 int n,m,q,cte;
 13 int ok[N][N],id[N][N],head[maxn],dis[maxn],inq[maxn],st[4
]; 14 int xx[]={-1,0,1,0}; 15 int yy[]={0,1,0,-1}; 16 struct EDGE{ 17 int to,nxt,val; 18 }edge[maxn*500]; 19 struct node{ 20 int x,y; 21 }; 22 node ins(int s1,int s2) {node k;k.x=s1;k.y=s2;return k;} 23 bool check(int x,int y) 24 { 25 if(x<1||y<1||x>n||y>m||!ok[x][y]) return
false; 26 else return true; 27 } 28 void edge_add(int u,int v,int w) 29 { 30 cte++; 31 edge[cte].to = v; 32 edge[cte].nxt=head[u]; 33 edge[cte].val= w; 34 head[u]=cte; 35 } 36 void bfs(int ex,int ey,int sx,int sy,int d) 37 { 38 queue<node>q; 39 memset(dis,-1
,sizeof(dis)); 40 dis[id[ex][ey]]=0,dis[id[sx][sy]]=1; 41 q.push(ins(ex,ey)); 42 int x,y,fx,fy; 43 while(!q.empty()) 44 { 45 node k=q.front();q.pop(); 46 x=k.x,y=k.y; 47 for(int i=0;i<4;i++) 48 { 49 fx=x+xx[i],fy=y+yy[i]; 50 if(!check(fx,fy)||dis[id[fx][fy]]!=-1) continue; 51 dis[id[fx][fy]]=dis[id[x][y]]+1; 52 q.push(ins(fx,fy)); 53 } 54 } 55 for(int i=0;i<4;i++) 56 { 57 fx=sx+xx[i],fy=sy+yy[i]; 58 if(!check(fx,fy)||dis[id[fx][fy]]==-1||i==d) continue; 59 edge_add(id[sx][sy]+d,id[sx][sy]+i,dis[id[fx][fy]]); 60 } 61 edge_add(id[sx][sy]+d,id[ex][ey]+(d+2)%4,1); 62 } 63 void bfsS(int ex,int ey,int sx,int sy) 64 { 65 queue<node>q; 66 memset(dis,-1,sizeof(dis)); 67 q.push(ins(ex,ey)),dis[id[ex][ey]]=0; 68 int x,y,fx,fy; 69 while(!q.empty()) 70 { 71 node k=q.front();q.pop(); 72 x=k.x,y=k.y; 73 for(int i=0;i<4;i++) 74 { 75 fx=x+xx[i],fy=y+yy[i]; 76 if(!check(fx,fy)||dis[id[fx][fy]]!=-1||(fx==sx&&fy==sy)) continue; 77 dis[id[fx][fy]]=dis[id[x][y]]+1; 78 q.push(ins(fx,fy)); 79 } 80 } 81 } 82 void build_edge() 83 { 84 int ct=1; 85 for(int i=1;i<=n;i++) 86 for(int j=1;j<=m;j++) 87 { 88 if(ok[i][j]) id[i][j]=ct,ct+=4; 89 } 90 for(int i=1;i<=n;i++) 91 for(int j=1;j<=m;j++) 92 { 93 if(!check(i,j)) continue; 94 if(check(i-1,j)) bfs(i-1,j,i,j,0); 95 if(check(i,j+1)) bfs(i,j+1,i,j,1); 96 if(check(i+1,j)) bfs(i+1,j,i,j,2); 97 if(check(i,j-1)) bfs(i,j-1,i,j,3); 98 } 99 } 100 void spfa(int S) 101 { 102 queue<int>q; 103 memset(dis,0x3f,sizeof(dis)); 104 q.push(S),dis[S]=0,inq[S]=1; 105 while(!q.empty()) 106 { 107 int u=q.front();q.pop(); 108 for(int j=head[u];j!=-1;j=edge[j].nxt) 109 { 110 int v=edge[j].to; 111 if(dis[v]>dis[u]+edge[j].val) 112 { 113 dis[v]=dis[u]+edge[j].val; 114 if(!inq[v]) q.push(v),inq[v]=1; 115 } 116 } 117 inq[u]=0; 118 } 119 } 120 121 int solve(int ex,int ey,int sx,int sy,int tx,int ty) 122 { 123 bfsS(ex,ey,sx,sy); 124 if(sx==tx&&sy==ty) return 0; 125 for(int i=0;i<4;i++) st[i]=dis[id[sx+xx[i]][sy+yy[i]]]; 126 int fx,fy,ans=inf; 127 for(int i=0;i<4;i++) 128 { 129 fx=sx+xx[i],fy=sy+yy[i]; 130 if(!check(fx,fy)||st[i]==-1) continue; 131 spfa(id[sx][sy]+i); 132 for(int j=0;j<4;j++) 133 { 134 if(!check(tx+xx[j],ty+yy[j])||dis[id[tx][ty]+j]==-1) continue; 135 ans=min(ans,st[i]+dis[id[tx][ty]+j]); 136 } 137 } 138 if(ans>=inf) return -1; 139 else return ans; 140 } 141 142 int main() 143 { 144 freopen("aa.in","r",stdin); 145 scanf("%d%d%d",&n,&m,&q); 146 for(int i=1;i<=n;i++) 147 for(int j=1;j<=m;j++) 148 scanf("%d",&ok[i][j]); 149 memset(head,-1,sizeof(head)); 150 build_edge(); 151 int ex,ey,sx,sy,tx,ty; 152 while(q--) 153 { 154 scanf("%d%d%d%d%d%d",&ex,&ey,&sx,&sy,&tx,&ty); 155 printf("%d\n",solve(ex,ey,sx,sy,tx,ty)); 156 } 157 return 0; 158 }

---恢復內容開始---

題目大意:類似於華容道遊戲,但所有棋子都是1技術分享圖片1大小的,棋盤上僅有一個空格,其它位置填滿了棋子,有些棋子是固定的不能被移動,每個格子內只能有一個棋子,棋子只能移動到相鄰的空格上,每移動任意棋子一次視為一次操作,求讓指定棋子移動到指定位置的最小操作數

據說爆搜能拿70pts....

觀察爆搜為什麽會T,主要是因為搜出了許多無用的狀態

考慮優化這個過程,空格在指定棋子的上下左右分別記為一種狀態,相當於一個位置記錄了4個狀態

為了方便轉移,直接把每個格子的四個狀態記成連續的,這樣方便轉移

每種狀態,它能轉移到其它狀態的情況分為兩種:

<1>指定格子的某個方向有空格,轉移到同一個格子的另一個方向有空格,指定格子上的棋子的位置不動

這種情況,我們對於每個格子的四個方向分別bfs即可,因為除了空格子,其它位置都被填滿了,所以空格子在被填滿的棋盤裏移動,和棋子在空棋盤裏的移動方式是一樣的(貌似是一句廢話但我仍然思考了很久...)

<2>指定棋子朝著這個狀態指向的空格移動

這種情況,顯然操作次數為1,轉移的狀態呢,則是從 當前格子的指向 轉移到 目標格子的反方向 (就是那個(d+2)%4)

雖然有500次詢問,但圖都是一樣的,所以建一次圖就行

那麽,如果一開始空格子不在指定棋子邊上呢?

我們把空格子移動到指定棋子的四個方向就行了,註意不要讓空格子移動到指定棋子的格子上

再分別跑一次SPFA最短路即可

總結:之前還有一道棋盤建圖題(荷葉塘),但那道題建圖更簡單更好想,思路都是先bfs建邊,然後跑最短路

註意特判移動次數為0的情況!

#include <cstdio>
#include <algorithm>
#include <cstring>
#include <queue>
#define inf 0x3f3f3f3f
#define N 35
#define maxn 5000
#define mod 1000000007
#define ll long long 
using namespace std;

int n,m,q,cte;
int ok[N][N],id[N][N],head[maxn],dis[maxn],inq[maxn],st[4];
int xx[]={-1,0,1,0};
int yy[]={0,1,0,-1};
struct EDGE{
    int to,nxt,val;
}edge[maxn*500];
struct node{
    int x,y;
};
node ins(int s1,int s2) {node k;k.x=s1;k.y=s2;return k;}
bool check(int x,int y)
{
    if(x<1||y<1||x>n||y>m||!ok[x][y]) return false;
    else return true;
}
void edge_add(int u,int v,int w)
{
    cte++;
    edge[cte].to = v;
    edge[cte].nxt=head[u];
    edge[cte].val= w;
    head[u]=cte;
}
void bfs(int ex,int ey,int sx,int sy,int d)
{
    queue<node>q;
    memset(dis,-1,sizeof(dis));
    dis[id[ex][ey]]=0,dis[id[sx][sy]]=1;
    q.push(ins(ex,ey));
    int x,y,fx,fy;
    while(!q.empty())
    {
        node k=q.front();q.pop();
        x=k.x,y=k.y;
        for(int i=0;i<4;i++)
        {
            fx=x+xx[i],fy=y+yy[i];
            if(!check(fx,fy)||dis[id[fx][fy]]!=-1) continue;
            dis[id[fx][fy]]=dis[id[x][y]]+1;
            q.push(ins(fx,fy));
        }
    }
    for(int i=0;i<4;i++)
    {
        fx=sx+xx[i],fy=sy+yy[i];
        if(!check(fx,fy)||dis[id[fx][fy]]==-1||i==d) continue;
        edge_add(id[sx][sy]+d,id[sx][sy]+i,dis[id[fx][fy]]);
    }    
    edge_add(id[sx][sy]+d,id[ex][ey]+(d+2)%4,1);
}
void bfsS(int ex,int ey,int sx,int sy)
{
    queue<node>q;
    memset(dis,-1,sizeof(dis));
    q.push(ins(ex,ey)),dis[id[ex][ey]]=0;
    int x,y,fx,fy;
    while(!q.empty())
    {
        node k=q.front();q.pop();
        x=k.x,y=k.y;
        for(int i=0;i<4;i++)
        {
            fx=x+xx[i],fy=y+yy[i];
            if(!check(fx,fy)||dis[id[fx][fy]]!=-1||(fx==sx&&fy==sy)) continue;
            dis[id[fx][fy]]=dis[id[x][y]]+1;
            q.push(ins(fx,fy));
        }
    }
}
void build_edge()
{
    int ct=1;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
        {
            if(ok[i][j]) id[i][j]=ct,ct+=4;
        } 
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
        {
            if(!check(i,j)) continue;
            if(check(i-1,j)) bfs(i-1,j,i,j,0);
            if(check(i,j+1)) bfs(i,j+1,i,j,1);
            if(check(i+1,j)) bfs(i+1,j,i,j,2);
            if(check(i,j-1)) bfs(i,j-1,i,j,3);
        } 
}
void spfa(int S)
{
    queue<int>q;
    memset(dis,0x3f,sizeof(dis));
    q.push(S),dis[S]=0,inq[S]=1;
    while(!q.empty())
    {
        int u=q.front();q.pop();
        for(int j=head[u];j!=-1;j=edge[j].nxt)
        {
            int v=edge[j].to;
            if(dis[v]>dis[u]+edge[j].val)
            {
                dis[v]=dis[u]+edge[j].val;
                if(!inq[v]) q.push(v),inq[v]=1;
            }
        }
        inq[u]=0;
    }
}

int solve(int ex,int ey,int sx,int sy,int tx,int ty)
{
    bfsS(ex,ey,sx,sy);
    if(sx==tx&&sy==ty) return 0;
    for(int i=0;i<4;i++) st[i]=dis[id[sx+xx[i]][sy+yy[i]]];
    int fx,fy,ans=inf;
    for(int i=0;i<4;i++)
    {
        fx=sx+xx[i],fy=sy+yy[i];
        if(!check(fx,fy)||st[i]==-1) continue;
        spfa(id[sx][sy]+i);
        for(int j=0;j<4;j++)
        {
            if(!check(tx+xx[j],ty+yy[j])||dis[id[tx][ty]+j]==-1) continue;
            ans=min(ans,st[i]+dis[id[tx][ty]+j]);
        } 
    }
    if(ans>=inf) return -1;
    else return ans;
}

int main()
{
    freopen("aa.in","r",stdin);
    scanf("%d%d%d",&n,&m,&q);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            scanf("%d",&ok[i][j]);
    memset(head,-1,sizeof(head));
    build_edge();
    int ex,ey,sx,sy,tx,ty;
    while(q--)
    {
        scanf("%d%d%d%d%d%d",&ex,&ey,&sx,&sy,&tx,&ty);
        printf("%d\n",solve(ex,ey,sx,sy,tx,ty));
    }
    return 0;
}

---恢復內容結束---

NOIP2013 華容道 (棋盤建圖+spfa最短路)