bzoj 3112(zjoi 2013 防守战线) 线性规划+网络流
来源:互联网 发布:什么是数据挖掘 编辑:程序博客网 时间:2024/03/29 19:04
题目链接http://www.lydsy.com/JudgeOnline/problem.php?id=3112
话说这道题真心变态。。。终于见识到了浙江省选的难度了。。。
我们根据这道题的样例容易写出这样的不等式组:
x2+x3>=1
x1+x2+x3+x4+x5>=4
x3+x4+x5>=2
所求为min{x1+5*x2+6*x3+3*x4+4*x5},其中x[i]表示在i位置修建塔的个数。
首先打开脑洞,我们容易得出(容易你妹!)要用对偶图的性质来做。http://cxjyxx.me/?p=261该博客中有提到为什么要转成对偶图。那么现在问题就来了,怎样将一堆不等式转为另一堆不等式。本人开始看了很多题解都不懂,然后在看到这张图片时顿悟了= =||
根据这张表,我们可以把样例表示成
其中c为min{x1+5*x2+6*x3+3*x4+4*x5}的系数,对于矩形每一列,0或1的取值表示在原不等式中,每一项x[i]的系数。那么把矩形横着来看,每一排的0或1的取值表示在新不等式中,每一项的系数。根据此表,我们可以得到新的不等式组,其中y为新变量。
y2<=1
y1+y2<=5
y1+y2+y3<=6
y2+y3<=3
y2+y3<=4
所求即为max{y1+5*y2+2*y3}。
我们添加变量,将原式放缩为
y2+y4=1
y1+y2+y5=5
y1+y2+y3+y6=6
y2+y3+y7=3
y2+y3+y8=4
然后我们依次用上一个式子减去下一个式子,就可以得到
y2+y4=1
y1+y5-y4=4
y3+y6-y5=1
y7-y1-y6=-3
y8-y7=1
-y2-y3-y8=-4
然后我们仔细观查这一堆不等式,我们就(bu)可(neng)以得到每个变量在其中出现2次,且一正一负。然后通过我们惊人的联想能力(尼玛)就不难想到,网络流。因为网络流满足流入流量等于流出流量。将每个不等式看做一个点,等式右边的看做容量且费用为0,加入两个虚拟节点,一个为源点,另一个为汇点。对于每个变量x从它为正的式子向为负的式子连一条容量为正无穷,权值为题中读入的边,再跑一次最小费用流即可。
后面的代码是多亏了这位神犇的博客,表示感谢!这是链接 http://www.cnblogs.com/zig-zag/archive/2013/04/23/3036970.html,我觉得他讲的比我好。。。
#include<cstdio>#include<cstring>#include<iostream>using namespace std;#define inf (1e9)int N,M;int a[1010][10010],next[10010];void input(){scanf("%d%d",&N,&M);for(int i=1;i<=N;i++)scanf("%d",&a[i][0]);for(int i=1;i<=M;i++){int x,y;scanf("%d%d%d",&x,&y,&a[0][i]);for(int j=x;j<=y;j++){a[j][i]=1;}}}void change(int x,int y){int last=-1;for(int i=0;i<=M;i++){if(a[x][i]){next[i]=last;last=i;}}for(int i=0;i<=N;i++){if(a[i][y]==0||x==i)continue;for(int j=last;j!=-1;j=next[j]){if(j==y)continue;a[i][j]-=a[i][y]*a[x][j];}a[i][y]=-a[i][y];}}int ask(){while(1){int sai=0;for(int i=1;i<=M;i++){if(a[0][i]>0){sai=i;break;}}if(sai==0)return -a[0][0];int temp,min_num=inf;for(int i=1;i<=N;i++){if(a[i][sai]>0&&a[i][0]<min_num){temp=i; min_num=a[i][0];}}change(temp,sai);}}void solve(){int ans=ask();printf("%d",ans); }int main(){input();solve();return 0;}
- bzoj 3112(zjoi 2013 防守战线) 线性规划+网络流
- [BZOJ 3112][ZJOI 2013]防守战线
- 线性规划单纯型法(bzoj 3112: [Zjoi2013]防守战线)
- BZOJ 3112 [Zjoi2013]防守战线 线性规划
- bzoj 3112 从 线性规划 到 网络流 浙江oi2013 战线防守
- BZOJ 3112 防守战线
- BZOJ 3112 [Zjoi2013]防守战线
- [LP对偶费用流] BZOJ 3112 [Zjoi2013]防守战线
- bzoj 3112: [Zjoi2013]防守战线 (单纯形)
- BZOJ 3112 Zjoi2013 防守战线 单纯形
- bzoj 3112: [Zjoi2013]防守战线 单纯形
- bzoj 3112: [Zjoi2013]防守战线 单纯形
- [BZOJ 3112][Zjoi2013]防守战线:单纯形
- BZOJ 3112([Zjoi2013]防守战线-单纯形)
- [线性规划 费用流]BZOJ1061 志愿者招募 && BZOJ3112防守战线
- [单纯形 || 差分费用流 || 辅助变量费用流] BZOJ 3112 [Zjoi2013]防守战线
- BZOJ 3265 志愿者招募(BZOJ 1061)加强版 && BZOJ 3112 [Zjoi2013]防守战线 单纯形
- [BZOJ3112][ZJOI2013]防守战线(单纯形)
- Flume,Kafka与Storm整合
- Zookeeper API 备忘(持续更新)
- SVN
- Minitab 推出过程改善新的 Monte Carlo 仿真软件
- Android Activity和Intent机制学习笔记
- bzoj 3112(zjoi 2013 防守战线) 线性规划+网络流
- Android学习:onClick +findViewById+Toast基础
- Vim 配色(color scheme)与随机配色
- 配置Spring需要的jar包!
- Nginx安装部署
- diff 和 patch
- python数组的使用
- Google Hacking
- uva11235 RMQ