bzoj1061 [Noi2008]志愿者招募 单纯形（模板）

题目大意：
有n天，m种志愿者。
每天需要ai个志愿者，雇佣每种志愿者的费用为ci，每个志愿者可以在Si天到Ti天工作。
问最小花费。

题目分析：
题解PoPoQQQ.html
单纯形wyfcyx.pdf

据说这道题网络流神题，单纯形裸题

首先我们根据题意
设第i天需要xi个志愿者。
$A{_i,_j}$代表第j个人可以在第i天工作。
$b{_i}$代表第i天需要的人数。
$c{_i}$代表租用第i类志愿者需要的花费。
$\sum{A{_i,_j}x{_j}}\geq b{_i}$
$Minisize \sum{_i^m}{c{_i}*x{_i}}$

这个不是单纯形的标准形式，对偶原理搞一波就变成了
$\sum{A{^T_i,_j}*y{_j}}\leq c{_i}$
$Maxisize \sum{_i^n}{b{_i}*y{_i}}$

这样就变成了单纯形裸题QWQ

（向Po姐学习！！！）
代码如下：

#include <cstdio>
#include <cmath>
#define N 1010
#define M 10100
using namespace std;
const double INF=1e10;
const double eps=1e-7;
int n,m;
double A[M][N],b[M],c[N],v;
void Pivot(int l,int e)
{
    b[l]/=A[l][e];
    for(int i=1;i<=n;i++)
        if(i!=e)
            A[l][i]/=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]-=A[i][e]*b[l];
            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()
{
    int i,l,e;
    while(true)
    {
        for(i=1;i<=n;i++)
            if(c[i]>eps) break;
        if((e=i)==n+1) return v;
        double tmp=INF;
        for(i=1;i<=m;i++)
            if(A[i][e]>eps && b[i]/A[i][e]<tmp)
                tmp=b[i]/A[i][e],l=i;
        if(tmp==INF) return INF;
        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,x,y;i<=m;i++)
    {
        scanf("%d%d%lf",&x,&y,&b[i]);
        for(int j=x;j<=y;j++)
            A[i][j]=1;
    }
    double ans=Simplex();
    printf("%d\n",int(ans+0.5));
    return 0;
}