HYSBZ 1061,单纯形&对偶原理

志愿者招募

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

方案一：单纯形
思路：设第i类人招募了xi个，aij表示第i个人第j天是否工作，工作为1，不工作为0，则应满足对于所有j，${\sum }_i^n{a}_{ij}{x}_i>=A_j$，需要最小化$z={\sum }_1^n{C}_i{x}_i$，根据对偶原理将矩阵转置后变为最大化即可

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int MAXN=10005,MAXM=1005,INF=1e9;
const double eps=1e-6;
int n,m;
double a[MAXN][MAXM],b[MAXN],c[MAXM],v;
//a为约束矩阵，有n行m列，c最大化z的系数，b为限制
inline void pivot(int l,int e)
{
    b[l]/=a[l][e];
    for(int j=1;j<=m;j++) if(j!=e) a[l][j]/=a[l][e];
    a[l][e]=1/a[l][e];
    for(int i=1;i<=n;i++) if(i!=l&&fabs(a[i][e])>0)
    {
            b[i]-=a[i][e]*b[l];
            for(int j=1;j<=m;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 j=1;j<=m;j++) if(j!=e) c[j]-=c[e]*a[l][j];
    c[e]=-c[e]*a[l][e];
}

inline double simplex()
{
    while(true)
    {
        int e=0,l=0;
        for(e=1;e<=m;e++) if(c[e]>eps) break;
        if(e==m+1) return v;
        double mn=INF;
        for(int i=1;i<=n;i++)
            if(a[i][e]>eps&&mn>b[i]/a[i][e]) mn=b[i]/a[i][e],l=i;
        if(mn==INF) return INF;//无穷大
        pivot(l,e);
    }
}

int main()
{
    while(~scanf("%d%d",&m,&n))
    {
        for(int i=1;i<=m;++i) scanf("%lf",&c[i]);
        int s,t;
        for(int i=1;i<=n;++i)
        {
            scanf("%d%d",&s,&t);
            for(int j = 1;j < s;++j) a[i][j] = 0;
            for(int j = s;j <= t;++j) a[i][j] = 1;
            for(int j = t+1;j <= m;++j) a[i][j] = 0;
            scanf("%lf",&b[i]);
        }
        printf("%.0f",simplex());
    }
}

方案二：最小费用流
https://blog.csdn.net/xing_mo/article/details/103986110