摘取作物

题目

Feather的农场里有N*M块地，排列成N行，每行M块地。Feather在每块地里种植了不同的农作物。现在这些农作物都成熟了，可以摘取下来出售了。其中第i行第j列的地里的农作物的价值为W[i,j]。
JackRabbit是Feather的好友，平时经常为Feather的农作物除草除虫。为了答谢JackRabbit，Feather决定把一部分农作物送给JackRabbit。JackRabbit很高兴，恨不得一下子把农场里的农作物摘空。
为了防止JackRabbit把农作物摘空，Feather提出了两个条件：
1.每行最多选取两块地；
2.每列最多选取两块地。
这下子把JackRabbit难住了。如何在满足这两个条件的前提下，使得摘取的农作物的价值之和最大呢？

分析

费用流的模型题，
先考虑构图，

对于每一行，用一个点表示，因为每行最多选两块地，所以从源点s向这n个行点连一条容量为2，费用为0的边。
对于每一列，同样用一个点表示，因为每列最多选两块地，所以从这m个列点向汇点t连一条容量为2，费用为0的边。
最后，因为在同一行、同一列最多只有一块的，所以，从n个行点向m个列点连一条容量为1，费用为W[i,j]的边。
做一遍费用流就可以了。

#include <cmath>
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <queue>
const int maxlongint=2147483647;
const int mo=1000000007;
const int N=50005;
using namespace std;
int f[400][400],cost[400][400],next[5000],dis[5000],d[1000000],n,m,ans;
bool v[5000];
bool spfa()
{
    for(int i=0;i<=n+m+1;i++)
        dis[i]=-maxlongint;
    memset(next,0,sizeof(next));
    memset(v,true,sizeof(v));
    int head=0,tail=1,k;
    dis[0]=0;
    d[0]=0;
    while(head<tail)
    {
        k=d[++head];
        v[k]=true;
        for(int i=0;i<=n+m+1;i++)
        {
            if(f[k][i] && dis[i]<dis[k]+cost[k][i])
            {
                dis[i]=dis[k]+cost[k][i];
                next[i]=k;
                if(v[i])
                {
                    d[++tail]=i;
                    v[i]=false;
                }
            }
        }
    }
    if(dis[n+m+1]<=0) return false;
        else return true; 
}
int solve()
{
    int sum=maxlongint;
    for(int i=n+m+1;i;i=next[i])
    {
        sum=min(sum,f[next[i]][i]);
    }
    for(int i=n+m+1;i;i=next[i])
    {
        f[next[i]][i]-=sum;
        f[i][next[i]]+=sum;
        ans+=cost[next[i]][i]*sum;
    }
}
int main()
{
    freopen("pick.in","r",stdin);
    freopen("pick.out","w",stdout);
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++) f[0][i]=2;
    for(int i=1;i<=m;i++) f[n+i][n+m+1]=2;
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=m;j++)
        {
            scanf("%d",&cost[i][n+j]);
            f[i][n+j]=1;
            cost[i][n+j]=cost[i][n+j];
            cost[n+j][i]=-cost[i][n+j];
        }
    }
    ans=0;
    while(spfa()) solve();
    printf("%d",ans);
}