bzoj1570 [CEOI2008]order

最小割 想法比较简单
我们对于每个任务 由S向他连容量为收益的边
对于每台机器 由它向T连容量为价格的边
然后对于每个任务向它的每个工序所需的机器连容量为租用机器价格的边
那么答案即为所有收益减去最(zui)小(da)割(liu)

这道题本身很水 但是反映了最小割的一种思维方法
对于一个任务 如果放弃 那么就要割掉源点向它连去的边 损失掉它的收益
如果要做 那么要么购买机器 割掉机器连到汇点的边 损失机器的价格 要么租用机器 损失掉租用机器的成本

最小割的思维方式很有意思对不对OAO

#include <cstdio>
#include <algorithm>
#include <cstring>

using namespace std;

const int N=3000,M=3000005,S=N-1,T=N-2,INF=1e9;

int n=0,m=0;
int head[N],next[M],to[M],c[M],f[M],edge=0;
int sum=0;
int front=0,back=0,q[M],d[N],cur[N];

inline void addEdge(int u,int v,int ac) {
    to[edge]=v,c[edge]=ac,next[edge]=head[u],head[u]=edge++;
    to[edge]=u,c[edge]=0,next[edge]=head[v],head[v]=edge++; 
}

inline void build() {
    for (int i=1;i<=n;++i) {
        int profit=0,k=0;
        scanf("%d %d",&profit,&k);
        sum+=profit;
        addEdge(S,i,profit);
        for (int j=1;j<=k;++j) {
            int num=0,price=0;
            scanf("%d %d",&num,&price);
            addEdge(i,n+num,price);
        }
    }
    for (int i=1;i<=m;++i) {
        int p=0;
        scanf("%d",&p);
        addEdge(n+i,T,p);
    }
}

inline int bfs() {
    front=back=0;
    memset(d,-1,sizeof(d)); 
    q[back++]=S;
    d[S]=1;
    while (front<back) {
        int x=q[front++];
        for (int e=head[x];~e;e=next[e]) {
            int v=to[e];
            if (d[v]==-1 && c[e]>f[e]) {
                d[v]=d[x]+1;
                q[back++]=v;
            }
        }
    }
    return d[T]!=-1;
}

int dfs(int x,int t) {
    if (x==T || t==0)
        return t;
    int flow=0;
    for (int e=cur[x];~e;e=next[e]) {
        int v=to[e];
        if (c[e]>f[e] && d[v]==d[x]+1) {
            int f0=dfs(v,min(t,c[e]-f[e]));
            flow+=f0;
            t-=f0;
            f[e]+=f0;
            f[e^1]-=f0;
            if (c[e]>f[e])
                cur[x]=e;
            if (!t)
                return flow;
        }

    }
    if (t)
        d[x]=-1;
    return flow;
}

inline int dinic() {
    int flow=0;
    while (bfs()) {
        for (int i=1;i<=S;++i)
            cur[i]=head[i];
        flow+=dfs(S,INF);
    }
    return flow;    
}

int main(void) {
    scanf("%d %d",&n,&m);
    memset(head,-1,sizeof(head));
    build();
    printf("%d\n",sum-dinic());
    return 0;
}

不常在bzoj做题 0msTLE我还以为是我自己没写当前弧优化把这题给T了
加了优化还是T
后来wxy神犇和我说 bzoj0msT其实是MLE导致的
我就看了一下黄学长数组开了多大 抄了过来
果然直接过了= =#

我要是再手算网络流边数我就是狗 我这句话就放在这了= =