1. 程式人生 > >BZOJ_1391_[Ceoi2008]order_最大權閉合子圖

BZOJ_1391_[Ceoi2008]order_最大權閉合子圖

desc tmp scrip break 機器 HR LG name ret

BZOJ_1391_[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
權值有正有負,很容易想到最大權閉合子圖。 S連工序-利潤代表割這個就不選這個任務 任務連機器-租的費用代表這個付出租的代價 機器連T-購買的費用代表付出購買的費用。 代碼:
#include <stdio.h>
#include <string.h>
#include <algorithm>
using namespace std;
#define N 2500
#define M 3600050
#define S (n+m+1)
#define T (n+m+2)
#define inf 100000000
int head[N],to[M],nxt[M],flow[M],cnt=1,dep[N],Q[N],l,r,sum,n,m,cur[N];
inline void add(int u,int v,int f) {
    to[++cnt]=v; nxt[cnt]=head[u]; head[u]=cnt; flow[cnt]=f;
    to[++cnt]=u; nxt[cnt]=head[v]; head[v]=cnt; flow[cnt]=0;
}
bool bfs() {
    int i;
    memset(dep,0,sizeof(dep)); l=r=0;
    Q[r++]=S; dep[S]=1;
    while(l<r) {
        int x=Q[l++];
        for(i=head[x];i;i=nxt[i]) {
            if(!dep[to[i]]&&flow[i]) {
                dep[to[i]]=dep[x]+1;
                if(to[i]==T) return 1;
                Q[r++]=to[i];
            }
        }
    }
    return 0;
}
int dfs(int x,int mf) {
    if(x==T) return mf;
    int nf=0,i;
    for(i=cur[x];i;i=nxt[i]) {
        if(dep[to[i]]==dep[x]+1&&flow[i]) {
            int tmp=dfs(to[i],min(mf-nf,flow[i]));
            if(!tmp) dep[to[i]]=0;
            nf+=tmp;
            flow[i]-=tmp;
            if(flow[i]) cur[x]=i;
            flow[i^1]+=tmp;
            if(nf==mf) break;
        }
    }
    return nf;
}
void dinic() {
    int ans=sum,f,i;
    while(bfs()) {
        for(i=1;i<=T;i++) cur[i]=head[i];
        while(f=dfs(S,inf)) ans-=f;
    }
    printf("%d\n",ans);
}
int main() {
    scanf("%d%d",&n,&m);
    int i,x,y,z,w;
    for(i=1;i<=n;i++) {
        scanf("%d%d",&x,&y);
        add(S,i,x);
        sum+=x;
        while(y--) {
            scanf("%d%d",&z,&w);
            add(i,z+n,w);
        }
    }
    for(i=1;i<=m;i++) {
        scanf("%d",&x);
        add(i+n,T,x);
    }
    dinic();
}

BZOJ_1391_[Ceoi2008]order_最大權閉合子圖