• bzoj 2055: 80人环游世界【有上下界有源汇最小费用最大流】


    连有上下界的边(ss,i,(0,m),0),(i',t,(0,m),0),表示从任意点开始和结束
    连(i,j,(0,m),d[i][j]),表示可以买票飞过去
    连(i,i',(v[i],v[i]),0),表示这条边会流固定v[i]个人
    连(s,ss,(m,m),0),表示限制m个人
    然后跑有上下界有源汇最小费用最大流(……)即可,新图的建立和最大流一样

    #include<iostream>
    #include<cstdio>
    #include<queue>
    using namespace std;
    const int N=305;
    int n,m,a[N],h[N],cnt=1,s,t,fr[N],dis[N],d[N],ans;
    bool v[N];
    struct qwe
    {
    	int ne,no,to,va,c;
    }e[N*N];
    int read()
    {
    	int r=0,f=1;
    	char p=getchar();
    	while(p>'9'||p<'0')
    	{
    		if(p=='-')
    			f=-1;
    		p=getchar();
    	}
    	while(p>='0'&&p<='9')
    	{
    		r=r*10+p-48;
    		p=getchar();
    	}
    	return r*f;
    }
    void add(int u,int v,int w,int c)
    {
    	cnt++;
    	e[cnt].ne=h[u];
    	e[cnt].no=u;
    	e[cnt].to=v;
    	e[cnt].va=w;
    	e[cnt].c=c;
    	h[u]=cnt;
    }
    void ins(int u,int v,int w,int c)
    {
    	add(u,v,w,c);
    	add(v,u,0,-c);
    }
    void wk(int u,int v,int l,int r,int c)
    {
    	d[v]+=l,d[u]-=l;
    	ans+=l*c;
    	ins(u,v,r-l,c);
    }
    bool spfa()
    {
    	queue<int>q;
    	for(int i=s;i<=t;i++)
    		dis[i]=1e9;
    	dis[s]=0;
    	v[s]=1;
    	q.push(s);
    	while(!q.empty())
    	{
    		int u=q.front();
    		q.pop();
    		v[u]=0;
    		for(int i=h[u];i;i=e[i].ne)
    			if(e[i].va>0&&dis[e[i].to]>dis[u]+e[i].c)
    			{
    				dis[e[i].to]=dis[u]+e[i].c;
    				fr[e[i].to]=i;
    				if(!v[e[i].to])
    				{
    					v[e[i].to]=1;
    					q.push(e[i].to);
    				}
    			}
    	}
    	return dis[t]<1e9;
    }
    void mcf()
    {
    	int x=1e9;
    	for(int i=fr[t];i;i=fr[e[i].no])
    		x=min(x,e[i].va);//cerr<<x<<endl;
    	for(int i=fr[t];i;i=fr[e[i].no])
    	{
    		e[i].va-=x;
    		e[i^1].va+=x;
    		ans+=x*e[i].c;
    	}
    }
    int main()
    {
    	n=read(),m=read();
    	for(int i=1;i<=n;i++)
    		a[i]=read(),wk(i,i+n,a[i],a[i],0);
    	for(int i=1;i<=n;i++)
    		for(int j=1;j<=n-i;j++)
    		{
    			int x=read();
    			if(x!=-1)
    				wk(i+n,i+j,0,m,x);
    		}
    	for(int i=1;i<=n;i++)
    		wk(n+n+1,i,0,m,0),wk(i+n,n+n+2,0,m,0);
    	wk(n+n+3,n+n+1,m,m,0);
    	wk(n+n+2,n+n+3,0,1e9,0);
    	s=0,t=n+n+4;
    	for(int i=1;i<=n+n+3;i++)
    	{
    		if(d[i]>0)
    			ins(s,i,d[i],0);
    		else
    			ins(i,t,-d[i],0);
    	}
    	while(spfa())
    		mcf();
    	printf("%d
    ",ans);
    	return 0;
    }
    
  • 相关阅读:
    vba create rectangle 矩形
    (zz)nnotation学习笔记2:关于Annotation与符号,还是要读arcmap帮助!!!
    backgroundworker
    C# , mongodb ,study
    工作中常用的dos命令 数据处理经验总结
    你需要明白的SQL SERVER书签查找(Bookmark Lookup)
    使用CalendarExtender
    js控制5秒后页面自动跳转
    HTML按秒数的页面跳转
    序列化与反序列化读取保存xml文件
  • 原文地址:https://www.cnblogs.com/lokiii/p/10779366.html
Copyright © 2020-2023  润新知