hdu 3879 hdu 3917 构造最大权闭合图 俩经典题

hdu3879  base station : 各一个无向图，点的权是负的，边的权是正的。自己建一个子图，使得获利最大。

一看，就感觉按最大密度子图的构想：选了边那么连接的俩端点必需选，于是就以边做点，轻轻松松构造了最大权闭合图。简单题。分分钟搞定。

hdu3917 :road  constructions :这题题目看了半天没理解。。。感觉描述的不好。。。一个有向图，每条路有响应公司承保，若选了该公司，那么该公司的路必需全部选，还有，该公司的承保的路的下面的一条路对应公司也要选，求最大获利。构图：开始在原图上瞎折腾，其实不然，用公司做结点，在原图上跑出关系，（若a公司选，则对应哪些公司 选），构造新图：最大权闭合之即可。。。。中午没睡醒。。迷迷糊糊的，因为1A的，数组开小了。。最讨厌用俩套链前来构图了。晕。。

  再重说一次最大权闭合图解法：源点向正权点连边对应点权值，原图边容量inf，负权点向汇点连边权值绝对值，ans=正权和-最大流。

#include<iostream> //hdu3879#include<queue>#include<cstdio>#include<cstring>#include<set>#include<vector>using namespace std;const int inf=0x3f3f3f3f;const int maxv=61000,maxe=300000;int nume=0;int head[maxv];int e[maxe][3];void inline adde(int i,int j,int c){    e[nume][0]=j;e[nume][1]=head[i];head[i]=nume;    e[nume++][2]=c;    e[nume][0]=i;e[nume][1]=head[j];head[j]=nume;    e[nume++][2]=0;}int ss,tt,n,m;int vis[maxv];int lev[maxv];bool bfs(){    for(int i=0;i<maxv;i++)      vis[i]=lev[i]=0;    queue<int>q;    q.push(ss);    vis[ss]=1;    while(!q.empty())    {        int cur=q.front();        q.pop();        for(int i=head[cur];i!=-1;i=e[i][1])        {            int v=e[i][0];            if(!vis[v]&&e[i][2]>0)            {                lev[v]=lev[cur]+1;                vis[v]=1;                q.push(v);            }        }    }    return vis[tt];}int dfs(int u,int minf){    if(u==tt||minf==0)return minf;    int sumf=0,f;    for(int i=head[u];i!=-1&&minf;i=e[i][1])    {        int v=e[i][0];        if(lev[v]==lev[u]+1&&e[i][2]>0)        {            f=dfs(v,minf<e[i][2]?minf:e[i][2]);            e[i][2]-=f;e[i^1][2]+=f;            sumf+=f;minf-=f;        }    }    if(!sumf) lev[u]=-1;    return sumf;}int dinic(){    int sum=0;    while(bfs())sum+=dfs(ss,inf);    return sum;};int cpy[maxv];int sumz=0;void read_build(){    for(int i=1;i<=n;i++)    {        scanf("%d",&cpy[i]);        adde(i,tt,cpy[i]);    }    int aa,bb,cc;    int numv=n+1;    for(int i=0;i<m;i++)     {        scanf("%d%d%d",&aa,&bb,&cc);        sumz+=cc;        adde(numv,aa,inf);        adde(numv,bb,inf);        adde(ss,numv++,cc);     }  /*  for(int i=1;i<=m+2;i++)      for(int j=head[i];j!=-1;j=e[j][1])      {           if(j%2==0)          printf("%d->%d:%d\n",i,e[j][0],e[j][2]);      }*/}void init(){    nume=0;sumz=0;    ss=n+m+1;tt=ss+1;    for(int i=0;i<maxv;i++)      {       head[i]=-1;      }}int main(){      while(~scanf("%d%d",&n,&m))     {         init();        read_build();       int ans;       ans=sumz-dinic();      printf("%d\n",ans);    }    return 0;}

#include<iostream>//3917#include<queue>#include<cstdio>#include<cstring>#include<set>#include<vector>using namespace std;const int inf=0x3f3f3f3f;const int maxv=5100,maxe=100000;int nume=0;int head[maxv];int e[maxe][3];void inline adde(int i,int j,int c){    e[nume][0]=j;e[nume][1]=head[i];head[i]=nume;    e[nume++][2]=c;    e[nume][0]=i;e[nume][1]=head[j];head[j]=nume;    e[nume++][2]=0;}int ss,tt,n,m,k;int vis[maxv];int lev[maxv];bool bfs(){    for(int i=0;i<maxv;i++)      vis[i]=lev[i]=0;    queue<int>q;    q.push(ss);    vis[ss]=1;    while(!q.empty())    {        int cur=q.front();        q.pop();        for(int i=head[cur];i!=-1;i=e[i][1])        {            int v=e[i][0];            if(!vis[v]&&e[i][2]>0)            {                lev[v]=lev[cur]+1;                vis[v]=1;                q.push(v);            }        }    }    return vis[tt];}int dfs(int u,int minf){    if(u==tt||minf==0)return minf;    int sumf=0,f;    for(int i=head[u];i!=-1&&minf;i=e[i][1])    {        int v=e[i][0];        if(lev[v]==lev[u]+1&&e[i][2]>0)        {            f=dfs(v,minf<e[i][2]?minf:e[i][2]);            e[i][2]-=f;e[i^1][2]+=f;            sumf+=f;minf-=f;        }    }    if(!sumf) lev[u]=-1;    return sumf;}int dinic(){    int sum=0;    while(bfs())sum+=dfs(ss,inf);    return sum;};int nume2=0;int head2[maxv];int e2[maxe][3];void inline adde2(int i,int j,int c){    e2[nume2][0]=j; e2[nume2][1]=head2[i]; head2[i]=nume2;    e2[nume2++][2]=c;}int cpy[maxv];int sumz=0;void read_build(){    for(int i=1;i<=m;i++)    {        scanf("%d",&cpy[i]);    }    scanf("%d",&k);    int aa,bb,cc,dd;    for(int i=0;i<k;i++)     {        scanf("%d%d%d%d",&aa,&bb,&cc,&dd);        cpy[cc]-=dd;        adde2(aa,bb,cc);     }    for(int i=1;i<=n;i++)      for(int j=head2[i];j!=-1;j=e2[j][1])      {          for(int vj=head2[e2[j][0]];vj!=-1;vj=e2[vj][1])          {              if(e2[vj][2]!=e2[j][2])                 adde(e2[j][2],e2[vj][2],inf);          }      }    for(int i=1;i<=m;i++)    {        if(cpy[i]>0)          {              adde(ss,i,cpy[i]);              sumz+=cpy[i];          }        else          adde(i,tt,-cpy[i]);    }  /*  for(int i=1;i<=m+2;i++)      for(int j=head[i];j!=-1;j=e[j][1])      {           if(j%2==0)          printf("%d->%d:%d\n",i,e[j][0],e[j][2]);      }*/}void init(){    nume=0;sumz=0;nume2=0;    ss=m+1;tt=ss+1;    for(int i=0;i<maxv;i++)      {         head2[i]=head[i]=-1;      }}int main(){      while(~scanf("%d%d",&n,&m)&&(n||m))     {         init();        read_build();       int ans;       ans=sumz-dinic();       if(ans<0)ans=0;      printf("%d\n",ans);    }    return 0;}