hdu 3917 (최대 권 폐쇄 도)

제목: 정 부 는 일부 길, m 개 회사 가 수리 하고 모든 회사 가 도 로 를 수리 할 때 정부 에 세금 을 내야 한다. 도로 정 부 는 회사 비용 을 지불 하고 정부 가 얻 을 수 있 는 최대 이윤 을 구 해 야 한다. 만약 에 한 회 사 를 선택 하여 프로젝트 를 맡 으 면 이 회사 가 맡 은 다른 프로젝트 도 반드시 그 가 책임 지고 관련 이 있 는 회사 도 선택해 야 한다.
건 도: 이윤 = 모든 회사 의 세금 - (선택 하지 않 은 회사 의 세금 + 회사 의 보 조 를 선택) 후 두 가지 와 가장 작은 상황 에서 이윤 이 가장 크다. 외환 연결 회사 의 용량 은 회사 의 세금 이 고 회사 와 외환 연결 용량 은 이 회사 에 주 는 보조 이 며 회사 간 의존 관계 가 있 는 연결 용량 은 무한대 이다. 후 두 가 지 는 이 네트워크 의 최소 분할 이다.

#include<stdio.h>
#include<string.h>
const int N=5100;
const int inf=0x3fffffff;
int dis[N],gap[N],start,end,ans,head[N],num,n,m;
struct edge
{
    int st,ed,next,flow;
}e[N*50];
struct node
{
    int x,y,c,w;
}P[3010];
void addedge(int x,int y,int w)
{
    e[num].st=x;e[num].ed=y;e[num].flow=w;e[num].next=head[x];head[x]=num++;
    e[num].st=y;e[num].ed=x;e[num].flow=0;e[num].next=head[y];head[y]=num++;
}
int dfs(int u,int minflow)
{
    if(u==end)return minflow;
    int i,v,f,min_dis=ans-1,flow=0;
    for(i=head[u];i!=-1;i=e[i].next)
    {
        v=e[i].ed;
        if(e[i].flow<=0)continue;
        if(dis[v]+1==dis[u])
        {
            f=dfs(v,e[i].flow>minflow-flow?minflow-flow:e[i].flow);
            e[i].flow-=f;
            e[i^1].flow+=f;
            flow+=f;
            if(flow==minflow)break;
            if(dis[start]>=ans)return flow;
        }
        min_dis=min_dis>dis[v]?dis[v]:min_dis;
    }
    if(flow==0)
    {
        if(--gap[dis[u]]==0)
            dis[start]=ans;
        dis[u]=min_dis+1;
        gap[dis[u]]++;
    }
    return flow;
}
int isap()
{
    int maxflow=0;
    memset(gap,0,sizeof(gap));
    memset(dis,0,sizeof(dis));
    gap[0]=ans;
    while(dis[start]<ans)
        maxflow+=dfs(start,inf);
    return maxflow;
}
int main()
{
    int i,j,x,sum,k,cot[N];
    while(scanf("%d%d",&n,&m),n||m)
    {
        memset(head,-1,sizeof(head));
        num=0;start=0;end=m+1;ans=end+1;sum=0;
        memset(cot,0,sizeof(cot));
        for(i=1;i<=m;i++)
        {
            scanf("%d",&x);
            addedge(start,i,x);
            sum+=x;
        }
        scanf("%d",&k);
        for(i=0;i<k;i++)
        {
            scanf("%d%d%d%d",&P[i].x,&P[i].y,&P[i].c,&P[i].w);
            cot[P[i].c]+=P[i].w;
        }
        for(i=0;i<k;i++)
        {
            for(j=i+1;j<k;j++)
            {
                if(P[i].c==P[j].c)continue;
                if(P[j].y==P[i].x)addedge(P[j].c,P[i].c,inf);
                else if(P[i].y==P[j].x)addedge(P[i].c,P[j].c,inf);
            }
        }
        for(i=1;i<=m;i++)
            addedge(i,end,cot[i]);
        printf("%d
",sum-isap());
    }
    return 0;
}