HDU2686 Matrix——最大费用最大流

点这里

题意： 有一个n*n的矩阵，矩阵中每一个方格都有一个权值，现要求你从左上角的点走到右下角的点再走回左上角，除了左上角的点，其他的点均不能重复经过，收益为所有经过的点的权值。求最大收益。
题解： 就是最小费用最大流变成了最大费用，还是能很容易看出端倪的。因为每个点只能经过一次，而每个点可能连接着2或4个点；所以为了限流，将每个点x拆分成x和x’，x连接这他左上的点，x与x’连接，x’连接右下的点。这样起到了控制每个点只经过一次的作用。但是左上角的点和右下角的点是不止经过一次的，所以额外增加一条费用为0、流量为1的边。

---

过程中犯的错：

• 没有拆点： 做太急了，都忘记需要拆点限流了。
• 额外边： 前面有解释原因。

---

#include<algorithm>
#include<iostream>		
#include<cstring>
#include<vector>
#include<queue>
#include<cmath>
using namespace std;
const int N = 2010;
const int inf = 0x3f3f3f3f;

int n, a[40][40];
int inq[N], dis[N], pre[N], preve[N];	//inq记录是否入队、dis记录起点到i的最小距离、pre记录前驱节点、preve最短路上连接这个点的边的编号
struct edge{
	int to, cost, cap, rev;	//rev用于记录反边的编号
	edge(int a, int b, int c, int d){	to = a; cost = b; cap = c; rev = d;}
};
vector<edge> e[N];
void addedge(int u, int v, int cost, int cap){
	e[u].push_back(edge(v, cost, cap, e[v].size()));
	e[v].push_back(edge(u, -cost, 0, e[u].size() - 1));	//加反边
}
bool spfa(int s, int t, int cnt){
	for(int i = 0; i <= cnt; i++){	pre[i] = -1; dis[i] = -1; inq[i] = 0;}
	queue<int> Q;
	Q.push(s); inq[s] = 1; dis[s] = 0;
	while(!Q.empty()){
		int u = Q.front(); Q.pop();	inq[u] = 0;
		for(int i = 0; i < e[u].size(); i++)
			if(e[u][i].cap > 0){
				int v = e[u][i].to, cost = e[u][i].cost;
				if(dis[v] < dis[u] + cost){
					dis[v] = dis[u] + cost;
					pre[v] = u, preve[v] = i;
					if(!inq[v]){	inq[v] = 1;	Q.push(v);}
				}
			}
	}
	return dis[t] != -1;
}
int mincost(int s, int t, int cnt){
	int cost = 0;
	while(spfa(s, t, cnt)){
		int v = t, flow = inf;
		while(pre[v] != -1){	//回溯整个最短路径
			int u = pre[v], i = preve[v];	//u是v的前驱节点，u的第i个边连接v
			flow = min(flow, e[u][i].cap);
			v = u; 
		}
		v = t;
		while(pre[v] != -1){
			int u = pre[v], i = preve[v];
			e[u][i].cap -= flow;
			e[v][e[u][i].rev].cap += flow;	//反边更新残留容量
			v = u; 
		}
		cost += dis[t] * flow;	//费用累加
	}
	return cost;
}
int main(){
	while(~scanf("%d", &n)){
		for(int i = 0; i < N; i++)	e[i].clear();
		for(int i = 1; i <= n; i++)	for(int j = 1; j <= n; j++)	scanf("%d", a[i] + j);
		for(int i = 1; i <= n; i++)
			for(int j = 1; j <= n; j++){
				int fir = (i - 1) * n + j, sec = fir + n * n;
				addedge(fir, sec, a[i][j], 1);
				if(j + 1 <= n)
					addedge(sec, fir + 1, 0, 1);
				if(i + 1 <= n)
					addedge(sec, fir + n, 0, 1);
			}
		addedge(1, 1 + n * n, 0, 1);
		addedge(n * n, 2 * n * n, 0, 1);
		printf("%d\n", mincost(1, 2 * n * n, 2 * n * n));
	}
	return 0;
}