单纯型法Ⅱ（bzoj 1061: [Noi2008]志愿者招募）

线性规划单纯型法：http://blog.csdn.net/jaihk662/article/details/78050666

C是一个n的向量，B是一个m的向量，在满足约束情况下最大化Cx

#include int n, m; double ans, c[1005], b[10005], a[10005][1005]; void Pviot(int l, int e) { int i, j; a[l][e] = 1/a[l][e]; b[l] *= a[l][e]; for(i=1;i<=n;i++) { if(i^e) a[l][i] *= a[l][e]; } for(i=1;i<=m;i++) { if((i^l) && a[i][e]) { for(j=1;j<=n;j++) { if(j^e) a[i][j] -= a[i][e]*a[l][j]; } b[i] -= a[i][e]*b[l],a[i][e]*=-a[l][e]; } } for(i=1;i<=n;i++) { if(i^e) c[i] -= c[e]*a[l][i]; } ans += c[e]*b[l]; c[e] *= -a[l][e]; } double SimpleX() { int p, i, temp; double k; while(1) { for(p=1;p<=n;p++) { if(c[p]>0) break; } if(p==n+1) return ans; k = 1e10, temp = 0; for(i=1;i<=m;i++) { if(a[i][p]>0 && b[i]/a[i][p]

1061: [Noi2008]志愿者招募

Time Limit: 20 Sec Memory Limit: 162 MB
Submit: 5052 Solved: 3012
[ Submit][ Status][ Discuss]

Description

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

Input

　　第一行包含两个整数N, M，表示完成项目的天数和可以招募的志愿者的种类。接下来的一行中包含N 个非负整数，表示每天至少需要的志愿者人数。接下来的M 行中每行包含三个整数Si, Ti, Ci，含义如上文所述。为了方便起见，我们可以认为每类志愿者的数量都是无限多的。

Output

　　仅包含一个整数，表示你所设计的最优方案的总费用。

Sample Input

3 3

2 3 4

1 2 2

2 3 5

3 3 2

Sample Output

设需要第i个志愿者xi个，对于样例可以得出不等式组：

①x1>=2；②x1+x2>=3；③x2+x3>=4

其中的G() = 2*x1+5*x2+2*x3，求G()的最小值

方法1：

添加辅助变量(x4, x5, x6)：①x1-x4=2；②x1+x2-x5=3；③x2+x3-x6=4：④0=0

拿所有的式子减去它上一个式子：①x1-x4=2；②x2+x4-x5=1；③x3-x1+x5-x6=1；④x6-x2-x3=-4

然后就可以最小费用最大流了，建图方法看上面链接（因为对应系数矩阵为列连续，所以合法）

方法2：

看系数矩阵：

2 5 2(求min)

1 2（>=）

1 1 3

1 1 4

对偶矩阵（翻转90°）：

2 3 4(求max)

1 1 2（<=）

1 1 5

1 2

这样就可以套模板了，AC代码就是上面那个