codevs1227（方格取数费用流）

给出一个n*n的矩阵,每一格有一个非负整数Aij,(Aij <= 1000)现在从(1,1)出发,可以往右或者往下走,最后到达(n,n),每达到一格,把该格子的数取出来,该格子的数就变成0,这样一共走K次,现在要求K次所达到的方格的数的和最大

控制流量，每一个格子最多走一次，使费用最大，可以发现这是费用流。

首先对于每一个格子，我们拆成两个点（因为要限制流量）；每一个格子可以取一次，但是每一个格子是可以走多次的，那么我们在两个点中建两种边，

1：费用为权值，流量为一，2：费用为0，流量为inf。分别代表两种情况。

又写炸了。。细节注释中

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cstdlib>
#include<cmath>
#include<iostream>
#include<queue>
#define debug(x) cout<<#x<<"="<<x<<endl
using namespace std;
const int N=5050;
const int inf=0x3f3f3f3f;
struct aa
{
	int pre,to,cap,flow,w;
}edge[N*20];
int head[N],tot,id1[55][55],id2[55][55];//计算好数组大小
int n,k,s,t,v[55][55],dis[N],now[N],pre[N];
bool b[N];
void addedge(int x,int y,int z,int w)
{
	edge[++tot].to=y;edge[tot].cap=z;edge[tot].w=w;
	edge[tot].pre=head[x];head[x]=tot;
	edge[++tot].to=x;edge[tot].cap=0;edge[tot].w=-w;
	edge[tot].pre=head[y];head[y]=tot;
}
bool spfa()
{
	memset(dis,-inf,sizeof(dis));//跑最大费用最大流，初始化为最小
	memset(b,false,sizeof(b));
	memset(now,0,sizeof(now));
	memset(pre,0,sizeof(pre));
	b[s]=true;dis[s]=0;
	queue<int> q;q.push(s);
	while (!q.empty())
	{
		int u=q.front();q.pop();
		for (int i=head[u];i;i=edge[i].pre)
		if (edge[i].cap>edge[i].flow&&
		dis[edge[i].to]<dis[u]+edge[i].w)
		{
			dis[edge[i].to]=dis[u]+edge[i].w;
			pre[edge[i].to]=u;now[edge[i].to]=i;
			if (!b[edge[i].to])
			{
				b[edge[i].to]=true;
				q.push(edge[i].to);
			}
		}
		b[u]=false;
	}
	if (dis[t]>0) return true;//必须有大于0的路径才ok
	return false;
}
int work()
{
	int ans=0;
	while (spfa())
	{
		int flow=inf;
		for (int i=t;i!=s;i=pre[i])
		flow=min(flow,edge[now[i]].cap-edge[now[i]].flow);
		ans+=flow*dis[t];
		for (int i=t;i!=s;i=pre[i])
		edge[now[i]].flow+=flow,
		edge[((now[i]-1)^1)+1].flow-=flow;//刚开始，把now【i】打成i了，查了半天
	}
	return ans;
}
int main()
{
	scanf("%d%d",&n,&k);
	s=0;t=n*n*2+1;
	int x,bj=0;
	for (int i=1;i<=n;i++)
		for (int j=1;j<=n;j++) 
		id1[i][j]=++bj,id2[i][j]=++bj;
	
	for (int i=1;i<=n;i++)
		for (int j=1;j<=n;j++) 
		{
			scanf("%d",&x);
			addedge(id1[i][j],id2[i][j],1,x);
			addedge(id1[i][j],id2[i][j],inf,0);
			if (i>1) addedge(id2[i-1][j],id1[i][j],inf,0);
			if (j>1) addedge(id2[i][j-1],id1[i][j],inf,0);
		}
	addedge(s,id1[1][1],k,0);
	addedge(id2[n][n],t,k,0);
	printf("%d",work());
	return 0;
} 

总结

1：费用流当中，在两个点之间加两种不同的边，来代表不同的情况，且两种边的费用不同，由此又能限制流的顺序，从而限制情况选择的顺序。

这是费用流当中，关于两个点之间不同方案的常用方法

本题当中是在，分点的入点和出点之间的两种情况建的边