HDU

来源:互联网 发布:pop服务器端口 编辑:程序博客网 时间:2024/06/04 21:18

2017青岛网络赛的题。。当时刚学了连通图,比赛的时候看到了割边,怎么想怎么是连通图,明明是最小割= =?题目就写着最小割??好吧我还是知道的太少了。。

现在学网络流看到了最小割的题目,从HDU的problem里面从后往前翻了好久,终于找到了。。

最小割最少割边的裸题。。

建边的时候每条边的边权w=*(m+1)+1

最大流maxflow*/(m+1) 最少割边数maxflow%(m+1)

乘(m+1)是为了保证边数叠加后依然是余数,不至于影响最小割的结果。

只要模上一个比m大的数就可以。。。。


#include<bits/stdc++.h>using namespace std;const int N=234;const int M=500010;const int INF=0x3f3f3f3f;int n,m;struct node{    int u,v,next,cap;}eage[M];int source,sink;int cur[N],dep[N],gap[N];int head[N];int top;int S[N];void Add(int u,int v,int w){    eage[top].u=u;    eage[top].v=v;    eage[top].cap=w;    eage[top].next=head[u];    head[u]=top++;    eage[top].u=v;    eage[top].v=u;    eage[top].cap=0;    eage[top].next=head[v];    head[v]=top++;}void BFS(){    queue<int>q;    memset(dep,-1,sizeof(dep));    memset(gap,0,sizeof(gap));    gap[0]=1;    dep[sink]=0;    q.push(sink);    while(!q.empty())    {        int u=q.front();        q.pop();        for(int i=head[u];i!=-1;i=eage[i].next)        {            int v=eage[i].v;            if(dep[v]==-1)            {                q.push(v);                dep[v]=dep[u]+1;                gap[dep[v]]++;            }        }    }}int Sap(){    BFS();    memcpy(cur,head,sizeof(head));    int tot=0;    int u=source;    int ans=0;    while(dep[source]<n)    {        if(u==sink)        {            int Min=INF;            int inser;            for(int i=0;i<=tot-1;i++)            {                if(Min>eage[S[i]].cap)                {                    Min=eage[S[i]].cap;                    inser=i;                }            }            for(int i=0;i<=tot-1;i++)            {                eage[S[i]].cap-=Min;                eage[S[i]^1].cap+=Min;            }            ans+=Min;            tot=inser;            u=eage[S[tot]].u;        }        if(u!=sink&&gap[dep[u]-1]==0)          break;        int v;        int i;        for(i=cur[u];i!=-1;i=eage[i].next)        {            v=eage[i].v;            if(eage[i].cap&&dep[v]+1==dep[u])            {                break;            }        }        if(i!=-1)        {            cur[u]=i;            S[tot++]=i;            u=v;            continue;        }        int Min=n;        for(int i=head[u];i!=-1;i=eage[i].next)        {            int v=eage[i].v;            if(eage[i].cap&&dep[v]<Min)            {                Min=dep[v];                cur[u]=i;            }        }        gap[dep[u]]--;        dep[u]=Min+1;        gap[dep[u]]++;        if(u!=source)u=eage[S[--tot]].u;    }    return ans;}int main(){    int t;    scanf("%d",&t);    while(t--)    {        memset(head,-1,sizeof(head));        top=0;        scanf("%d%d",&n,&m);        scanf("%d%d",&source,&sink);        for(int i=1;i<=m;i++)        {            int u,v,w;            scanf("%d%d%d",&u,&v,&w);            Add(u,v,w*(m+1)+1);        }        int ans=Sap();        printf("%d\n",ans%(m+1));    }    return 0;}