火星探险问题

这道题不就是一个拆点后的少连边输出方案的方格取数嘛。方案可以用栈来输出。
有一点很妙的地方是，反向边走过以后就有值了，可以这样来打印方案。

为什么我打印方案那么多注释？周五做题真的有毒。。。我方案打印不对。。改了一年，距离写着题6个小时后，我发现。。我SPFA打错了。。没有vis等于0。永远不要确信你某个地方不会写错就不去检查它，即使一个算法你打的再熟练，理解再深刻，你也有手滑打错的可能。而这种可能，你忽视它的话，可能你整个考试就炸了。

#include<bits/stdc++.h>using namespace std;const int MAXN=7e5+5;const int base=3e5;const int INF=1e9+7;struct edge{    int u,to,next,w,c;}e[MAXN<<1];int car,p,q,s,t;int head[MAXN],cnt=1;inline void add(int u,int v,int w,int cost){    e[++cnt]=(edge){u,v,head[u],w,cost},head[u]=cnt;    e[++cnt]=(edge){v,u,head[v],0,-cost},head[v]=cnt;}int ma[205][205];queue<int>qq;int dis[MAXN],pre[MAXN];bool vis[MAXN];bool SPFA(int x){    memset(pre,0,sizeof(pre));    memset(vis,0,sizeof(vis));    while(qq.size())qq.pop();    for(int i=s;i<=t;i++)dis[i]=-INF;    qq.push(x);    dis[x]=0;vis[x]=1;    while(qq.size()){        int u=qq.front();qq.pop();vis[u]=0;        for(int i=head[u];i;i=e[i].next){            int v=e[i].to,w=e[i].c;            if(e[i].w){                if(dis[u]+w>dis[v]){                    dis[v]=dis[u]+w;pre[v]=i;                    if(!vis[v]){                        qq.push(v);vis[v]=1;                    }                }            }        }    }    if(dis[t]==-INF)return 0;    return 1;}int MC(){    int ans=0;    while(SPFA(s)){        int tem=INF;        for(int i=pre[t];i;i=pre[e[i].u]){            tem=min(tem,e[i].w);        }        for(int i=pre[t];i;i=pre[e[i].u]){            ans+=e[i].c*tem;            e[i].w-=tem;e[i^1].w+=tem;        }    }    return ans;}int fx[MAXN*10],fxtop=0,times=0;void dfs(int u,int fa,int pre){//  cout<<u<<endl;//    if(times==car)exit(0);  /*  if(u==t-1){        times++;        for(int i=1;i<=fxtop;i++)printf("%d %d\n",times,fx[i]);        return;    }*/    if(u!=1&&u<=p*q){         if(pre+1==u)printf("%d 1\n",times);//you         else printf("%d 0\n",times);         pre=u;    }  //  else if(fa+base==)    if(u==p*q)return;    for(int i=head[u];i;i=e[i].next){        int v=e[i].to,w=e[i].w;        if(v==fa)continue;          if(e[i^1].w&&i%2==0){          //  if(u-base+1==v)fx[++fxtop]=1;//出到入 入栈          //  else if(u-base+p==v)fx[++fxtop]=0;            dfs(v,u,pre);            e[i^1].w--;            break;        }    }}void bfs(int u, int fa) {    if(fa-base+1==u)printf("%d 1\n",times);//you    else if(fa-base+p==u)printf("%d 0\n",times);    //if (x == tot) return;    for (int i=head[u];i;i=e[i].next) {        if(e[i].to==fa)continue;        if (e[i^1].w&&i%2==0) { //只搜反向边             bfs(e[i].to, u);            e[i^1].w--;            return;        }    }}int main(){    int tem;    scanf("%d%d%d",&car,&p,&q);// p列  q 行    s=0,t=p*q+base+1;    for(int i=1;i<=q;i++){//拆点 出点+base        for(int j=1;j<=p;j++){ //1 2 3 4            int now=(i-1)*p+j;            scanf("%d",&ma[i][j]);                   //5 6 7 8 同行+1              add(now,now+base,INF,0);            if(ma[i][j]==2) add(now,now+base,1,1);        }    }    for(int i=1;i<=q;i++)        for(int j=1;j<=p;j++){            int now=(i-1)*p+j;            int lin_totodown=now+p;            int lin_totoright=now+1;            if(i!=q&&ma[i+1][j]!=1)add(now+base,lin_totodown,INF,0);            if(j!=p&&ma[i][j+1]!=1)add(now+base,lin_totoright,INF,0);        }    add(s,1,car,0);    add(p*q+base,t,car,0);    MC();    for(int i=1;i<=car;i++){        //if(times==car)break;      //  fxtop=0;memset(fx,0,sizeof(fx));        times++;        dfs(1,0,1);     }     return 0;}