zoukankan      html  css  js  c++  java
  • bzoj1834[ZJOI2010]网络扩容

    传送门

    最小费用最大流,这么垃圾的题我想的居然挂了,正解就差一点点,我就假装自己秒了这个题吧,思路就是先跑一边最大流,然后在原图的残量网络中加一些边,就是对于原本的每条边((x,y))连一条容量为inf,费用为(cost_{(x,y)})的边,最后多建一个超级源点,向原源点连一条容量为k,费用为0的边就好了,sb题

    #include<cstdio>
    #include<iostream>
    #include<cstring>
    #include<algorithm>
    using namespace std;
    #define min(a,b) (a<b?a:b)
    #define rg register
    void read(int &x) {
    	char ch; bool ok;
    	for(ok=0,ch=getchar(); !isdigit(ch); ch=getchar()) if(ch=='-') ok=1;
    	for(x=0; isdigit(ch); x=x*10+ch-'0',ch=getchar()); if(ok) x=-x;
    }
    int ss,tt,tot,sum,n,m,k,id[200001],pre[200001],nxt[200001],h[100001],sla[100001],s,t,v[200001],dis[100001],w[200001],cnt=1,inf=1e9,ans,f[200001];bool vis[100001];
    void add(int x,int y,int z,int u)
    {
    	pre[++cnt]=y,nxt[cnt]=h[x],h[x]=cnt,v[cnt]=z,w[cnt]=u,id[cnt]=1;
    	pre[++cnt]=x,nxt[cnt]=h[y],h[y]=cnt,v[cnt]=0,w[cnt]=-u;
    }
    int dfs(int x,int flow)
    {
    	if(x==t){ans+=dis[t]*flow,sum+=flow;return flow;}
    	int f=flow;vis[x]=1;
    	for(rg int i=h[x];i;i=nxt[i])
    		if(!vis[pre[i]]&&v[i])
    		{
    			if(!(dis[x]+w[i]-dis[pre[i]]))
    			{
    				int y=dfs(pre[i],min(v[i],f));
    				f-=y,v[i]-=y,v[i^1]+=y;
    				if(!f)return flow;
    			}
    			else sla[pre[i]]=min(sla[pre[i]],dis[x]+w[i]-dis[pre[i]]);
    		}
    	return flow-f;
    }
    bool aug(int s,int t)
    {
    	int mn=inf;
    	for(rg int i=s;i<=t;i++)if(!vis[i])mn=min(mn,sla[i]),sla[i]=inf;
    	if(mn==inf)return 1;
    	for(rg int i=s;i<=t;i++)if(!vis[i])dis[i]+=mn;
    	return 0;
    }
    void zkw(int s,int t)
    {
    	memset(sla,63,sizeof sla);
    	while(1)
    	{
    		while(1){memset(vis,0,sizeof vis);if(!dfs(s,inf))break;}
    		if(aug(s,t))break;
    	}
    }
    int main()
    {
    	read(n),read(m),read(k);
    	for(rg int i=1,x,y,z,zz;i<=m;i++)read(x),read(y),read(z),add(x,y,z,0),read(zz),f[cnt^1]=zz;
    	s=1,t=n,zkw(s,t),printf("%d ",sum);
    	for(rg int i=1;i<=n;i++)
    		for(rg int j=h[i];j;j=nxt[j])if(id[j])add(i,pre[j],inf,f[j]);
    	ss=0,add(ss,s,k,0),zkw(ss,t),printf("%d
    ",ans);
    }
    
  • 相关阅读:
    .NET设计模式系列文章[转]
    asp.net MVC+NHibernate项目开发日志
    NHibernate 中关于 unsavedvalue 的问题
    领域驱动设计DDD之读书笔记
    让人郁闷的“DesktopCompatible”
    实战DDD(DomainDriven Design)领域驱动设计
    手机软件开发中的P/Invoke
    如何配置NHibernate的XML配置文件的智能感知(IntelliSense)功能 ?
    转]WinCE中的未公开函数
    Nhibernate+MVC开发日志
  • 原文地址:https://www.cnblogs.com/lcxer/p/10223451.html
Copyright © 2011-2022 走看看