[noip2013]华容道 题解

本文介绍了一道结合广度优先搜索(BFS)与最短路径快速算法(SPFA)的编程题。通过预处理地图数据,利用BFS解决部分问题,并通过SPFA求解最短路径。代码实现细节丰富,包括地图读取、状态转移等。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

才说完都很短,然后就一道码农题233。

首先70分的BFS是很好想到的,然后就可以愉快地开始码这道题了。

但是,这是远远不够的,最后的点如过呢。

我们可以预处理一下,然后spfa,从每个状态到它可以转移的状态就好了。

最后,这种题真心恶心呀。

#include<bits/stdc++.h>
#define N 40000
#define INF 0x3f3f3f3f
using namespace std;
int dx[4]={1,-1,0,0},dy[4]={0,0,1,-1};
int first[N+5],nxt[N+5],to[N+5],val[N+5],siz,tot,dis[N+5];
int a[35][35],num[35][35][5];
int ex,ey,sx,sy,tx,ty,n,m,Q;
bool vis[35][35],f[N+5];
struct H{
    int x,y,step;
}; 
inline char nc()
{
    static char buf[100000],*p1=buf,*p2=buf;
    return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}
inline int read()
{
    int x=0,b=1;
    char c=nc();
    for(;!(c<='9'&&c>='0');c=nc())if(c=='-')b=-1;
    for(;c<='9'&&c>='0';c=nc())x=x*10+c-'0';
    return x*b;
}
void add(int x,int y,int z)
{
    to[siz]=y;
    nxt[siz]=first[x];
    first[x]=siz;
    val[siz]=z;
    siz++;
}
int bfs(int ax,int ay,int bx,int by,int cx,int cy)
{
    if(ax==bx&&ay==by) return 0;
    memset(vis,0,sizeof(vis));
    queue <H> q;
    q.push((H){ax,ay,0});
    vis[ax][ay]=1;
    while(!q.empty())
    {
        H k=q.front();q.pop();
        if(k.x==bx&&k.y==by) return k.step;
        for(int i=0;i<4;i++)
        {
            int x=k.x+dx[i],y=k.y+dy[i];
            if(x>=1&&x<=n&&y>=1&&y<=m&&a[x][y]&&(x!=cx||y!=cy)&&!vis[x][y])
            {
                q.push((H){x,y,k.step+1});
                vis[x][y]=1;
                if(x==bx&&y==by) return k.step+1;
            }
        }
    }
    return INF;
} 
int spfa()
{
    queue <int> q;
    if(sx==tx&&sy==ty) return 0;
    for(int i=1;i<=tot;i++) dis[i]=INF;
    for(int k=0;k<4;k++)
     if(num[sx][sy][k]){
        dis[num[sx][sy][k]]=bfs(ex,ey,sx+dx[k],sy+dy[k],sx,sy);
        q.push(num[sx][sy][k]); 
    }
    while(!q.empty())
    {
        int k=q.front();q.pop();f[k]=0;
        for(int i=first[k];i!=-1;i=nxt[i])
        {
            int t=to[i];
            if(dis[t]<=dis[k]+val[i])continue;
            dis[t]=dis[k]+val[i];
            if(!f[t]) q.push(t),f[t]=1;
        }
    }
    int ans=INF;
    for(int k=0;k<4;k++)
     if (num[tx][ty][k])
     ans=min(ans,dis[num[tx][ty][k]]);
    return ans==INF?-1:ans;
}
void init()
{
    memset(first,-1,sizeof(first));
    for(int i=1;i<=n;i++)
     for(int j=1;j<=m;j++)
      for(int k=0;k<4;k++)
       if(a[i][j]&&a[i+dx[k]][j+dy[k]])num[i][j][k]=++tot;
    for(int i=1;i<=n;i++)
     for(int j=1;j<=m;j++)
      for(int k=0;k<4;k++)
       if(num[i][j][k])
        add(num[i][j][k],num[i+dx[k]][j+dy[k]][k^1],1);
    for(int i=1;i<=n;i++)
     for(int j=1;j<=m;j++)
      for(int k=0;k<4;k++)
       for(int l=0;l<4;l++)
        if(l!=k&&num[i][j][k]&&num[i][j][l])
         add(num[i][j][k],num[i][j][l],bfs(i+dx[k],j+dy[k],i+dx[l],j+dy[l],i,j));    
}
int main()
{
    freopen("in.txt","r",stdin);
    n=read(),m=read(),Q=read();
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            a[i][j]=read();
    init();
    while(Q--)
    {
        ex=read(),ey=read(),sx=read(),sy=read(),tx=read(),ty=read();
        printf("%d\n",spfa());
    }
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值