【HDU 2014 Multi-University Training Contest 1 1002】/【HDU 4862】Jump

多校训练就这么华丽丽的到了 ，于是乎各种华丽丽的被虐也开始了。

这是多校的1002; 最小费用最大流。

题目大意：

有n*m个方格，每个方格都一个的十进制一位的数。你可以操作K次。

对于每一次操作，你可以选择一个出发点向下或向右Jump。跳的花费是|x1-x2|+|y1-y2|-1的能量 。如果你跳的这两个位置上数字相同，那么你就会获得数字表示的能量值。

对于每一次操作，你可以这样跳任意次 ，但是每个位置只能经过一次在这K次操作中。

初始能量值是0，当操作完成后，如果n*m个方格没有都经过过，输出“-1”，否则输出可以得到的最大能量值。

解题思路：

建立一个流量网络，一个二部图。X部分向Y部分链接的情况表示可以从一个点跳到另一个点，超级源点和超级汇点分别同X部分的点和Y部分的点链接。在X部分中多加一个点它与源点的流量是K费用是0，与Y部分所有点链接流量是1费用是0。这表示操作K次。

下面是代码：

#include <stdio.h>
#include <string.h>
#include <iostream>
#include <algorithm>
#include <vector>
#include <queue>
#include <set>
#include <map>
#include <string>
#include <math.h>
#include <stdlib.h>
#define clear(A, X, SIZE) memset(A, X, sizeof(A[0]) * (SIZE))
#define clearall(A, X) memset(A, X, sizeof(A))
#define memcopy1(A , X, SIZE) memcpy(A , X ,sizeof(X[0])*(SIZE))
#define memcopy1all(A, X) memcpy(A , X ,sizeof(X))
#define max( x, y )  ( ((x) > (y)) ? (x) : (y) )
#define min( x, y )  ( ((x) < (y)) ? (x) : (y) )

using namespace std;

struct node
{
    int u,v,c,f,next;
} edge[21000];
char s[12][15];
const int inf=1<<30;
int head[210],cnt,dis[210],pre[210],n,m,cost,flow;
bool vis[210];
bool spfa()
{
    int i,u;
    clearall(pre,-1);
    clearall(dis,0x3f3f3f);
    clearall(vis,false);
    queue <int>q;
    dis[n*m*2]=0;
    vis[n*m*2]=true;
    q.push(n*m*2);
    while(!q.empty())
    {
        u=q.front();
        q.pop();
        i=head[u];
        vis[u]=false;
        while(i!=-1)
        {
            if(edge[i].f>0&&dis[edge[i].v]>dis[u]+edge[i].c)
            {
                dis[edge[i].v]=dis[u]+edge[i].c;
                pre[edge[i].v]=i;
                if(!vis[edge[i].v])
                {
                    vis[edge[i].v]=true;
                    q.push(edge[i].v);
                }
            }
            i=edge[i].next;
        }
    }
    if(pre[2*n*m+1]==-1)return false;
    else return true;
}
void does()
{
    cost=0;
    flow=0;
    while(spfa())
    {
        int max1=inf;
        int p=pre[2*n*m+1];
        while(p!=-1)
        {
            max1=min(max1,edge[p].f);
            p=pre[edge[p].u];
        }
        p=pre[2*n*m+1];
        while(p!=-1)
        {
            edge[p].f-=max1;
            edge[p^1].f+=max1;
            cost+=max1*edge[p].c;
            p=pre[edge[p].u];
        }
        flow+=max1;
    }
}
void addedge(int u,int v,int f,int c)
{
    edge[cnt].u=u;
    edge[cnt].v=v;
    edge[cnt].f=f;
    edge[cnt].c=c;
    edge[cnt].next=head[u];
    head[u]=cnt++;
    edge[cnt].u=v;
    edge[cnt].v=u;
    edge[cnt].f=0;
    edge[cnt].c=-c;
    edge[cnt].next=head[v];
    head[v]=cnt++;
}
int main()
{
    int T,k,case1=1,u,v,f,c;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d%d",&n,&m,&k);
        for(int i=0; i<n; i++)
        {
            scanf("%s",s[i]);
        }
        clearall(head,-1);
        cnt=0;
        for(int i=0; i<n; i++)
        {
            for(int j=0; j<m; j++)
            {
                u=i*m+j;
                for(int k=j+1; k<m; k++) //right
                {
                    v=i*m+k+n*m;
                    c=k-j-1;
                    c=-c;
                    if(s[i][j]==s[i][k])c+=s[i][j]-'0';
                    addedge(u,v,1,-c);
                }
                for(int k=i+1; k<n; k++) //down
                {
                    v=k*m+j+n*m;
                    c=k-i-1;
                    c=-c;
                    if(s[i][j]==s[k][j])c+=s[i][j]-'0';
                    addedge(u,v,1,-c);
                }
            }
        }
        for(int i=0; i<n; i++)
        {
            for(int j=0; j<m; j++)
            {
                addedge(2*n*m,i*m+j,1,0);
                addedge(i*m+j+n*m,2*n*m+1,1,0);
                addedge(2*n*m+2,i*m+j+n*m,1,0);
            }
        }
        addedge(2*n*m,2*n*m+2,k,0);
        does();
        printf("Case %d : ",case1++);
        if(flow!=n*m)puts("-1");
        else printf("%d\n",-cost);
    }
    return 0;
}