zoukankan      html  css  js  c++  java
  • 51nod 1392 装盒子(费用流)

    如果权值为(1)就是最长反链。
    然而并不是。考虑用费用流。
    把每一个盒子(i)拆成i和(i+n)
    设源点为(S),汇点为(T)
    (S)向每一个i连容量为(1),费用为(L[i]*W[i])的边
    每一个(i)(T)连容量为(1),费用为(0)的边。
    每一个(i)(i+n)连容量为(1),费用为(0)的边。
    如果j能被i装下,就从i向j连一条容量为(1),费用为(-L[j]*W[j])的边。
    跑最小费用最大流即可。
    (PS:要特判用多个盒子完全一样的情况,在51nod提交时点运行代码给的样例自己画画就知道为什么了。)

    #include<iostream>
    #include<cstring>
    #include<cstdio>
    #include<cmath>
    #include<algorithm>
    #include<queue>
    using namespace std;
    #define int long long
    const int N=233;
    const int INF=1e9;
    int cnt,head[N*2];
    struct edge{
    	int to,nxt,flow,cost;
    }e[N*N*2+N*4];
    struct node{
    	int a,b;
    }c[N];
    bool cmp(node x,node y){
    	if(x.a==y.a)return x.b<y.b;
    	return x.a<y.a;
    }
    void add_edge(int u,int v,int flow,int cost){
    	cnt++;
    	e[cnt].nxt=head[u];
    	e[cnt].to=v;
    	e[cnt].flow=flow;
    	e[cnt].cost=cost;
    	head[u]=cnt;
    }
    int S,T,dis[N*2],vis[N*2],from[N*2],n,ans,a[N],b[N];
    bool spfa(){
    	for(int i=S;i<=T;i++)dis[i]=INF;
    	dis[S]=0;
    	queue<int> q;
    	q.push(S);vis[S]=1;
    	while(!q.empty()){
    		int u=q.front();
    		q.pop();
    		vis[u]=0;
    		for(int i=head[u];i;i=e[i].nxt){
    			int v=e[i].to;
    			if(dis[v]>dis[u]+e[i].cost&&e[i].flow){
    				dis[v]=dis[u]+e[i].cost;
    				from[v]=i;
    				if(vis[v]==0)vis[v]=1,q.push(v);
    			}
    		}
    	}
    	if(dis[T]==INF)return false;
    	int mn=INF;
    	for(int i=T;i!=S;i=e[from[i]^1].to){
    		mn=min(mn,e[from[i]].flow);
    	}
    	ans+=mn*dis[T];
    	for(int i=T;i!=S;i=e[from[i]^1].to){
    		e[from[i]].flow-=mn;
    		e[from[i]^1].flow+=mn;
    	}
    	return true;
    }
    int read(){
    	int sum=0,f=1;char ch=getchar();
    	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    	while(ch>='0'&&ch<='9'){sum=sum*10+ch-'0';ch=getchar();}
    	return sum*f;
    }
    signed main(){
    	n=read();
    	for(int i=1;i<=n;i++)c[i].a=read(),c[i].b=read();
    	sort(c+1,c+1+n,cmp);
    	int num=0;
    	for(int i=1;i<=n;i++)
    		if(c[i].a==c[i-1].a&&c[i].b==c[i-1].b)continue;
    		else a[++num]=c[i].a,b[num]=c[i].b,ans+=a[num]*b[num];
    	n=num;
    	cnt=1;
    	S=0;T=n*2+1;
    	for(int i=1;i<=n;i++){
    		add_edge(S,i,1,0),add_edge(i,S,0,0);
    		add_edge(i+n,T,1,0),add_edge(T,i+n,1,0);
    		for(int j=1;j<=n;j++){
    			if(i==j)continue;
    			if(a[j]<=a[i]&&b[j]<=b[i])
    				add_edge(i,j+n,1,-a[j]*b[j]),add_edge(j+n,i,0,a[j]*b[j]);
    		}
    	}
    	while(spfa());
    	printf("%lld",ans);
    	return 0;
    }
    
  • 相关阅读:
    正则表达式(四)--文本换行分割
    java加密类型和算法名称
    记事本记录日志
    DNS
    jstl--c:choose标签
    csv文本编辑引号问题
    JDBC----ReflectionUtils
    Jsp
    计算机网络 编程 总结:
    N颗骰子的问题
  • 原文地址:https://www.cnblogs.com/Xu-daxia/p/10436317.html
Copyright © 2011-2022 走看看