1061: [Noi2008]志愿者招募

题目链接

题目大意：给定n天，第i天需要ai个志愿者，有m类志愿者，每类志愿者工作时间为[l,r]，花费为ci，求最小花费

题解：对偶一下变成标准型，单纯形法

我的收获：2333

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
typedef long long ll;

const double eps=1e-8,INF=1e15;

int n,m;
double a[10005][1005];

void Pivot(int l,int e){
    double t=a[l][e];a[l][e]=1;
    for(int i=0;i<=n;i++) a[l][i]/=t;
    for(int i=0;i<=m;i++){
        if(i!=l&&abs(a[i][e])>eps){
            t=a[i][e];a[i][e]=0;
            for(int j=0;j<=n;j++) a[i][j]-=a[l][j]*t;
        }
    }
}

void simplex(){
    while(true){
        int l=0,e=0;double mn=INF;
        for(int j=1;j<=n;j++) if(a[0][j]>eps){e=j;break;}
        if(!e) return ;
        for(int i=1;i<=m;i++) if(a[i][e]>eps&&a[i][0]/a[i][e]<mn)
            mn=a[i][0]/a[i][e],l=i;
        if(!l) return ;
        Pivot(l,e);
    }
}

void work()
{
    simplex(); 
    printf("%.0lf\n",-a[0][0]);
}

void init()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++) scanf("%lf",&a[0][i]);
    for(int i=1;i<=m;i++){
        int x,y;scanf("%d%d%lf",&x,&y,&a[i][0]); 
        for(int j=x;j<=y;j++) a[i][j]=1;
    }
}

int main()
{
    init();
    work();
    return 0;
}