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; }
hdu 3879 hdu 3917 构造最大权闭合图 俩经典题,布布扣,bubuko.com
时间: 2024-10-07 06:38:48