zoj 3308 Move the Knights

最小费用流算法实践

最新推荐文章于 2022-02-25 19:38:17 发布

原创最新推荐文章于 2022-02-25 19:38:17 发布 · 340 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#iostream #struct #c

本文介绍了一个关于最小费用流算法的应用实例，通过解决一个具体的骑士移动问题来展示算法的实现细节。文章提供了一段完整的C++代码，详细解释了如何构建图、添加边以及进行最小费用最大流的计算。

构图，最小费用流,直接套用。

wa了一把，题目没看清，最后一个是max( P_(row1,col1) , P_(row2,col2) ) 当做了，max(P_(row1,col1) x P_(row2,col2), P_(row1,col1) + P_(row2,col2) ).悲剧!!!

#include <cstdio> #include <cstdlib> #include <iostream> #include <string.h> #include <queue> #include <limits.h> using namespace std; #define MAX 230 #define Max(a,b) ( (a>b) ? (a):(b) ) typedef struct NODE { int from,to,cap,cost; int next; } NODE; NODE node[MAX*MAX]; int head[MAX]; int R,C,N,K,cou,sufx,flow; int g[16][16]; int dir[8][2]={ {-2,-1},{-1,-2}, {1,-2},{2,-1}, {2,1},{1,2}, {-1,2},{-2,1} }; void Add(int u,int v,int cap,int cost) { node[cou].from=u; node[cou].to=v; node[cou].cap=cap; node[cou].cost=cost; node[cou].next=head[u]; head[u]=cou++; node[cou].from=v; node[cou].to=u; node[cou].cap=0; node[cou].cost=-cost; node[cou].next=head[v]; head[v]=cou++; } int MincostMaxflow(int s,int t) { queue<int> q; int inq[MAX],pre[MAX],dis[MAX],re[MAX]; int u,v,i,a,c=0,ind,cost,cap; while(1) { memset(inq,0,sizeof(inq)); for(i=s;i<=MAX;i++) dis[i] = INT_MAX; dis[s]=0; inq[s]=1; pre[s]=s; q.push(s); while( !q.empty() ) { u = q.front(); q.pop(); inq[u]=0; ind = head[u]; while(ind!=0) { u = node[ind].from; v = node[ind].to; cost = node[ind].cost; cap = node[ind].cap; if(cap>0 && dis[v]>dis[u]+cost) { dis[v] = dis[u] + cost; pre[v] = u; re[v] = ind; if(!inq[v]) { q.push(v); inq[v] = 1; } } ind = node[ind].next; } } if(dis[t]==INT_MAX) break; a = INT_MAX; for(u=t;u!=s;u=pre[u]) if(node[re[u]].cap<a) a = node[re[u]].cap; for(u=t;u!=s;u=pre[u]) { node[re[u]^1].cap += a; node[re[u]].cap -=a; } flow+=a; c += dis[t]*a; } return c; } void addEdge(int type,int x,int y) { Add(1,x*C+y+sufx,1,0);//从s'点到其中一个骑士所在点 int i,tmpx,tmpy,to,from; for(i=0;i<8;i++) { tmpx=x+dir[i][0]; tmpy=y+dir[i][1]; if(tmpx>=1 && tmpx<=R && tmpy>=1 && tmpy<=C) { from = x*C+y+sufx; to = tmpx*C+tmpy+sufx; if(head[to]==0) { Add(to,2,1,0); } switch(type) { case 1: { Add(from,to,1,g[x][y]*g[tmpx][tmpy]); break; } case 2: { Add(from,to,1,g[x][y]+g[tmpx][tmpy]); break; } default:{ Add(from,to,1,Max(g[x][y],g[tmpx][tmpy]) ); } } } } } void init() { memset(head,0,sizeof(head)); memset(node,'/0',sizeof(node)); sufx = 2-C; cou=2; flow=0; Add(0,1,K,0); } int main() { int i,j,ans,x,y,type; while(scanf("%d%d%d%d",&R,&C,&N,&K)!=EOF) { init(); for(i=1;i<=R;i++) for(j=1;j<=C;j++) scanf("%d",&g[i][j]); for(i=1;i<=N;i++) { scanf("%d%d%d",&type,&x,&y); addEdge(type,x,y); } ans = MincostMaxflow(0,2); if(flow!=K) printf("-1/n"); else printf("%d/n",ans); } return 0; }