【bzoj1391】[Ceoi2008]order

Description

有N个工作，M种机器，每种机器你可以租或者买过来. 每个工作包括若干道工序，每道工序需要某种机器来完成,你可以通过购买或租用机器来完成。 现在给出这些参数，求最大利润
Input

第一行给出 N,M(1<=N<=1200,1<=M<=1200) 下面将有N块数据，每块数据第一行给出完成这个任务能赚到的钱(其在[1,5000])及有多少道工序 接下来若干行每行两个数，分别描述完成工序所需要的机器编号及租用它的费用(其在[1,20000]) 最后M行，每行给出购买机器的费用(其在[1,20000])
Output

最大利润
Sample Input

2 3

100 2

1 30

2 20

100 2

1 40

3 80

50

80

110
Sample Output

50
HINT

题解
对于工作和机器建点，由S向每一个工作连边(S,u,a)a为完成工作的奖励，由一个工作向所需机器连边(u,v,b)b为租用机器费用，由每个机器向T连边(v,T,c)c为购买机器的费用。

　 求最小割，构图提供的选择有：不做该工作，租用一个工序需要的机器，购买一个机器满足所有需要。最小割保证ST不相连，即保证给所有工作分配了一种方案：完成所有工序或不做该工作，且该方案花费最少，最小割即为最小亏损。

代码isap

#include<bits/stdc++.h>
#define oo 100000000
#define ll long long
#define mo 1000000007
#define N 50005
using namespace std;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while (ch<'0'||ch>'9'){if (ch=='-') f=-1;ch=getchar();}
    while (ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
int tot,S,T,sum,ans,n,m;
int Head[2500],ret[3000001],Next[3000001],len[3000001],dis[2500],gap[2500];
inline void ins(int u,int v,int l)
{
    tot++;
    ret[tot]=v;len[tot]=l;
    Next[tot]=Head[u];Head[u]=tot;
}
int isap(int u,int aug)
{
    if (u==T) return aug;
    if (aug==0) return 0;
    int min_d=T,flow=0;
    for (int i=Head[u];i;i=Next[i])
    {
        int v=ret[i];
        if (len[i]>0)
        {
            if (dis[u]==dis[v]+1)
            {
                int t=isap(v,min(aug-flow,len[i]));
                len[i]-=t;
                len[i^1]+=t;
                flow+=t;
            }
            if (flow==aug) return flow;
            if (dis[S]==T) return flow;
            min_d=min(min_d,dis[v]+1);
        }
    }
    if (flow==0)
    {
        gap[dis[u]]--;
        if (gap[dis[u]]==0) dis[S]=T;
        dis[u]=min_d;
        gap[dis[u]]++;
    }
    return flow;
}
int main()
{
    n=read();m=read();S=n+m+1;T=S+1;tot=1;
    for (int i=1;i<=n;i++)
    {
        int val=read(),t=read();
        ins(S,i,val);ins(i,S,0);sum+=val;
        while (t--)
        {
            int x=read()+n;val=read();
            ins(i,x,val);ins(x,i,0);
        }
    }
    for (int i=1;i<=m;i++)
    {
        int val=read();
        ins(i+n,T,val);ins(T,i+n,0);
    }
    while (dis[S]!=T)
    {
        ans+=isap(S,oo);
    }
    printf("%d",sum-ans);
    return 0;
}