bzoj3550 [ONTAK2010]Vacation

最新推荐文章于 2018-12-06 13:57:59 发布

原创最新推荐文章于 2018-12-06 13:57:59 发布 · 186 阅读

0 ·

CC 4.0 BY-SA版权

线性规划专栏收录该内容

4 篇文章

订阅专栏

本文介绍了解决线性规划问题的一种有效方法——单纯形法，并通过一个具体实例展示了如何利用该方法求解最大化的线性规划问题。文中提供了一段C++代码实现，用于演示单纯形法的具体应用。

http://www.elijahqi.win/archives/2905
Description
有3N个数，你需要选出一些数，首先保证任意长度为N的区间中选出的数的个数<=K个，其次要保证选出的数的个数最大。

Input
第一行两个整数N，K。
第二行有3N个整数。

Output

一行一个整数表示答案。

Sample Input
5 3

14 21 9 30 11 8 1 20 29 23 17 27 7 8 35

Sample Output
195

HINT

【数据范围】

N<=200，K<=10。

Source
By Sbullet

打遍板子复习下单纯形题解就不写了，看蒟蒻我写的前几篇单纯形的题目

#include<cmath>
#include<cctype>
#include<cstdio>
#define N 1100
#define eps 1e-7
#define inf 0x3f3f3f3f
#include<algorithm>
using namespace std;
inline char gc(){
    static char now[1<<16],*S,*T;
    if (T==S){T=(S=now)+fread(now,1,1<<16,stdin);if (T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(!isdigit(ch)) {if (ch=='-') f=-1;ch=gc();}
    while(isdigit(ch)) x=x*10+ch-'0',ch=gc();
    return x*f;
}
int n,m,k;
double a[N][N];
inline void pivot(int l,int e){
    double t=a[l][e];a[l][e]=1;
    for (int i=0;i<=n;++i) a[l][i]/=t;
    for (int i=0;i<=m;++i){
        if (i==l||abs(a[i][e])<eps) continue;
        t=a[i][e];a[i][e]=0;
        for (int j=0;j<=n;++j) a[i][j]-=t*a[l][j];
    }
}
inline void simplex(){
    while(1){
        int l=0,e=0;
        for (int i=1;i<=n;++i) if (a[0][i]>eps) {e=i;break;}
        if (!e) return;double mn=inf;
        for (int i=1;i<=m;++i) if (a[i][e]>eps&&mn>a[i][0]/a[i][e]) mn=a[i][0]/a[i][e],l=i;
        if (!l) return;pivot(l,e);
    }
}
int main(){
    freopen("bzoj3550.in","r",stdin);
    n=read();m=n;k=read();
    for (int i=1;i<=n*3;++i) a[0][i]=read();
    for (int i=1;i<=n*3-m+1;++i){
        for (int j=1;j<=m;++j) a[i][i+j-1]=1;
        a[i][0]=k;
    }m=n*3-m+1;for (int i=1;i<=n*3;++i) a[++m][0]=1,a[m][i]=1;n*=3;
    simplex();printf("%d\n",-(int)a[0][0]);
    return 0;
}