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。

Source

By Sbullet

[ Submit][ Status][ Discuss]

分析：
单纯形算法
因为每个数字只能选择一次，所以我们要在所有条件后加上$x_i<=1$

因为所求是最小值，约定条件都是<=
我们直接上模板即可

tip

一共有$5*n+1$个约定条件
不要搞错了

#include<cstdio>
#include<cstring>
#include<iostream>
#include<cmath>

using namespace std;

const double INF=1e10;
const double eps=1e-7;
const int N=203;
double a[5*N][3*N],b[5*N],c[3*N],ans;
int n,K,m;

void priov(int l,int e)
{
    b[l]/=a[l][e];
    for (int i=1;i<=m;i++) if (i!=e)
        a[l][i]/=a[l][e];
    a[l][e]=1/a[l][e];
    for (int i=1;i<=n;i++) if (i!=l&&fabs(a[i][e])>eps)
    {
        b[i]-=b[l]*a[i][e];
        for (int j=1;j<=m;j++) if (j!=e)
            a[i][j]-=a[l][j]*a[i][e];
        a[i][e]=-a[i][e]*a[l][e];
    }
    ans+=c[e]*b[l];
    for (int i=1;i<=m;i++) if (i!=e)
        c[i]-=c[e]*a[l][i];
    c[e]=-c[e]*a[l][e];
}

void simple()
{
    ans=0.0;
    int i,l,e;
    while (1)
    {
        for (i=1;i<=m;i++)
            if (c[i]>eps) break;
        e=i;
        if (e==m+1) return;
        double t=INF;
        for (int i=1;i<=n;i++)
            if (a[i][e]>eps&&t>b[i]/a[i][e])
                t=b[i]/a[i][e],l=i;
        if (t==INF) return;
        priov(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++)
    {
        for (int j=i;j<i+n;j++) a[i][j]++;
        b[i]=K;
    }
    for (int i=2*n+2;i<5*n+2;i++) a[i][i-2*n-1]++,b[i]=1;
    m=3*n;
    n*=5; n++;
    simple();
    printf("%.0lf",ans);    
    return 0;
}