如果权值为\(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;
}

51nod 1392 装盒子(费用流)的更多相关文章

  1. 51nod 1392 装盒子

    有n个长方形盒子,第i个长度为Li,宽度为Wi,我们需要把他们套放.注意一个盒子只可以套入长和宽分别不小于它的盒子,并且一个盒子里最多只能直接装入另外一个盒子 (但是可以不断嵌套),例如1 * 1 可 ...

  2. 51nod 1392:装盒子 匈牙利+贪心

    1392 装盒子 基准时间限制:1 秒 空间限制:131072 KB 分值: 160 难度:6级算法题  收藏  关注 有n个长方形盒子,第i个长度为Li,宽度为Wi,我们需要把他们套放.注意一个盒子 ...

  3. 51nod 算法马拉松4 D装盒子(网络流 / 二分图最优匹配)

    装盒子   基准时间限制:1 秒 空间限制:131072 KB 分值: 160 有n个长方形盒子,第i个长度为Li,宽度为Wi,我们需要把他们套放.注意一个盒子只可以套入长和宽分别不小于它的盒子,并且 ...

  4. CSU 1948: 超级管理员(普通费用流&&zkw费用流)

    Description 长者对小明施加了膜法,使得小明每天起床就像马丁的早晨一样. 今天小明早上醒来发现自己成了一位仓管员.仓库可以被描述为一个n × m的网格,在每个网格上有几个箱子(可能没有).为 ...

  5. BZOJ2673 [Wf2011]Chips Challenge 费用流 zkw费用流 网络流

    https://darkbzoj.cf/problem/2673 有一个芯片,芯片上有N*N(1≤N≤40)个插槽,可以在里面装零件. 有些插槽不能装零件,有些插槽必须装零件,剩下的插槽随意. 要求装 ...

  6. NOIP 模拟 box - 费用流 / 匈牙利

    题目大意: 给出n(\(\le 200\))个盒子,第i个盒子长\(x_i\),宽\(y_i\),一个盒子可以放入长宽都大于等于它的盒子里,并且每个盒子里只能放入一个盒子(可以嵌套),嵌套的盒子的占地 ...

  7. hdu-5988 Coding Contest(费用流)

    题目链接: Coding Contest Time Limit: 2000/1000 MS (Java/Others)     Memory Limit: 65536/65536 K (Java/Ot ...

  8. POJ2195 Going Home[费用流|二分图最大权匹配]

    Going Home Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 22088   Accepted: 11155 Desc ...

  9. BZOJ3130: [Sdoi2013]费用流[最大流 实数二分]

    3130: [Sdoi2013]费用流 Time Limit: 10 Sec  Memory Limit: 128 MBSec  Special JudgeSubmit: 960  Solved: 5 ...

随机推荐

  1. 【BZOJ1125】【POI2008】poc - splay+哈希

    题意: Description n列火车,每条有l节车厢.每节车厢有一种颜色(用小写字母表示).有m次车厢交换操作.求:对于每列火车,在交换车厢的某个时刻,与其颜色完全相同的火车最多有多少. Inpu ...

  2. 自己对WEBGL坐标系的转换过程的理解【如图】

  3. js 操作table: insertRow(),deleteRow(),insertCell(),deleteCell()方法

    表格有几行: var trCnt = table.rows.length;  (table为Id ) 每行有几列:for (var i=0; i<trCnt; i++)              ...

  4. MVC架构之delegate

    Qt的MVC架构可以实现很多数据显示的功能,本次主要对代理进行一个总结: 重实现QStyledItemDelegate类,实现自定义类. (1)ComboxDelegate.h #ifndef COM ...

  5. 训练1-X

    输入n(n<100)个数,找出其中最小的数,将它与最前面的数交换后输出这些数. Input 输入数据有多组,每组占一行,每行的开始是一个整数n,表示这个测试实例的数值的个数,跟着就是n个整数.n ...

  6. Nginx 配置 Gzip 压缩

    打开配置文件 /etc/nginx/nginx.conf,取消掉以下的注释项: #gzip on; 取消后: gzip on; 在此配置后加上以下内容: gzip on; gzip_vary on; ...

  7. Centos文件查看命令字符

    文件(夹)查看类命令 ls--显示指定目录下内容 说明:ls 显示结果以不同的颜色来区分文件类别.蓝色代表目录,灰色代表普通文件,绿色代表可执行文件,红色代表压缩文件,浅蓝色代表链接文件. -a--- ...

  8. 使用c++Beep实现春节十二响蜂鸣程序

    直接编译运行即可 #include<bits/stdc++.h> #include<windows.h> using namespace std; char a[31][71] ...

  9. python 中的一些小命令

  10. 王立平--poser

    Poser是Metacreations公司推出的一款lemmaId=234814&ss_c=ssc.citiao.link" style="color:rgb(51,102 ...