其他分享
首页 > 其他分享> > P3159 [CQOI2012]交换棋子(费用流)

P3159 [CQOI2012]交换棋子(费用流)

作者:互联网

题目描述

有一个n行m列的黑白棋盘,你每次可以交换两个相邻格子(相邻是指有公共边或公共顶点)中的棋子,最终达到目标状态。要求第i行第j列的格子只能参与mi,j次交换。

输入输出格式

输入格式:

 

第一行包含两个整数n,m(1<=n, m<=20)。以下n行为初始状态,每行为一个包含m个字符的01串,其中0表示黑色棋子,1表示白色棋子。以下n行为目标状态,格式同初始状态。以下n行每行为一个包含m个0~9数字的字符串,表示每个格子参与交换的次数上限。

 

输出格式:

 

输出仅一行,为最小交换总次数。如果无解,输出-1。

 

输入输出样例

输入样例#1: 复制
3 3
110
000
001
000
110
100
222
222
222
输出样例#1: 复制
4





好像所有的没有思路的题都能用毒瘤的网络流来搞一下,
这题真的毒瘤,第一次见拆三个点的
看到一个点最多参与几次交换,就应该想到拆点限流,然后最后问你最先的代价,应该可以联想到费用这一块,
那就向费用流这块搞啊,这题的建模是真的难想,所以说就不详细写了,细节挺多的。
然后就是不用担心一个1到达目标的时候,打乱的别的1,因为这时通过交换两个点的先后顺序来避免这个问题。
放个网址: https://www.luogu.org/problemnew/solution/P3159



  1 #include<queue>
  2 #include<cstdio>
  3 #include<cstring>
  4 #include<iostream>
  5 #include<algorithm>
  6 #define MAXN 2010
  7 #define MAXM 64010 
  8 #define INF 0x3f3f3f3f
  9 #define fpop(x) x.front();x.pop()
 10 using namespace std;
 11 
 12 int pre_node[MAXN],pre_edge[MAXN];
 13 
 14 char ch,mp_bg[25][25],mp_ed[25][25];
 15 
 16 int n,m,cnt=-1,dis[MAXN],vis[MAXN],flow[MAXN],head[MAXN],maxf[25][25];
 17 
 18 struct edge{
 19     int nxt,to,w,f;
 20 }e[MAXM];
 21 
 22 inline int _bg(int x,int y){return n*m*0+(x-1)*m+y;}//起始点[x,y]的编号
 23 inline int _ed(int x,int y){return n*m*1+(x-1)*m+y;}//目标点[x,y]的编号 
 24 inline int _inn(int x,int y){return n*m*2+(x-1)*m+y;}//棋盘[x,y]的Inn点编号 
 25 inline int _mid(int x,int y){return n*m*3+(x-1)*m+y;}//棋盘[x,y]的mid点标号
 26 inline int _out(int x,int y){return n*m*4+(x-1)*m+y;}//棋盘[x,y]的Out点编号 
 27 
 28 inline bool in_map(int x,int y){
 29     return 1<=x && x<=n && 1<=y && y<=m;
 30 }//判断是否越界 
 31 
 32 inline void add_edge(int from,int to,int flw,int val){
 33     e[++cnt].nxt=head[from];
 34     e[cnt].to=to;
 35     e[cnt].f=flw;
 36     e[cnt].w=val;
 37     head[from]=cnt;
 38 }
 39 
 40 queue<int>que;
 41 
 42 inline bool spfa(int s,int t){
 43     memset(vis,0,sizeof(vis));
 44     memset(dis,0x3f,sizeof(dis));
 45     memset(flow,0x3f,sizeof(flow));
 46     que.push(s); vis[s]=true; dis[s]=0;
 47     while(!que.empty()){
 48         int u=fpop(que);
 49         for(int i=head[u];~i;i=e[i].nxt){
 50             int v=e[i].to;
 51             if(dis[v]>dis[u]+e[i].w && e[i].f){
 52                 dis[v]=dis[u]+e[i].w;
 53                 flow[v]=min(flow[u],e[i].f);
 54                 pre_node[v]=u;
 55                 pre_edge[v]=i;
 56                 vis[v]=true;que.push(v);
 57             }
 58         }
 59     }
 60     return dis[t]!=INF;
 61 }
 62 
 63 int mv[8][2]={{1,0},{-1,0},{0,1},{0,-1},{1,1},{1,-1},{-1,1},{-1,-1}};
 64 
 65 int main(){
 66     memset(head,-1,sizeof(head));
 67     scanf("%d%d",&n,&m);
 68     for(int i=1;i<=n;++i){
 69         for(int j=1;j<=m;++j){
 70             scanf(" %c",&mp_bg[i][j]);
 71         }
 72     }
 73     for(int i=1;i<=n;++i){
 74         for(int j=1;j<=m;++j){
 75             scanf(" %c",&mp_ed[i][j]);
 76         }
 77     }
 78     for(int i=1;i<=n;++i){
 79         for(int j=1;j<=m;++j){
 80             scanf(" %c",&ch);
 81             maxf[i][j]=ch-'0';
 82             //最大经过次数 
 83         }
 84     }
 85     //输入起始态和目标态棋盘
 86     int s=0,t=n*m*5+1;
 87     int cnt_1=0,cnt_2=0;
 88     for(int i=1;i<=n;++i){
 89         for(int j=1;j<=m;++j){
 90             if(mp_bg[i][j]==mp_ed[i][j]){
 91                 add_edge(_inn(i,j),_mid(i,j),maxf[i][j]/2,0);
 92                 add_edge(_mid(i,j),_inn(i,j),000000000000,0);
 93 
 94                 add_edge(_mid(i,j),_out(i,j),maxf[i][j]/2,0);
 95                 add_edge(_out(i,j),_mid(i,j),000000000000,0);
 96             }else{
 97                 if(mp_bg[i][j]=='1'){
 98                     add_edge(_inn(i,j),_mid(i,j),(maxf[i][j]+0)/2,0);
 99                     add_edge(_mid(i,j),_inn(i,j),000000000000,0);
100 
101                     add_edge(_mid(i,j),_out(i,j),(maxf[i][j]+1)/2,0);
102                     add_edge(_out(i,j),_mid(i,j),000000000000,0);
103                 }
104                 if(mp_ed[i][j]=='1'){
105                     add_edge(_inn(i,j),_mid(i,j),(maxf[i][j]+1)/2,0);
106                     add_edge(_mid(i,j),_inn(i,j),000000000000,0);
107 
108                     add_edge(_mid(i,j),_out(i,j),(maxf[i][j]+0)/2,0);
109                     add_edge(_out(i,j),_mid(i,j),000000000000,0);
110                 }
111             }
112 
113             if(mp_bg[i][j]=='1'){
114                 ++cnt_1;
115                 //连接源点到初始点     f=1 w=0;
116                 add_edge(s,_mid(i,j),1,0);
117                 add_edge(_mid(i,j),s,0,0);
118                 //连接起始点到棋盘 
119               //  add_edge(_bg(i,j),_mid(i,j),1,0);
120               //  add_edge(_mid(i,j),_bg(i,j),0,0);
121             }
122             if(mp_ed[i][j]=='1'){
123                 ++cnt_2;
124                 //连接终结点到汇点     f=1 w=0; 
125                 add_edge(_mid(i,j),t,1,0);
126                 add_edge(t,_mid(i,j),0,0);
127                 //连接棋盘到终结点 
128               //  add_edge(_mid(i,j),_ed(i,j),1,0);
129                // add_edge(_ed(i,j),_mid(i,j),0,0);
130             }
131             //棋盘的八连通边  f=INF w=1;
132             for(int k=0;k<8;++k){
133                 int ni=i+mv[k][0];
134                 int nj=j+mv[k][1];
135                 if(in_map(ni,nj)){
136                     //从点[i,j]的out连接点[ni,nj]的inn 
137                     add_edge(_out(i,j),_inn(ni,nj),INF,+1);
138                     add_edge(_inn(ni,nj),_out(i,j),000,-1);
139                 }
140             }
141         }
142     }
143     //棋子数变动->No solution 
144     if(cnt_1!=cnt_2){
145         puts("-1");
146         return 0;
147     }
148     //然后跑费用流 
149     int max_flow=0,min_cost=0;
150     while(spfa(s,t)){
151         max_flow+=flow[t];
152         min_cost+=flow[t]*dis[t];
153         int u=t;
154         while(u!=s){
155             e[pre_edge[u]^0].f-=flow[t];
156             e[pre_edge[u]^1].f+=flow[t];
157             u=pre_node[u];
158         }
159     }
160     if(max_flow!=cnt_1){
161         puts("-1");
162         return 0;
163     }
164     printf("%d\n",min_cost);
165 }

 




标签:pre,25,return,int,P3159,flow,棋子,CQOI2012,dis
来源: https://www.cnblogs.com/zhangbuang/p/10638102.html