【BZOJ1061】[Noi2008]志愿者招募【单纯形法】
来源:互联网 发布:手机阿里云系统 编辑:程序博客网 时间:2024/05/23 19:53
双倍经验题,BZOJ3265。
先用对偶原则转换成求对偶问题的解,这样直接转化成了标准型,然后跑Simplex就好了。
下面是对样例的一个计算过程。
/* Footprints In The Blood Soaked Snow */#include <cstdio>#include <cmath>typedef double DB;typedef long long LL;const int maxn = 1005, maxm = 10005;const DB inf = 0x3f3f3f3f3f3f3f3f, eps = 1e-7;int n, m;DB b[maxm], c[maxn], cof[maxm][maxn], ans;inline void pivot(int id, int pos) {b[id] /= cof[id][pos];cof[id][pos] = 1 / cof[id][pos];for(int i = 1; i <= n; i++) if(i != pos) cof[id][i] *= cof[id][pos];for(int i = 1; i <= m; i++) if(i != id && fabs(cof[i][pos]) > eps) {b[i] -= cof[i][pos] * b[id];for(int j = 1; j <= n; j++) if(j != pos) cof[i][j] -= cof[i][pos] * cof[id][j];cof[i][pos] = -cof[i][pos] * cof[id][pos];}ans += c[pos] * b[id];for(int i = 1; i <= n; i++) if(i != pos) c[i] -= c[pos] * cof[id][i];c[pos] = -c[pos] * cof[id][pos];}inline DB simplex() {while(1) {int pos, id;for(pos = 1; pos <= n; pos++) if(c[pos] > eps) break;if(pos == n + 1) return ans;DB tmp = inf;for(int i = 1; i <= m; i++) if(cof[i][pos] > eps && b[i] / cof[i][pos] < tmp)tmp = b[i] / cof[i][pos], id = i;if(tmp == inf) return inf;pivot(id, pos);}}int main() {scanf("%d%d", &n, &m);for(int i = 1; i <= n; i++) scanf("%lf", &c[i]);for(int i = 1; i <= m; i++) {int x, y; scanf("%d%d", &x, &y);for(int j = x; j <= y; j++) cof[i][j] = 1;scanf("%lf", &b[i]);}printf("%lld\n", LL(simplex() + 0.5));return 0;}
0 0
- 【BZOJ1061】[Noi2008]志愿者招募【单纯形法】
- Bzoj1061:[Noi2008]志愿者招募:单纯形算法
- [单纯形+对偶] BZOJ1061: [Noi2008]志愿者招募
- [bzoj1061][线性规划][单纯形][Noi2008]志愿者招募]
- [BZOJ1061][NOI2008]志愿者招募 单纯形模板
- BZOJ1061: [Noi2008]志愿者招募 && 单纯形学习笔记
- bzoj1061 [Noi2008]志愿者招募 单纯形(模板)
- 【费用流|单纯形】BZOJ1061 [Noi2008]志愿者招募
- [BZOJ1061][Noi2008]志愿者招募
- bzoj1061【NOI2008】志愿者招募
- 【NOI2008】【BZOJ1061】志愿者招募
- bzoj1061: [Noi2008]志愿者招募
- 【NOI2008】BZOJ1061志愿者招募
- [bzoj1061][NOI2008]志愿者招募
- BZOJ1061: [Noi2008]志愿者招募
- BZOJ1061: [Noi2008]志愿者招募
- 【bzoj1061】[Noi2008]志愿者招募
- bzoj1061: [Noi2008]志愿者招募
- ios之启动图片适配
- leetcode之Repeated DNA Sequences
- Hibernate学习篇01——Hibernate基础
- Android事件分发机制完全解析,带你从源码的角度彻底理解(下)
- 设计TCP服务器的规则
- 【BZOJ1061】[Noi2008]志愿者招募【单纯形法】
- 12C ORA-错误汇总20
- new/delete和malloc/free的区别,sizeof的理解和struct内存对齐
- 从零开始学习EasyDarwin(概述篇)
- EM35x Development Kit开发环境搭建
- cookie 和session 的区别详解
- Servlet处理Json请求数据包
- cin.ignore()的使用
- Hibernate HQL之org.hibernate.QueryException:could not resolve