bzoj1391 [Ceoi2008]order 网络流——最大权闭合子图

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

题目分析：
最大权闭合子图模型。
从源点向工作连容量为收益的边。
从机器向汇点连容量为购买的边。
从工作向机器连容量为租用的边。
用总收益减去最小割（即最大流）即可。

这题不仅卡时，而且卡内存（我会告诉你们我交了3遍RE，3遍MLE，3遍TLE么orz……
读入优化，当前弧优化，时间戳……我尽力了Orz
代码如下：

#include <cstdio>
#include <iostream>
#define N 1600000
using namespace std;
namespace IStream
{
    char get_char()
    {
        const int L=1<<15;
        static char buffer[L];
        static char *C,*mempool;
        if(C==mempool) mempool=(C=buffer)+fread(buffer,1,L,stdin);
        return *C++;
    }
    int get_int()
    {
        int re=0;
        char c;
        do c=get_char();while(c<'0' || c>'9');
        while(c>='0' && c<='9') re=(re<<1)+(re<<3)+c-'0',c=get_char();
        return re;
    }
}
const int INF=0x3f3f3f3f;
inline int Min(int x,int y) { return x<y?x:y; }
int n,m,ans=0,S,T;
int Time_Stamp=0;
int fir[2500],nes[N<<1],v[N<<1],q[N<<1],tot=1;
int d[2500],dl[2500],cur[2500],tim[2500];
void edge(int x,int y,int z)
{
    v[++tot]=y;
    q[tot]=z;
    nes[tot]=fir[x];
    fir[x]=tot;
}
#define edge(x,y,z) edge(x,y,z),edge(y,x,0)
bool bfs()
{
    static int c;
    Time_Stamp++;
    int l=1,r=1;
    d[S]=1; dl[1]=S;
    tim[S]=Time_Stamp;
    while(l<=r)
    {
        c=dl[l++];
        for(int t=fir[c];t;t=nes[t])
        {
            if(!q[t] || tim[v[t]]==Time_Stamp) continue;
            tim[v[t]]=Time_Stamp;
            d[v[t]]=d[c]+1;
            dl[++r]=v[t];
            if(v[t]==T) return true;
        }
    }
    return false;
}
int dfs(int c,int flow)
{
    if(c==T || flow==0) return flow;
    int ans=0;
    for(int& t=cur[c];t;t=nes[t])
    {
        if(!q[t] || tim[v[t]]!=Time_Stamp || d[v[t]]!=d[c]+1) continue;
        int tmp=dfs(v[t],Min(flow,q[t]));
        q[t]-=tmp; q[t^1]+=tmp;
        flow-=tmp; ans+=tmp;
        if(!flow) break;
    }
    if(!ans) d[c]=-1;
    return ans;
}
int dinic()
{
    int ans=0;
    while(bfs())
    {
        for(int i=0;i<=T;i++) cur[i]=fir[i];
        ans+=dfs(S,INF);
    }
    return ans;
}
int main()
{
    using namespace IStream;
    cin>>n>>m;
    S=0;T=n+m+1;
    for(int i=1,x,y,cost;i<=n;i++)
    {
        x=get_int(); y=get_int();
        edge(S,i,x);
        ans+=x;
        while(y--)
        {
            x=get_int(); cost=get_int();
            edge(i,n+x,cost);
        }
    }
    for(int i=1,x;i<=m;i++)
    {
        x=get_int();
        edge(n+i,T,x);
    }
    ans-=dinic();
    printf("%d\n",ans);
    return 0;
}