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=1005,N=605;
const double INF=1e15,eps=1e-8;
inline int read(){
    char c=getchar();int x=0,f=1;
    while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
    while(c>='0'&&c<='9'){x=x*10+c-'0';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]=1;
    for(int j=0;j<=n;j++) a[l][j]/=t;
    int p=0;
    for(int j=0;j<=n;j++) if(abs(a[l][j])>eps) q[++p]=j;
    for(int i=0;i<=m;i++) if(i!=l && abs(a[i][e])>eps){
        double t=a[i][e]; a[i][e]=0;
        for(int j=1;j<=p;j++) a[i][q[j]]-=t*a[l][q[j]];
    }
}
void simplex(){
    while(true){
        int l=0,e=0; double mn=INF;
        for(int j=1;j<=n;j++) if(a[0][j]>eps) {e=j;break;}
        if(!e) return;
        for(int i=1;i<=m;i++) 
            if(a[i][e]>eps && a[i][0]/a[i][e]<mn) l=i,mn=a[i][0]/a[i][e];
        if(!l) return;//unbounded
        Pivot(l,e);
    }
}
int main(){
    freopen("in","r",stdin);
    n=read();k=read();
    m=2*n+1;int _=3*n;
    for(int i=1;i<=_;i++) a[0][i]=read(); 
    for(int i=1;i<=m;i++){
        a[i][0]=k; 
        for(int j=i;j<=i+n-1;j++) a[i][j]=1;
    }
    m=5*n+1;
    for(int i=2*n+2;i<=m;i++)
        a[i][i-(n<<1)-1]=1,a[i][0]=1;
    n=3*n;
    simplex();
    printf("%.0lf",-a[0][0]);
}