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[][]);
}
上一篇:eclipse 鲜为人知的调试技巧,你用过多少


下一篇:迷宫 (BFS)