【bzoj3550】[ONTAK2010]Vacation 单纯形

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。

题解
单纯形裸题。
根据条件构造即可。

代码

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<set>
#include<complex>
#define ll long long
#define mod 10000007
using namespace std;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if (ch=='-') f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
#define eps (1e-10)
const int N=1005;
int n,m,k;
double a[N][N],b[N],c[N],ans;
inline void pivot(int l,int e)
{
    b[l]/=a[l][e];
    for (int i=1;i<=n;i++) if (i!=e) a[l][i]/=a[l][e];
    a[l][e]=1.0/a[l][e];
    for (int i=1;i<=m;i++)if (i!=l&&fabs(a[i][e])>eps)
    {
        b[i]-=a[i][e]*b[l];
        for (int j=1;j<=n;j++) if (j!=e) a[i][j]-=a[i][e]*a[l][j];
        a[i][e]=-a[i][e]*a[l][e];
    }
    ans+=c[e]*b[l];
    for (int i=1;i<=n;i++)if (i!=e) c[i]-=c[e]*a[l][i];
    c[e]=-c[e]*a[l][e];
}
void simplex()
{
    int l,e;
    while (1)
    {
        for (e=1;e<=n;e++) if (c[e]>eps) break;
        if (e==n+1) break;
        double t=10000000000;
        for (int i=1;i<=m;i++)if (a[i][e]>eps&&b[i]/a[i][e]<t)
        {
            t=b[i]/a[i][e];
            l=i;
        }
        pivot(l,e);
    }
}
int main()
{
    scanf("%d%d",&n,&k);
    for (int i=1;i<=3*n;++i) scanf("%lf",&c[i]);
    for (int i=1;i<=2*n+1;++i)
    {
        b[i]=k+0.0;
        for (int j=i;j<=i+n-1;++j)
            a[i][j]=1.0;
    }
    for (int i=2*n+2;i<=5*n+1;++i)
        b[i]=1,a[i][i-2*n-1]=1;
    m=5*n+1;n*=3;
    simplex();
    printf("%.0lf",ans+eps);
    return 0;
}