UOJ.179.线性规划(单纯形)
这写得还不错:http://www.cnblogs.com/zzqsblog/p/5457091.html
引入基变量\(x_{i+n}\),将约束\(\sum_{i=1}^m a_{ij}x_j\leq b_i\)改写为$$x_{i+n}=b_i-\sum_{i=1}^m a_{ij}x_j$$。
目标函数为\(\sum_{i=1}^n C_ix_i\)。当存在\(r,c\)满足\(C_c>0\),\(B_r>0\),\(a_{rc}>0\),对第\(r\)个限制中的\(x_c\)做代换,即$$x_c=B_r-\sum_{j!=c}a_{rj}x_j-x_{r+n}$$(\(x_c\)成为基变量,\(x_{r+n}\)成为非基变量),然后代入目标函数中,非基变量取0,就一定可以使目标函数增大。这一步通过\(Pivot(r,c)\)(转轴)实现,同时要把其它约束中的\(x_c\)替换掉。
当所有\(B_r\geq 0\)时,所有非基变量取0可以得到一个基本解(零解),即一定存在解。若存在\(B_r<0\),在限制\(r\)中找一个\(a_{rc}<0\)的\(x_c\)做代换,就可以使\(B_r>0\)。
当然前提是任意\(x_i>0,i\in [1,n+m]\)。
//0ms 520kb
#include <cmath>
#include <cstdio>
#include <cctype>
#include <algorithm>
#define gc() getchar()
#define eps 1e-8
const int N=25;
const double INF=1e9;
int n,m,id[50];
double A[N][N],Ans[N];
inline int read()
{
int now=0,f=1;register char c=gc();
for(;!isdigit(c);c=='-'&&(f=-1),c=gc());
for(;isdigit(c);now=now*10+c-'0',c=gc());
return now*f;
}
void Pivot(int r,int c)//r:Basic varivle c:Nonbasic variable
{//交换基变量与非基变量
std::swap(id[r+n],id[c]);
double t=A[r][c]; A[r][c]=1;//
for(int i=0; i<=n; ++i) A[r][i]/=t;
for(int i=0; i<=m; ++i)//在其它等式中换掉基变量
if(fabs(A[i][c])>eps && i!=r)
{
t=A[i][c]; A[i][c]=0;//
for(int j=0; j<=n; ++j) A[i][j]-=t*A[r][j];
}
}
bool Init()
{
for(int r,c; ; )
{
r=c=0;
for(int i=1; i<=m; ++i)//B[r]<0
if(A[i][0]<-eps && (!r || rand()&1)) r=i;
if(!r) return 1;
for(int i=1; i<=n; ++i)//A[r][c]<0
if(A[r][i]<-eps && (!c || rand()&1)) c=i;
if(!c) return puts("Infeasible"),0;
Pivot(r,c);
}
}
bool Simplex()
{
for(int r,c; ; )
{
r=c=0;
for(int i=1; i<=n; ++i)//C[c]>0
if(A[0][i]>eps) {c=i; break;}
if(!c) return 1;
double mn=INF;//找一个系数为正且约束最紧的A[r][c]
for(int i=1; i<=m; ++i)
if(A[i][c]>eps && A[i][0]/A[i][c]<mn) r=i, mn=A[i][0]/A[i][c];
if(!r) return puts("Unbounded"),0;//无约束
Pivot(r,c);
}
}
int main()//x[i+n]=B[i]-∑a[i][j]*x[j]
{
srand(20180724);
n=read(), m=read(); int type=read();
for(int i=1; i<=n; ++i) A[0][i]=read();//目标函数系数C[i]
for(int i=1; i<=m; ++i)
{
for(int j=1; j<=n; ++j) A[i][j]=read();
A[i][0]=read();//B[i]
}
for(int i=1; i<=n; ++i) id[i]=i;
if(Init() && Simplex())
{
printf("%.8lf\n",-A[0][0]);//代换的时候Bi系数是负的s
if(type)
{
for(int i=1; i<=m; ++i) Ans[id[i+n]]=A[i][0];//成为基变量的xi取值即为bi,非基变量上的xi取0.
for(int i=1; i<=n; ++i) printf("%.8lf ",Ans[i]);
}
}
return 0;
}
UOJ.179.线性规划(单纯形)的更多相关文章
- UOJ#179. 线性规划[模板]
传送门 http://uoj.ac/problem/179 震惊,博主竟然还不会线性规划! 单纯形实在学不会啊……背个板子当黑盒用…… 学(chao)了NanoApe dalao的板子 #includ ...
- UOJ#179. 线性规划(线性规划)
描述 提交 自定义测试 这是一道模板题. (这个题现在标程挂了..哪位哥哥愿意提供一下靠谱的标程呀?) 本题中你需要求解一个标准型线性规划: 有 nn 个实数变量 x1,x2,…,xnx1,x2,…, ...
- uoj#179 线性规划
这是一道模板题. 本题中你需要求解一个标准型线性规划: 有nn个实数变量x1,x2,⋯,xnx1,x2,⋯,xn和mm条约束,其中第ii条约束形如∑nj=1aijxj≤bi∑j=1naijxj≤bi. ...
- 【UOJ #179】线性规划 单纯形模板
http://uoj.ac/problem/179 终于写出来了单纯性算法的板子,抄的网上大爷的qwq 辅助线性规划找非基变量时要加个随机化才能A,我也不知道为什么,卡精度吗? 2017-3-6UPD ...
- 【UOJ#179】线性规划 单纯形
题目链接: http://uoj.ac/problem/179 Solution 就是单纯形模板题,这篇博客就是存一下板子. Code #include<iostream> #includ ...
- 【UOJ 179】 #179. 线性规划 (单纯形法)
http://uoj.ac/problem/179 补充那一列修改方法: 对于第i行: $$xi=bi-\sum Aij*xj$$ $$=bi-\sum_{j!=e} Aij*xj-Aie*xe ...
- UVA 10498 Happiness(线性规划-单纯形)
Description Prof. Kaykobad has given Nasa the duty of buying some food for the ACM contestents. Nasa ...
- Note -「线性规划」学习笔记
\(\mathcal{Definition}\) 线性规划(Linear Programming, LP)形式上是对如下问题的描述: \[\operatorname{maximize}~~~~z= ...
- 线性规划VB求解
线性规划VB求解 Rem 定义动态数组 Dim a() As Single, c() As Single, b() As Single, cb() As Single Dim aa() As Sing ...
随机推荐
- dynamic
dynamic的特性很多,好像和反射也有关,不过这里先介绍一个特性,关于反射的再补充. 我们来看一个方法: public virtual ActionResult Insert(T info) 有一个 ...
- 【原创】backbone1.1.0源码解析之View
作为MVC框架,M(odel) V(iew) C(ontroler)之间的联系是必不可少的,今天要说的就是View(视图) 通常我们在写逻辑代码也好或者是在ui组件也好,都需要跟dom打交道,我们 ...
- 用python处理文本,本地文件系统以及使用数据库的知识基础
主要是想通过python之流的脚本语言来进行文件系统的遍历,处理文本以及使用简易数据库的操作. 本文基于陈皓的:<程序员技术练级攻略> 一.Python csv 对于电子表格和数据库导出文 ...
- JD m站自我解析理解
第一步:从首页着手 文档部分:应用的是H5默认文档开头 即:<!DOCUMENT html> head部分:放了一些相关的JS,title描述,然后就是meta表述了.比较有参考的如下 & ...
- DDLog设置方法
CHENYILONG Blog DDLog设置方法 本文永久地址为http://www.cnblogs.com/ChenYilong/p/3984246.html,转载请注明出处. 201 ...
- 【转】XMPP_3920_最靠谱的中文翻译文档
CHENYILONG Blog XMPP_3920_最靠谱的中文翻译文档 Fullscreen © chenyilong. Powered by Postach.io Blog
- 经典设计模式-iOS的实现
最近看了<HeadFirst 设计模式>这本书,给组内伙伴准备一次分享,把这次分享记录下来,有需要的可以看看. 这本书主要介绍了四人帮23种经典设计模式中的的14种,也是常用的几种.看完这 ...
- [转]Restrict关键字
0 定义 C99中新增加的用于修饰指针的关键字,用于表示该指针所指向的内存,只有通过该指针访问得到(如下ptr指向的内存单元只能通过ptr访问得到).从而可以让编译器对代码进行优化,生成更有效率的汇编 ...
- struts与ognl结合【重要】
-----------------------------ognl表达式------------------------ OGNL:对象视图导航语言. ${user.addr.name} 这种写法就 ...
- python3之模块psutil系统性能信息
psutil是个跨平台库,能够轻松实现获取系统运行的进程和系统利用率,包括CPU.内存.磁盘.网络等信息. 它主要应用于信息监控,分析和限制系统资源及进程的管理.它实现了同等命令命令行工具提供的功能, ...