BZOJ 1061 Noi2008 志愿者招募 单纯形

wyfcyx的《线性规划与单纯形算法》

http://wenku.baidu.com/view/ce5784754a7302768f99391d

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#define maxn 1010
#define maxm 10010
#define inf 2147483647
using namespace std;
long long a[maxn][maxm],next[maxm];
int n,m;

void pivot(int l,int e)
{
    int last=-1;
    for (int i=0;i<=m;i++)
        if (a[l][i])
        {
            next[i]=last;
            last=i;
        }
    for (int i=0;i<=n;i++)
    {
        if (a[i][e]==0||i==l) continue;
        for (int j=last;j!=-1;j=next[j])
        {
            if (j==e) continue;
            a[i][j]-=a[i][e]*a[l][j];
        }
        a[i][e]=-a[i][e];
    }
}

long long simplex()
{
    while (1)
    {
        int now=0;
        for (int i=1;i<=m;i++) 
            if (a[0][i]>0) { now=i; break; }
        if (now==0) return -a[0][0];
        int tmp,mi=inf;
        for (int i=1;i<=n;i++)
        {
            if (a[i][now]>0&&a[i][0]<mi)
            {
                tmp=i;
                mi=a[i][0];
            }
        }
        pivot(tmp,now);
    }
}
int main()
{
    scanf("%d%d", &n, &m);
    for(int i=1;i<=n;i++) scanf("%lld", &a[0][i]);
    int l, r, d;
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d%lld", &l, &r, &d);
        a[i][0] = d;
        for(int j=l;j<=r;j++) a[i][j] = 1;
    }
    long long  ans = simplex();
    printf("%lld\n", ans);
    return 0;
}