这道题不就是一个拆点后的少连边输出方案的方格取数嘛。方案可以用栈来输出。
有一点很妙的地方是,反向边走过以后就有值了,可以这样来打印方案。
为什么我打印方案那么多注释?周五做题真的有毒。。。我方案打印不对。。改了一年,距离写着题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;
}