【网络流24题23】火星探险问题

题面戳我
不想放题面了。。。

sol

首先还是费用流建模,建模方式参见【网络流24题20】深海机器人问题。这题的题解我还没有写。。。
针对输出方案:在残余网络上搞一搞可以求得每个位置\((x,y)\)被经过的次数。然后做k(探测车的数量)次,每次从\((1,1)\)跑到\((n,m)\),把经过的点上面的次数全部减一。这样走肯定是对的。

code

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
const int inf = 1e9;
const int N = 40;
struct edge{int to,next,cap,w,cost;}a[N*N<<5];
int k,n,m,num,s,t,mmp[N][N],P[N][N],head[2*N*N],cnt=1,dis[2*N*N],vis[2*N*N],pe[2*N*N],tot[N][N];
queue<int>Q;
int gi()
{
    int x=0,w=1;char ch=getchar();
    while ((ch<'0'||ch>'9')&&ch!='-') ch=getchar();
    if (ch=='-') w=0,ch=getchar();
    while (ch>='0'&&ch<='9') x=(x<<3)+(x<<1)+ch-'0',ch=getchar();
    return w?x:-x;
}
void link(int u,int v,int w,int cost)
{
    a[++cnt]=(edge){v,head[u],w,w,cost};
    head[u]=cnt;
    a[++cnt]=(edge){u,head[v],0,0,-cost};
    head[v]=cnt;
}
bool spfa()
{
    memset(dis,63,sizeof(dis));
    dis[s]=0;Q.push(s);
    while (!Q.empty())
    {
        int u=Q.front();Q.pop();
        for (int e=head[u];e;e=a[e].next)
        {
            int v=a[e].to;
            if (a[e].w&&dis[v]>dis[u]+a[e].cost)
            {
                dis[v]=dis[u]+a[e].cost;pe[v]=e;
                if (!vis[v]) vis[v]=1,Q.push(v);
            }
        }
        vis[u]=0;
    }
    if (dis[t]==dis[0]) return false;
    int sum=inf;
    for (int i=t;i!=s;i=a[pe[i]^1].to)
        sum=min(sum,a[pe[i]].w);
    for (int i=t;i!=s;i=a[pe[i]^1].to)
        a[pe[i]].w-=sum,a[pe[i]^1].w+=sum;
    return true;
}
int main()
{
    k=gi();m=gi();n=gi();
    for (int i=1;i<=n;i++)
        for (int j=1;j<=m;j++)
        {
            mmp[i][j]=gi();
            P[i][j]=++num;
        }
    for (int i=1;i<=n;i++)
        for (int j=1;j<=m;j++)
        {
            if (mmp[i][j]==2)
                link(P[i][j],P[i][j]+num,1,-1);
            if (mmp[i][j]!=1)
                link(P[i][j],P[i][j]+num,inf,0);
            if (i<n)
                link(P[i][j]+num,P[i+1][j],inf,0);
            if (j<m)
                link(P[i][j]+num,P[i][j+1],inf,0);
        }
    s=num*2+1;link(s,P[1][1],k,0);
    t=s+1;link(P[n][m]+num,t,k,0);
    while (spfa());
    for (int i=1;i<=n;i++)
        for (int j=1;j<=m;j++)
            for (int e=head[P[i][j]];e;e=a[e].next)
                if (a[e].to==P[i][j]+num)
                    tot[i][j]+=a[e].cap-a[e].w;
    for (int t=1;t<=k;t++)
    {
        int x=1,y=1;
        while (x!=n||y!=m)
        {
            if (tot[x][y+1])
                ++y,--tot[x][y],printf("%d 1\n",t);
            else
                ++x,--tot[x][y],printf("%d 0\n",t);
        }
    }
    return 0;
}

相关推荐