bzoj 1061: [Noi2008]志愿者招募 单纯形

题意

申奥成功后，布布经过不懈努力，终于成为奥组委下属公司人力资源部门的主管。布布刚上任就遇到了一个难题：为即将启动的奥运新项目招募一批短期志愿者。经过估算，这个项目需要N 天才能完成，其中第i 天至少需要Ai 个人。 布布通过了解得知，一共有M 类志愿者可以招募。其中第i 类可以从第Si 天工作到第Ti 天，招募费用是每人Ci 元。新官上任三把火，为了出色地完成自己的工作，布布希望用尽量少的费用招募足够的志愿者，但这并不是他的特长！于是布布找到了你，希望你帮他设计一种最优的招募方案。
1 ≤ N ≤ 1000，1 ≤ M ≤ 10000，题目中其他所涉及的数据均 不超过2^31-1。

分析

裸的单纯形法。
推荐教程

代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<cmath>
#define inf 1e10
#define eps 1e-7
using namespace std;

const int N=1005;
const int M=10005;

int n,m;
double c[N],a[M][N],b[M],v;

void pivot(int l,int e)
{
    for (int i=1;i<=n;i++) if (i!=e) a[l][i]/=a[l][e];
    b[l]/=a[l][e];a[l][e]=1/a[l][e];
    for (int i=1;i<=m;i++)
        if (i!=l&&fabs(a[i][e])>eps)
        {
            b[i]-=b[l]*a[i][e];
            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];
        }
    v+=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];
}

double simplex()
{
    while (1)
    {
        int e=0,l;
        for (int i=1;i<=n;i++)
            if (c[e=i]>eps) break;
        if (c[e]<eps) return v;
        double tmp=inf;
        for (int i=1;i<=m;i++)
            if (a[i][e]>eps&&b[i]/a[i][e]<tmp) tmp=b[i]/a[i][e],l=i;
        pivot(l,e);
    }
}

int main()
{
    scanf("%d%d",&n,&m);
    for (int i=1;i<=n;i++) scanf("%lf",&c[i]);
    for (int i=1;i<=m;i++)
    {
        int x,y,z;
        scanf("%d%d%d",&x,&y,&z);
        for (int j=x;j<=y;j++) a[i][j]=1;
        b[i]=z;
    }
    double ans=simplex();
    printf("%d",(int)(ans+0.1));
    return 0;
}