POJ 3422 Kaka's Matrix Travels（最小费用最大流）

http://poj.org/problem?id=3422

题意 ： 给你一个N*N的方格，每个格子有一个数字，让你从左上角开始走，只能往下往右走，走过的数字变为0，走K次，问最大能是多大，累加的。

思路 ：http://blog.csdn.net/qq172108805/article/details/7857503，这道题挺难的，我也是看的题解。

#include <iostream>
#include <stdio.h>
#include <string.h>
#include <queue>
#include <stdlib.h>
#include <algorithm>
 
using namespace std;
 
const int maxn =  ;
const int INF =  ;
const int maxm =  ;
 
struct node
{
    int u,v,c,f,b,next ;
}map[maxm] ;
bool flag[maxn] ;
int cap[][],dist[maxn],pre[maxn],head[maxn] ;
int s,t,n,k,ans,cnt ;
int del,val ;
 
void addedge(int u,int v,int c,int b)
{
    map[cnt].u = u ;
    map[cnt].v = v ;
    map[cnt].c = c ;
    map[cnt].b = b ;
    map[cnt].f =  ;
    map[cnt].next = head[u] ;
    head[u] = cnt++ ;
 
    map[cnt].u = v ;
    map[cnt].v = u ;
    map[cnt].c =  ;
    map[cnt].b = -b ;
    map[cnt].f =  ;
    map[cnt].next = head[v] ;
    head[v] = cnt++ ;
}
 
void spfa()
{
    memset(flag,,sizeof(flag)) ;
    queue<int >Q ;
    Q.push(s) ;
    flag[s] = true ;
    while(!Q.empty())
    {
        int x = Q.front() ;
        Q.pop() ;
        flag[x] = false ;
        for(int i = head[x] ; i+ ; i = map[i].next)
        {
            int y = map[i].v ;
            if(map[i].c > map[i].f && dist[y] < dist[x] + map[i].b)
            {
                dist[y] = dist[x]+map[i].b ;
                pre[y] = i ;
                if(!flag[y])
                {
                    flag[y] = true ;
                    Q.push(y) ;
                }
            }
        }
    }
}
 
void mcmf()
{
    for( ; ; )
    {
        memset(dist,-,sizeof(dist)) ;
        memset(pre,-,sizeof(pre)) ;
        dist[s] =  ;
        spfa() ;
        if(dist[t] == -)
        break ;
        int minn = INF ;
        for(int i = pre[t] ; i+ ; i = pre[map[i].u])
        minn = min(minn,map[i].c-map[i].f) ;
        for(int i = pre[t] ; i+ ; i = pre[map[i].u])
        {
            map[i].f += minn ;
            map[i^].f -= minn ;
        }
        ans += dist[t] ;
    }
}
 
void Init()
{
    cnt =  ;
    memset(head,-,sizeof(head)) ;
    s =  ;
    t = *n*n+ ;
    del = n*n ;
    ans =  ;
}
int main()
{
    scanf("%d %d",&n,&k) ;
    Init() ;
    for(int i =  ; i <= n ; i++)
    {
        for(int j =  ; j <= n ; j++)
        scanf("%d",&cap[i][j]) ;
    }
    addedge(s,,k,) ;
    addedge(n*n+del,t,k,) ;
    for(int i =  ; i <= n ; i++)
    {
        for(int j =  ; j <= n ; j++)
        {
            val = (i-)*n+j ;
            addedge(val,val+del,,cap[i][j]) ;
            addedge(val,val+del,INF,) ;
            if(j != n) addedge(val+del,val+,INF,) ;
            if(i < n) addedge(val+del,val+n,INF,) ;
        }
    }
    mcmf() ;
    printf("%d\n",ans) ;
    return ;
}