首页 > 代码库 > hdu 2987最大权闭合图模板类型题
hdu 2987最大权闭合图模板类型题
/* 最大权闭合图模板类型的题,考验对知识概念的理解。 题意:现在要辞退一部分员工,辞退每一个员工可以的到一部分利益(可以是负的),并且辞退员工,必须辞退他的下属,求最大利益和辞退的最小人数。 最大权闭合图模板类型。 求出最大权后沿着源点s,dfs到的点就为最小的人数。 证明/* 转载:利用一个经典的trick:多关键字 > 建图前,对所有b[i],执行变换b[i]=b[i]*10000-1,然后,会惊异地发现, > 此时最大流所对应的方案就是满足辞退最少人数的了。 > 为什么?显然,变换后的流量r2除以10000后再取整就等于原来的流量,但是 > r2的后四位却蕴含了辞退人数的信息:每多辞退一个人,流量就会少1。 > > 剩下的就是如何根据一个网络流输出方案。 > 我的做法:从源点开始沿着残余网络dfs(只走没有满载的边), > 能dfs到的点对应的人就是需要辞退的。*/ */ #include"stdio.h" #include"string.h" #include"string" #include"queue" #define ll __int64 #define N 9999 #define inf 0x3fffffff using namespace std; struct node { ll u,v,w,next; }bian[N*50]; ll fee[N]; ll head[N],yong,s,t,dis[N]; void init(){ yong=0; memset(head,-1,sizeof(head)); memset(dis,-1,sizeof(dis)); } void addedge(ll u,ll v,ll w) { bian[yong].u=u; bian[yong].v=v; bian[yong].w=w; bian[yong].next=head[u]; head[u]=yong++; } void add(ll u,ll v,ll w) { addedge(u,v,w); addedge(v,u,0); } void bfs() { ll u,v,i; queue<ll>q; q.push(t); dis[t]=0; while(!q.empty()) { u=q.front(); q.pop(); for(i=head[u];i!=-1;i=bian[i].next) { v=bian[i].v; if(dis[v]==-1) { dis[v]=dis[u]+1; q.push(v); } } } return ; } ll ISAP() { ll sum=0; bfs(); ll gap[N],cur[N],stac[N],top,i; memset(gap,0,sizeof(gap)); for(i=s;i<=t;i++) { gap[dis[i]]++; cur[i]=head[i]; } ll k=s; top=0; while(dis[s]<t+1) { if(k==t) { ll minn=inf,index; for(i=0;i<top;i++) { ll e=stac[i]; if(minn>bian[e].w) { minn=bian[e].w; index=i; } } for(i=0;i<top;i++) { ll e=stac[i]; bian[e].w-=minn; bian[e^1].w+=minn; } sum+=minn; top=index; k=bian[stac[top]].u; } for(i=cur[k];i!=-1;i=bian[i].next) { ll v=bian[i].v; if(bian[i].w&&dis[k]==dis[v]+1) { cur[k]=i; k=v; stac[top++]=i; break; } } if(i==-1) { ll m=t+1; for(i=head[k];i!=-1;i=bian[i].next) if(m>dis[bian[i].v]&&bian[i].w) { m=dis[bian[i].v]; cur[k]=i; } if(--gap[dis[k]]==0)break; gap[dis[k]=m+1]++; if(k!=s) k=bian[stac[--top]].u; } } return sum; } ll cou=0,vis[N]; void dfs(ll u){//,ll pre) {会形成环越界 ll i,j; cou++; vis[u]=1; for(i=head[u];i!=-1;i=bian[i].next) { j=bian[i].v; if(bian[i].w&&!vis[j]) { dfs(j); } } return ; } int main() { ll n,m,i,j,sum; while(scanf("%I64d%I64d",&n,&m)!=EOF) { init(); s=0;t=n+1;sum=0; for(i=1;i<=n;i++) { scanf("%I64d",&fee[i]); if(fee[i]>0) { add(s,i,fee[i]); sum+=fee[i]; } else add(i,t,-fee[i]); } while(m--) { scanf("%I64d%I64d",&i,&j); add(i,j,inf); } cou=0; sum=sum-ISAP(); memset(vis,0,sizeof(vis)); dfs(0); printf("%I64d %I64d\n",cou-1,sum); } return 0; }
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。