bzoj3550: [ONTAK2010]Vacation（单纯形法+线性规划）

传送门

直接暴力把线性规划矩阵给打出来然后单纯形求解就行了
简单来说就是每个数记一个\(d_i\)表示选或不选，那么就是最大化\(\sum d_ic_i\)，并满足一堆限制条件
然后不要忘记限制每个数最多选一次
（据说还可以费用流然而实在不会啊……）

//minamoto
#include<bits/stdc++.h>
#define R register
#define inf 1e18
#define fp(i,a,b) for(R int i=a,I=b+1;i<I;++i)
#define fd(i,a,b) for(R int i=a,I=b-1;i>I;--i)
#define go(u) for(int i=head[u],v=e[i].v;i;i=e[i].nx,v=e[i].v)
using namespace std;
char buf[1<<21],*p1=buf,*p2=buf;
inline char getc(){return p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++;}
int read(){R int res,f=1;R char ch;while((ch=getc())>'9'||ch<'0')(ch=='-')&&(f=-1);for(res=ch-'0';(ch=getc())>='0'&&ch<='9';res=res*10+ch-'0');return res*f;
}
const int N=2005;const double eps=1e-8;
double a[N][N];int n,k,m;
void pivot(int l,int e){double t=a[l][e];a[l][e]=1;fp(i,0,m)a[l][i]/=t;fp(i,0,n)if(i!=l&&fabs(a[i][e])>eps){t=a[i][e],a[i][e]=0;fp(j,0,m)a[i][j]-=t*a[l][j];}
}
void simplex(){while(true){int l=0,e=0;double mn=inf;fp(i,1,m)if(a[0][i]>eps){e=i;break;}if(!e)break;fp(i,1,n)if(a[i][e]>eps&&a[i][0]/a[i][e]<mn)mn=a[i][0]/a[i][e],l=i;pivot(l,e);}
}
int main(){
//  freopen("testdata.in","r",stdin);n=read(),k=read();fp(i,1,n*3)a[0][i]=read();fp(i,1,n*3)a[i][i]=a[i][0]=1;fp(i,1,n*2+1){fp(j,0,n-1)a[i+n*3][i+j]=1;a[i+n*3][0]=k;}m=3*n,n=5*n+1;simplex();printf("%d\n",(int)(-a[0][0]+0.5));return 0;
}

转载于:https://www.cnblogs.com/bztMinamoto/p/10103013.html

bzoj3550: [ONTAK2010]Vacation（单纯形法+线性规划）相关推荐

Bzoj3550 [ONTAK2010]Vacation
Time Limit: 10 Sec Memory Limit: 96 MB Submit: 286 Solved: 209 Description 有3N个数,你需要选出一些数,首先保证任意长度 ...
线性规划问题及单纯形法-线性规划问题的求解方法
线性规划问题的求解方法两种方法 1.图解法(两个变量使用直角坐标.三个变量使用立体坐标) 2.单纯形法(适用于任意变量,但需将一般形式编程标准形式) 2图解法建立直角坐标(x1,x2>=0) ...
线性规划问题及单纯形法-线性规划变标准形
2.线性规划变标准形线性规划模型的标准形式 (1)目标函数为求极大值 (2)所有功能约束条件(非负条件除外),都是等式 (3)右端常数项为非负 (4)决策变量为非负标准形转换方法 (1)目标函数值 ...
bzoj3550: [ONTAK2010]Vacation1283: 序列
给出一个长度为的正整数序列Ci,求一个子序列,使得原序列中任意长度为的子串中被选出的元素不超过K(K,M<=100) 个,并且选出的元素之和最大. 据说是什么经典区间带权限制问题? 有两种写 ...
【运筹学】线性规划人工变量法 ( 单纯形法总结 | 人工变量法引入 | 人工变量法原理分析 | 人工变量法案例 )
文章目录一.单纯形法总结二.人工变量法引入三.人工变量法案例四.线性规划标准型五.人工变量法六.人工变量法解分析一.单纯形法总结求解线性规划 , 使用的是单纯形法 ; 迭代转化 : 其 ...
【运筹学】线性规划人工变量法 ( 人工变量法案例 | 初始单纯形表 | 检验数计算 | 入基变量 | 出基变量 )
文章目录一.生成初始单纯形表二.计算非基变量检验数三.最优解判定四.选择入基变量五.选择出基变量六.更新单纯形表上一篇博客 [运筹学]线性规划人工变量法 ( 单纯形法总结 | 人工变量 ...
最优化——分析线性规划的对偶问题的等价性
文章目录最优化-对偶原理与对偶单纯形法线性规划的对偶原理原问题为何与对偶问题等价前提1 前提2 证明等价性最优化-对偶原理与对偶单纯形法线性规划的对偶原理对于标准线性规划问题: max⁡ ...
【运筹学】线性规划人工变量法阶段总结 ( 使用人工变量法求解线性规划全过程详细解析 ) ★
文章目录一.人工变量法案例二.线性规划标准型三.添加人工变量四.初始单纯形表五.初始单纯形表 : 计算非基变量检验数六.初始单纯形表 : 最优解判定七.初始单纯形表 : 选择入基变量八 ...
单纯形法的C语言实现
//haust //momentum 摘要:在运筹学中利用经典的单纯形法人工求解线性规划问题较为复杂,但利用计算机强大的计算能力便能够快速实现.本文利用经典的单纯形法的算法,借助传统的C语言进行实现, ...
运筹学第一章：线性规划【复习自用】
标题~ 对于以后期末复习线性规划写在前面:名词解释及关系线性规划(LP)问题可行解最优解基本解可行基/基本可行解定理:基本解与基本可行解最优基/最优基本可行解最优性判别定理退化/ ...

bzoj3550: [ONTAK2010]Vacation（单纯形法+线性规划）

bzoj3550: [ONTAK2010]Vacation（单纯形法+线性规划）相关推荐

最新文章

热门文章