bzoj3112 [Zjoi2013]防守战线（单纯形）

题目描述 Description

战线可以看作一个长度为 n 的序列，现在需要在这个序列上建塔来防守敌
兵，在序列第 i号位置上建一座塔有 Ci 的花费，且一个位置可以建任意多的塔
费用累加计算。有 m个区间[L1, R1], [L2, R2], …, [Lm, Rm]，在第 i 个区间
的范围内要建至少 Di座塔。求最少花费。

输入描述 Input Description

第一行为两个数n,m。
接下来一行，有 n个数，描述 C数组。
接下来 m行，每行三个数 Li,Ri,Di，描述一个区间。

输出描述 Output Description

仅包含一行，一个数，为最少花费。

样例输入 Sample Input

5 3
1 5 6 3 4
2 3 1
1 5 4
3 5 2

样例输出 Sample Output

11

数据范围及提示 Data Size & Hint

位置 1建 2个塔，位置 3建一个塔，位置 4建一个塔。花费 1*2+6+3=11。

对于 20%的数据，n≤20，m≤20。
对于 50%的数据（包括上部分的数据），Di 全部为1。
对于 70%的数据（包括上部分的数据），n≤100，m≤1000。
对于 100%的数据，n≤1000，m≤10000，1≤Li≤Ri≤n，其余数据均≤10000。

分析：
单纯形算法
因为约束条件都是大于等于，目标函数是最小值，所以我们要把线性规划转化成对偶图

在构建对偶图的时候要想明白了
A数组应该是一个n*m（n行m列）的矩阵 （这里说的m和n与题目中的意义不同）
变成对偶图之后，就会变成一个m*n的矩阵
B数组表示约束条件的常数
C数组表示目标函数的系数
构建对偶图，B和C数组要互换

tip

数据范围

#include<cstdio>
#include<cstring>
#include<iostream>
#include<cmath> 

using namespace std;

const double eps=1e-8;
const double INF=1e10;
const int N=1005;
const int M=10005;
double a[N][M],b[N],c[M],ans;
int n,m;

void priov(int l,int e)
{
    b[l]/=a[l][e];
    for (int i=1;i<=m;i++) if (i!=e)
        a[l][i]/=a[l][e];
    a[l][e]=1/a[l][e];
    for (int i=1;i<=n;i++) if (i!=l&&fabs(a[i][e])>eps)   //一定要有变量e 
    {
        b[i]-=b[l]*a[i][e];
        for (int j=1;j<=m;j++) if (j!=e)
            a[i][j]-=a[l][j]*a[i][e];
        a[i][e]=-a[i][e]*a[l][e];
    }
    ans+=c[e]*b[l];
    for (int i=1;i<=m;i++) if (i!=e)
        c[i]-=c[e]*a[l][i];
    c[e]=-c[e]*a[l][e];
}

void simple()
{
    ans=0.0;
    int i,l,e;
    while (1)
    {
        for (i=1;i<=m;i++) 
            if (c[i]>eps) break;
        e=i;
        if (e==m+1) return;
        double t=INF;
        for (i=1;i<=n;i++)
            if ((a[i][e]>eps)&&(t>b[i]/a[i][e]))
                t=b[i]/a[i][e],l=i;
        if (t==INF) return;
        priov(l,e);
    }
}

int main()
{
    scanf("%d%d",&n,&m);
    for (int i=1;i<=n;i++) scanf("%lf",&b[i]);
    for (int i=1;i<=m;i++)
    {
        int l,r;
        scanf("%d%d",&l,&r); 
        for (int j=l;j<=r;j++) a[j][i]++;    //a[j][i] 对偶图 
        scanf("%lf",&c[i]);
    }
    simple();
    printf("%.0lf",ans);
    return 0;
}