题解

可能是世界上最裸的一个单纯形

(话说全幺模矩阵是啥我到现在都不知道)

假装我们已经看过了算导,或者xxx的论文,知道了单纯形是怎么实现的

扔一个blog走掉。。https://www.cnblogs.com/ECJTUACM-873284962/p/7097864.html

那么我们根据题意可以列出这样的方程

\(x_i\)表示第\(i\)类志愿者招募的个数

根据题目可列线性规划的式子

以样例为例

\(z = min 2x_1 + 5x_2 + 2x_3\)

\(x_1 + 0 + 0 >= 2\)

\(x_1 + x_2 + 0 >= 3\)

\(0 + x_2 + x_3 >= 4\)

\(x_1,x_2,x_3 >= 0\)

显然,标准型要求我们这些式子是小于号并且z要取max

好吧,反号?

不过我们有个很神奇的原理叫对偶原理

(我是真的不知道为啥……)

也就是

\(min c^T X\)

\(Ax = b\)

等价于

\(max b^T X\)

\(A^TX = c\)

好吧,这是我们喜欢的形式啊

然后我们的方程就可以写成这个样子

\(z = max 2x_1 + 3x_2 + 4x_3\)

\(x_1 + x_2 + 0 <= 2\)

\(0 + x_2 + x_3 <= 5\)

\(0 + 0 + x_3 <= 2\)

我们把这个方程转换成松弛型(也就是全是等于号)

\(x_1 + x_2 + 0 + x_4 = 2\)

\(0 + x_2 + x_3 + x_5 = 5\)

\(0 + 0 + x_3 + x_6 = 2\)

我们3个方程组,6个元,是消不出来的,我们这是一些作为基变量,剩下的非基变量都设成0,那样的话一定是单纯形上的一个顶点

这是一个\(m * (n + m)\)的矩阵,有点大

我们在处理的时候,初始设定所有的\(x_4,x_5,x_6\)作为基变量

每一个方程就是一个关于基变量的等式,我们找到一个替入变量,找到能使替入变量值最大的方程组的等式,将替入变量的位置当做替出变量的位置,矩阵就是\(nm\)的了

代码

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <queue>
#include <cmath>
#define enter putchar('\n')
#define mp make_pair
#define pb push_back
#define fi first
#define se second
#define eps 1e-7
//#define ivorysi
using namespace std;
typedef long long int64;
typedef double db;
template<class T>
void read(T &res) {
res = 0;char c = getchar();T f = 1;
while(c < '0' || c > '9') {
if(c == '-') f = -1;
c = getchar();
}
while(c >= '0' && c <= '9') {
res = res * 10 + c - '0';
c = getchar();
}
res *= f;
}
template<class T>
void out(T x) {
if(x < 0) {putchar('-');x = -x;}
if(x >= 10) {
out(x / 10);
}
putchar('0' + x % 10);
}
const int MAXN = 1005;
const int MAXM = 10005;
int N,M;
db b[MAXM],c[MAXN],a[MAXM][MAXN],ans;
void pivot(int id,int p) {
b[id] /= a[id][p];
for(int j = 1 ; j <= N ; ++j) {
if(j != p) a[id][j] /= a[id][p];
}
a[id][p] = 1 / a[id][p];
for(int i = 1 ; i <= M ; ++i) {
if(i != id) {
if(fabs(a[i][p]) < eps) continue;
for(int j = 1 ; j <= N ; ++j) {
if(j != p) a[i][j] -= a[i][p] * a[id][j];
}
b[i] -= a[i][p] * b[id];a[i][p] = -a[i][p] * a[id][p];
}
}
for(int j = 1 ; j <= N ; ++j) {
if(j != p) c[j] -= c[p] * a[id][j];
}
ans += c[p] * b[id];c[p] = -c[p] * a[id][p];
} void Init() {
read(N);read(M);
for(int i = 1 ; i <= N ; ++i) scanf("%lf",&c[i]);
int s,t;
for(int i = 1 ; i <= M ; ++i) {
scanf("%d %d %lf",&s,&t,&b[i]);
for(int j = s ; j <= t ; ++j) {
a[i][j] = 1.0;
}
}
}
void Solve() {
while(1) {
db t = -1;int p = 0;
for(int i = 1 ; i <= N ; ++i) {
if(c[i] > t) {
t = c[i];p = i;
}
}
if(t <= 0) {printf("%.0lf\n",ans);return;}
t = 0x5fffffff;int id = 0;
for(int j = 1 ;j <= M ; ++j) {
if(a[j][p] > 0 && b[j] / a[j][p] < t) {
t = b[j] / a[j][p];
id = j;
}
}
pivot(id,p);
}
}
int main() {
#ifdef ivorysi
freopen("f1.in","r",stdin);
#endif
Init();
Solve();
}

【BZOJ】1061: [Noi2008]志愿者招募的更多相关文章

  1. BZOJ 1061: [Noi2008]志愿者招募

    1061: [Noi2008]志愿者招募 Time Limit: 20 Sec  Memory Limit: 162 MBSubmit: 4064  Solved: 2476[Submit][Stat ...

  2. BZOJ 1061: [Noi2008]志愿者招募 [单纯形法]【学习笔记】

    1061: [Noi2008]志愿者招募 Time Limit: 20 Sec  Memory Limit: 162 MBSubmit: 3975  Solved: 2421[Submit][Stat ...

  3. BZOJ 1061: [Noi2008]志愿者招募 费用流

    1061: [Noi2008]志愿者招募 题目连接: http://www.lydsy.com/JudgeOnline/problem.php?id=1061 Description 申奥成功后,布布 ...

  4. BZOJ 1061: [Noi2008]志愿者招募【单纯形裸题】

    1061: [Noi2008]志愿者招募 Time Limit: 20 Sec  Memory Limit: 162 MBSubmit: 4813  Solved: 2877[Submit][Stat ...

  5. BZOJ 1061: [Noi2008]志愿者招募 [单纯形法]【学习笔记看另一篇吧】

    1061: [Noi2008]志愿者招募 Time Limit: 20 Sec  Memory Limit: 162 MBSubmit: 3975  Solved: 2421[Submit][Stat ...

  6. BZOJ.1061.[NOI2008]志愿者招募(线性规划 对偶原理 单纯形 / 费用流SPFA)

    题目链接 线性规划 用\(A_{ij}=0/1\)表示第\(i\)天\(j\)类志愿者能否被招募,\(x_i\)为\(i\)类志愿者招募了多少人,\(need_i\)表示第\(i\)天需要多少人,\( ...

  7. BZOJ 1061: [Noi2008]志愿者招募(线性规划与网络流)

    http://www.lydsy.com/JudgeOnline/problem.php?id=1061 题意: 思路: 直接放上大神的建模过程!!!(https://www.byvoid.com/z ...

  8. 【刷题】BZOJ 1061 [Noi2008]志愿者招募

    Description 申奥成功后,布布经过不懈努力,终于成为奥组委下属公司人力资源部门的主管.布布刚上任就遇到了一个难题:为即将启动的奥运新项目招募一批短期志愿者.经过估算,这个项目需要N 天才能完 ...

  9. BZOJ 1061 [Noi2008]志愿者招募(费用流)

    题目描述 申奥成功后,布布经过不懈努力,终于成为奥组委下属公司人力资源部门的主管.布布刚上任就遇到了一个难题:为即将启动的奥运新项目招募一批短期志愿者.经过估算,这个项目需要N 天才能完成,其中第i ...

  10. bzoj 1061 [Noi2008]志愿者招募(数学模型,MCMF)

    [题目链接] http://www.lydsy.com/JudgeOnline/problem.php?id=1061 [题意] 雇人满足每天至少需要的人数. [思路一] Byvoid的题解 clic ...

随机推荐

  1. day2 程序流程控制

    流程控制:1.调用方法.调用方法将导致控制流程离开当前方法,转移到被调用的方法 2.选择.java中有两种做出选择的机制:if/else语句和switch语句.三目运算符可以看作是if/else的一个 ...

  2. 关于构造IOCTL命令的学习心得

    在编写ioctl代码之前,需要选择对应不同命令的编号.为了防止对错误的设备使用正确的命令,命令号应该在系统范围内唯一,这种错误匹配并不是不会发生,程序可能发现自己正在试图对FIFO和audio等这类非 ...

  3. Windows下php,mysql,apache相关安装与配置,完善中…

    PHP 的安装 由于php是一个zip文件(非install版),安装较为简单解压就行.把解压的 php5.2.1-Win32重命名为 php5.并复制到安装盘目录下.例如安装路径为 c:\php5 ...

  4. bzoj千题计划133:bzoj3130: [Sdoi2013]费用流

    http://www.lydsy.com/JudgeOnline/problem.php?id=3130 第一问就是个最大流 第二问: Bob希望总费用尽量大,那肯定是把所有的花费加到流量最大的那一条 ...

  5. 贪心问题:区间覆盖 NYOJ 喷水装置(二)

    喷水装置(二) 时间限制:3000 ms  |  内存限制:65535 KB 难度:4 描述 有一块草坪,横向长w,纵向长为h,在它的橫向中心线上不同位置处装有n(n<=10000)个点状的喷水 ...

  6. angularJS $http $q $promise

    一天早晨,爹对儿子说:“宝儿,出去看看天气如何!” 每个星期天的早晨,爹都叫小宝拿着超级望远镜去家附近最高的山头上看看天气走势如何,小宝说没问题,我们可以认为小宝在离开家的时候给了他爹一个promis ...

  7. Java入门系列(三)面向对象三大特性之封装、继承、多态

    面向对象综述 封装 封装的意义,在于明确标识出允许外部使用的所有成员函数和数据项,或者叫接口. 有了封装,就可以明确区分内外,使得类实现者可以修改封装内的东西而不影响外部调用者:而外部调用者也可以知道 ...

  8. Java写的数据库连接池

    原文地址: http://lgscofield.iteye.com/blog/1820521 import java.sql.*; import java.util.Enumeration; impo ...

  9. JS日历,可获得指定日期周数及星期几

    需求来自一个朋友:编写一个简易日历.在文本框中输入要查找的日期,程序可以计算出这一天处在该年份的第几周,并且能判断出这一天到底是星期几. 应为要有交互,选择了Js来实现,也算是 结对编程 的初试吧. ...

  10. [hadoop]hadoop api 新版本与旧版本的差别

    突然现在对以后的职业方向有些迷茫,不知道去干什么,现在有一些语言基础,相对而言好的一些有Java和C,选来选去不知道该选择哪个方向,爬了好多网页后,觉得自己应该从java开始出发,之前有点心不在焉,不 ...