火星探险问题

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

为什么我打印方案那么多注释？周五做题真的有毒。。。我方案打印不对。。改了一年，距离写着题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;
}