网络流模板

dinic模板：

 

struct edge{
    int from,to,w,next;
}e[N<<1];
int cnt;

int dep[1111],head[2111];
//queue<int> que;


void init(){
    memset(head,-1,sizeof(head));
    cnt=0;
}
void add(int u,int v,int w)
{
    e[cnt].from=u;
    e[cnt].to=v;
    e[cnt].w=w;
    e[cnt].next=head[u];
    head[u]=cnt++;

     e[cnt].from=v;
    e[cnt].to=u;
    e[cnt].w=0;
    e[cnt].next=head[v];
    head[v]=cnt++;
}

int bfs(int s,int t)
{
//    while(!que.empty()) que.pop();
    queue<int> que;
    memset(dep,-1,sizeof(dep));

    dep[s]=0;
    que.push(s);

    while(!que.empty()){
        int u=que.front();
        que.pop();
//        if(u==t) return 1;
            for(int i=head[u];i!=-1;i=e[i].next){
                int w=e[i].w;
                int v=e[i].to;
            if(w&&dep[v]==-1){
                dep[v]=dep[u]+1;
                que.push(v);
                if(v==t) return 1;
            }
        }

    }
    return 0;
}

int dfs(int u,int mi,int t){
    if(u==t||mi==0) return mi;
    int tmp;
    int ret=0;
//    for(int v=1;v<=t;v++){
        for(int i=head[u];i!=-1;i=e[i].next){
            int w=e[i].w;
            int v=e[i].to;
        if(w&&dep[v]==dep[u]+1&&mi-ret>0){ ///mi-ret>0少一次递归，优化
           tmp=dfs(v,min(mi-ret,w),t) ;
           if(!tmp) dep[v]=0;
           e[i].w-=tmp;
            e[i^1].w+=tmp;
            ret+=tmp;
            if(ret==mi) return ret;
        }
    }
    return ret;
}

int dinic(int s,int t)
{
    int ans=0;
    while(bfs(s,t)){
            ans+=dfs(s,INF,t);
    }
    return ans;
}

 

最小费用最大流模板：

 

const int N=100;
const int INF=0x3f3f3f3f;

struct edge{
    int u,v,flow,cost,next;
    edge(){}
    edge(int _u,int _v,int _flow,int _cost,int _next){
        u=_u;v=_v;flow=_flow;cost=_cost;next=_next;
    }
}e[30000];

int head[N],cnt;
queue<int> que;

void init()
{
    memset(head,-1,sizeof(head));
    cnt=0;
}

void add(int x,int y,int z,int c)
{
    e[cnt]=edge(x,y,z,c,head[x]);
    head[x]=cnt++;
    e[cnt]=edge(y,x,0,-c,head[y]);
    head[y]=cnt++;
}

bool vis[N];
int dis[N],pre[N];

bool spfa(int s,int t)
{
    memset(dis,INF,sizeof(dis));
    memset(vis,0,sizeof(vis));
    memset(pre,-1,sizeof(pre));

    dis[s]=0;
    vis[s]=true;
    while(!que.empty()) que.pop();

    que.push(s);

    while(!que.empty())
    {
        int x=que.front();
        que.pop();
        vis[x]=false;

        for(int i=head[x];i!=-1;i=e[i].next){
            if(e[i].flow){
                int y=e[i].v;
                if(dis[x]+e[i].cost<dis[y]){
                    dis[y]=dis[x]+e[i].cost;
                    pre[y]=i;
                    if(!vis[y]){
                        que.push(y);
                        vis[y]=true;
                    }
                }
            }
        }

    }
     return pre[t]!=-1;

}

int costflow(int s,int t,int x)///源点x，汇点t，需要的总流x
{
    int ret=0;
    int max_flow=0;

    while(spfa(s,t))
    {
        int min_flow=INF;
        for(int i=pre[t];i!=-1;i=pre[e[i].u]){
            min_flow=min(e[i].flow,min_flow);
        }
        for(int i=pre[t];i!=-1;i=pre[e[i].u]){
            e[i].flow-=min_flow;
            e[i^1].flow+=min_flow;
        }
        ret+=dis[t]*min_flow;
        max_flow+=min_flow;
    }

    if(max_flow<x) return -1;///达不到满流
    return ret;
}