贴板子系列_1 网络流

网络流

//Dinic  
struct MAX_FLOW
{
    struct edge
    {
        int v,next,w;
    } vs[14000];
    int st[N],ee,h[N],cur[N];
    inline void addedge(int u,int v,int w)
    {
        vs[++ee].v=v;
        vs[ee].next=st[u];
        st[u]=ee;
        vs[ee].w=w;
    }
    inline void insert(int u,int v,int w)
    {
        addedge(u,v,w);
        addedge(v,u,0);
    }
    bool bfs()
    {
        memset(h,-1,sizeof(h));
        queue <int> q;
        q.push(S),h[S]=0;
        while(!q.empty())
        {
            int lx=q.front();
            q.pop();
            for(int i=st[lx]; i; i=vs[i].next)
                if(vs[i].w&&h[vs[i].v]==-1)
                {
                    h[vs[i].v]=h[lx]+1;
                    q.push(vs[i].v);
                }
        }
        return h[T]!=-1;
    }
    int dfs(int x,int f)
    {
        if(x==T) return f;
        int used=0,w;
        for(int i=cur[x]; i; i=vs[i].next)
            if(vs[i].w&&h[vs[i].v]==h[x]+1)
            {
                w=dfs(vs[i].v,min(f-used,vs[i].w));
                vs[i].w-=w;
                vs[i^1].w+=w;
                used+=w;
                if(vs[i].w) cur[x]=i;
                if(used==f) return used;
            }
        if(!used) h[x]=-1;
        return used;
    }
    void dinic()
    {
        while(bfs())
        {
            for(int i=S; i<=T; i++)
                cur[i]=st[i];
            ans+=dfs(S,inf);
        }
    }
}G;

费用流

struct edge
{
    int u,v,w;
    int next,c;
}vs[M];
int n,ee=1,st[N],from[N],S,T,m,tot;
int ans=0,vis[N],dis[N],p[N];
queue <int> q;
inline void addedge(int u,int v,int w,int c)
{
    vs[++ee].u=u;vs[ee].v=v;vs[ee].w=w;
    vs[ee].next=st[u];st[u]=ee;vs[ee].c=c;
}
inline void insert(int u,int v,int w,int c)
{
    addedge(u,v,w,c);addedge(v,u,0,-c);
}
bool spfa()
{
    for(int i=S;i<=T;i++) dis[i]=inf;
    while(!q.empty())q.pop();
    q.push(S);dis[S]=0;vis[S]=1;
    while(!q.empty())
    {
        int lx=q.front();q.pop();
        for(int i=st[lx];i;i=vs[i].next)
            if(vs[i].w&&dis[lx]+vs[i].c<dis[vs[i].v])
            {
                dis[vs[i].v]=dis[lx]+vs[i].c;
                from[vs[i].v]=lx;p[vs[i].v]=i;
                if(!vis[vs[i].v]) q.push(vs[i].v),vis[vs[i].v]=1;
            }
        vis[lx]=0;
    }
    return dis[T]!=inf;
}
inline void mcf()
{
    int ss=inf;
    for(int i=T;i!=S;i=from[i])
        ss=min(ss,vs[p[i]].w);
    for(int i=T;i!=S;i=from[i])
    {
        vs[p[i]].w-=ss;vs[p[i]^1].w+=ss;
        ans+=ss*vs[p[i]].c;
    }tot-=ss;
}