【bzoj1061】[NOI2008]志愿者招募 线性规划与费用流
题目描述
输入
输出
样例输入
3 3
2 3 4
1 2 2
2 3 5
3 3 2
样例输出
14
题解
线性规划与费用流
“事实上任何最大流,最小费用流,上下界网络流都是在解决一个线性规划问题”(摘自 wyfcyx 大爷的ppt)
具体地,在一个网络流问题中,把边看作变量,把点看作限制条件,就会得到一个线性规划问题。
当然,由于网络流的特性,这类线性规划问题只能是:每个变量出现次数为且仅为2次,且1次系数为+1,1次系数为-1。
如果我们能够将这类特殊的线性规划问题转化为网络流问题,就可以圆满地解决原问题。
考虑本题,限制条件就是每天大的志愿者数大于等于规定数量,求 ∑费用*数量 的最小值。
这样可能不是很清楚,这里用样例举例子:
设$x_i$表示第i类志愿者的人数
限制条件:$\begin{cases}x_i\ge0\\x_1\ge2\\x_1+x_2\ge3\\x_2+x_3\ge4\end{cases}$,要最小化$2x_1+5x_2+2x_3$
我们把不等关系转化为相等关系,可以得到新的限制条件:$\begin{cases}x_i\ge0\\y_i\ge0\\x_1=2+y_1\\x_1+x_2=3+y_2\\x_2+x_3=4+y_3\end{cases}$
这样尽管多了几个变量,却把不等关系转化为容易处理的相等关系,易于建图。
但是这样依然不满足我们使用网络流解决线性规划问题的条件,因为变量出现次数不全为2,且系数不为+1和-1。
此时我们想到差分,把0=0分别放到限制条件中相等关系的最上端与最下端,然后上下进行差分并移项,可以得到:$\begin{cases}x_i\ge0\\y_i\ge0\\x_1-y_1-2=0\\x_2+y_1-y_2-1=0\\x_3-x_1+y_2-y_3-1=0\\-x_2-x_3+y_3+4=0\end{cases}$
这样就把限制条件“神奇”地转化为能够使用网络流解决的问题。为什么?因为用到了题目中的“志愿者工作时间是si天到ti天”,即一定是连续的。我们差分的本质是:$x_i$系数为+1时表示第i类志愿者刚开始工作,$x_i$系数为-1时表示第i类志愿者刚结束工作。所以一定是两个时间点。
我们再回过头来看这个线性规划问题。如何来使用网络流来解决它?这用到了网络流“流量守恒”的性质,即除S和T外,流入流量=流出流量。
我们把系数为+1看作要流出的流量,系数为-1看作要流入的流量,那么显然正负相等,符合条件。对于同一个变量的值是固定的,所以应从+1流向-1;对于常数项,系数为+1则从该点流向T,系数为-1则从S流向该点。
由于限制条件是与变量$x_i$有关的,所以在变量$x_i$对应的边上加上费用,然后跑费用流即可出解。
总结一下建图方法:把题目中的n天转化为n个等式,差分得到n+1个,代表图中的点;对于第i类志愿者,加边si->ti+1,容量为inf,费用为ci;加边j-1->j,容量为inf,费用为0,代表限制条件中的y。对于每天需要的志愿者数量$a_i(1\le i\le n+1)$,如果$a_i>a_{i-1}$,则加边S->i,容量为$a_i-a_{i-1}$,费用为0(因为差分后右面系数为正,移项后为负,代表流入);否咋加边i->T,容量为$a_{i-1}-a_i$,费用为0.
然后跑最小费用最大流出解。
说了这么多其实代码还是很简单的。
#include <cstdio>
#include <cstring>
#include <queue>
#define N 5000
#define M 100000
#define inf 0x3f3f3f3f
using namespace std;
queue<int> q;
int a[N] , head[N] , to[M] , val[M] , cost[M] , next[M] , cnt = 1 , s , t , dis[N] , from[N] , pre[N];
void add(int x , int y , int v , int c)
{
to[++cnt] = y , val[cnt] = v , cost[cnt] = c , next[cnt] = head[x] , head[x] = cnt;
to[++cnt] = x , val[cnt] = 0 , cost[cnt] = -c , next[cnt] = head[y] , head[y] = cnt;
}
bool spfa()
{
int x , i;
memset(from , -1 , sizeof(from));
memset(dis , 0x3f , sizeof(dis));
dis[s] = 0 , q.push(s);
while(!q.empty())
{
x = q.front() , q.pop();
for(i = head[x] ; i ; i = next[i])
if(val[i] && dis[to[i]] > dis[x] + cost[i])
dis[to[i]] = dis[x] + cost[i] , from[to[i]] = x , pre[to[i]] = i , q.push(to[i]);
}
return ~from[t];
}
int mincost()
{
int ans = 0 , i , k;
while(spfa())
{
k = inf;
for(i = t ; i != s ; i = from[i]) k = min(k , val[pre[i]]);
ans += k * dis[t];
for(i = t ; i != s ; i = from[i]) val[pre[i]] -= k , val[pre[i] ^ 1] += k;
}
return ans;
}
int main()
{
int n , m , i , x , y , z;
scanf("%d%d" , &n , &m) , s = 0 , t = n + 2;
for(i = 1 ; i <= n ; i ++ ) scanf("%d" , &a[i]) , add(i + 1 , i , inf , 0);
for(i = 1 ; i <= m ; i ++ ) scanf("%d%d%d" , &x , &y , &z) , add(x , y + 1 , inf , z);
for(i = 1 ; i <= n + 1 ; i ++ )
{
if(a[i] - a[i - 1] > 0) add(s , i , a[i] - a[i - 1] , 0);
if(a[i] - a[i - 1] < 0) add(i , t , a[i - 1] - a[i] , 0);
}
printf("%d\n" , mincost());
return 0;
}
【bzoj1061】[NOI2008]志愿者招募 线性规划与费用流的更多相关文章
- BZOJ1061 NOI2008 志愿者招募 线性规划、费用流
传送门 一道思路很妙的线性规划网络流 设\(X_i\)表示第\(i\)天需要的人数,\(P_i\)表示第\(i\)种人雇佣的个数 那么我们可以列出一系列式子 比如说样例就可以列出三个式子: \(P_1 ...
- [BZOJ1061][Noi2008]志愿者招募 线性规划+费用流
题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=1061 根据题意列方程,然后用网络流解线性规划. 题解直接贴ByVoid的吧,太神了:htt ...
- BZOJ1061: [Noi2008]志愿者招募(线性规划)
Time Limit: 20 Sec Memory Limit: 162 MBSubmit: 5725 Solved: 3437[Submit][Status][Discuss] Descript ...
- 网络流解线性规划问题 BZOJ1061: [Noi2008]志愿者招募
线性规划定义: 在给定有限的资源和竞争约束情况下,很多问题都可以表述为最大化或最小化某个目标.如果可以把目标指定为某些变量的线性函数,而且如果可以将资源约束指定为这些变量的等式或不等式,则得到了一个线 ...
- [BZOJ1061][Noi2008]志愿者招募
[BZOJ1061][Noi2008]志愿者招募 试题描述 申奥成功后,布布经过不懈努力,终于成为奥组委下属公司人力资源部门的主管.布布刚上任就遇到了一个难 题:为即将启动的奥运新项目招募一批短期志愿 ...
- 【费用流】BZOJ1061: [Noi2008]志愿者招募(这题超好)
1061: [Noi2008]志愿者招募 Time Limit: 20 Sec Memory Limit: 162 MBSubmit: 5291 Solved: 3173[Submit][Stat ...
- 线性规划费用流解法(Bzoj1061: [Noi2008]志愿者招募)
题面 传送门 Sol 线性规划费用流解法用与求解未知数为非负数的问题 这道题可以列出一堆形如 \(x[i]+x[j]+x[k]+...>=a[p]\) 的不等式 我们强行给每个式子减去一个东西, ...
- [BZOJ1061] [Noi2008] 志愿者招募 (费用流)
Description 申奥成功后,布布经过不懈努力,终于成为奥组委下属公司人力资源部门的主管.布布刚上任就遇到了一个难 题:为即将启动的奥运新项目招募一批短期志愿者.经过估算,这个项目需要N 天才能 ...
- BZOJ.1061.[NOI2008]志愿者招募(线性规划 对偶原理 单纯形 / 费用流SPFA)
题目链接 线性规划 用\(A_{ij}=0/1\)表示第\(i\)天\(j\)类志愿者能否被招募,\(x_i\)为\(i\)类志愿者招募了多少人,\(need_i\)表示第\(i\)天需要多少人,\( ...
随机推荐
- spring cloud 学习之路由网关(zuul)
学习自方志朋的博客 http://blog.csdn.net/forezp/article/details/69939114 在微服务架构中,需要几个基础的服务治理组件,包括服务注册与发现.服务消费. ...
- pyqt5通过qt designer 设计方式连接多个UI图形界面
当我们通过pyqt开发时,eric6为我们提供了一个方便的工具:图形化的绘制UI工具--qtdesigner.我们可以通过它开发多个UI,然后利用信号-槽工具,将功能代码附着在上面.也可以将多个界面连 ...
- 利用PHP str_replace()函数替换符合百度MIP内容标准
了解过百度MIP的同学都知道(什么?你不知道什么是百度MIP?就是移动网页加速器鸭
- CentOS下安装pip
CentOS下安装pip 通常情况下使用命令: yum -y install pip 也有可能报错,无法安装.这是应该使用第二种方法. 1.首先需要先安装扩展源EPEL: yum -y install ...
- hive 学习系列三(表格的创建create-table)
表格创建: 语法 第一种建表的形式: 说明: temporary 临时表,在当前回话内,这张表有效,当回话结束,可以理解为程序结束,则程序终止. external 外部表, hdfs 上的表的文件,并 ...
- ruby net/http模块使用
ruby中的NET::HTTP:这里暂时先列出几个固定用法: 其中一,二不支持请求头设置(header取ruby默认值),只能用于基本的请求,不支持持久连接,如果您执行许多HTTP请求,则不推荐它们: ...
- [Cracking the Coding Interview] 4.6 Successor 后继节点
Write an algorithm to find the 'next' node(i.e. in-order successor) of a given node in a binary sear ...
- Android面试收集录 蓝牙与WiFi
1.打开手机中的蓝牙功能有哪些方法? 法1:使用Intent ==>new Intent(BluetoothAdaper.ACTION_REQUEST_ENABLE); startActivi ...
- Linq中dbSet 的查询
1.Find:按照关键字的ID号来查询(速度快) 如: ADShiTi aDShiTi = db.ADShiTis.Find(id); 2.FirstOrDefault:根据部分条件查询,显示最前的一 ...
- Java:Random函数及其种子的作用
伪随机(preundorandom):通过算法产生的随机数都是伪随机!! 只有通过真实的随机事件产生的随机数才是真随机!!比如,通过机器的硬件噪声产生随机数.通过大气噪声产生随机数 Random生成的 ...