BZOJ 3550: [ONTAK2010]Vacation [单纯形法]

有3N个数，你需要选出一些数，首先保证任意长度为N的区间中选出的数的个数<=K个，其次要保证选出的数的个数最大。

---

好像都是费用流...

单纯性裸题呀...

注意每个数最多选1次

#include <iostream>

#include <cstdio>

#include <cstring>

#include <algorithm>

#include <cmath>

using namespace std;

typedef long long ll;

const int M=,N=;

const double INF=1e15,eps=1e-;

inline int read(){

    char c=getchar();int x=,f=;

    while(c<''||c>''){if(c=='-')f=-;c=getchar();}

    while(c>=''&&c<=''){x=x*+c-'';c=getchar();}

    return x*f;

}

int n,m,k;

double a[M][N];

int q[N];

void Pivot(int l,int e){

    double t=a[l][e]; a[l][e]=;

    for(int j=;j<=n;j++) a[l][j]/=t;

    int p=;

    for(int j=;j<=n;j++) if(abs(a[l][j])>eps) q[++p]=j;

    for(int i=;i<=m;i++) if(i!=l && abs(a[i][e])>eps){

        double t=a[i][e]; a[i][e]=;

        for(int j=;j<=p;j++) a[i][q[j]]-=t*a[l][q[j]];

    }

}

void simplex(){

    while(true){

        int l=,e=; double mn=INF;

        for(int j=;j<=n;j++) if(a[][j]>eps) {e=j;break;}

        if(!e) return;

        for(int i=;i<=m;i++)

            if(a[i][e]>eps && a[i][]/a[i][e]<mn) l=i,mn=a[i][]/a[i][e];

        if(!l) return;//unbounded

        Pivot(l,e);

    }

}

int main(){

    freopen("in","r",stdin);

    n=read();k=read();

    m=*n+;int _=*n;

    for(int i=;i<=_;i++) a[][i]=read();

    for(int i=;i<=m;i++){

        a[i][]=k;

        for(int j=i;j<=i+n-;j++) a[i][j]=;

    }

    m=*n+;

    for(int i=*n+;i<=m;i++)

        a[i][i-(n<<)-]=,a[i][]=;

    n=*n;

    simplex();

    printf("%.0lf",-a[][]);

}