[BZOJ3112][ZJOI2013]防守战线-线性规划

防守战线

Description

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

Input

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

Output

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

Sample Input

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

Sample Output

11

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。

Source

ZJOI2013

---

当你不需要知道对偶的具体原理而是直接使用它的话，你会发现对偶是如此的方便……

---

思路:
可以轻松地建出一个线性规划模型:

最小化 $\sum\limits_{i=1}^n c_ix_i$
满足

$\forall j \in [1,m],\sum\limits_{i=l_j}^{r_j}x_i \leqslant d_j$

然后可以发现这是最小化，而不是单纯形需要的最大化。
于是无脑对偶一波，直接上单纯形板子，收工。
对于咱的这种保存信息的方式，对偶时直接令原本的所有$a[i][j]$放在$a[j][i]$就可以了，什么别的都不用管。

注意常数优化，不加优化会TLE，加了优化只要2s左右就能过……

#include<bits/stdc++.h>
using namespace std;

typedef double db;
const int N=1009;
const int M=10009;
const db eps=1e-7;

db a[N][M];
int n,m,id[N+M];

inline int read()
{
    int x=0;char ch=getchar();
    while(ch<'0' || '9'<ch)ch=getchar();
    while('0'<=ch && ch<='9')x=x*10+(ch^48),ch=getchar();
    return x;
}

namespace simplex
{
    int q[N];
    inline void pivot(int l,int e)
    {
        swap(id[e],id[n+l]);
        db t=a[l][e];
        a[l][e]=1.0;
        for(int i=0;i<=n;i++)
            a[l][i]/=t;

        int top=0;
        for(int i=0;i<=n;i++)
            if(fabs(a[l][i])>eps)
                q[++top]=i;
        for(int i=0;i<=m;i++)
            if(i!=l && fabs(a[i][e])>eps)
            {
                int t=a[i][e];
                a[i][e]=0.0;
                for(int j=1;j<=top;j++)
                    a[i][q[j]]-=t*a[l][q[j]];
            }
    }

    inline bool init()
    {
        while(1)
        {
            int l=0,e=0;
            db k=-eps;
            for(int i=1;i<=m;i++)
                if(a[i][0]<k)
                    k=a[i][0],l=i;
            if(!l)return 1;
            for(int i=1;i<=n && !e;i++)
                if(a[l][i]<-eps)
                    e=i;
            if(!e)return 0;
            pivot(l,e);
        }
    }

    inline bool simplex()
    {
        while(1)
        {
            int l=0,e=0;
            for(int i=1;i<=n && !e;i++)
                if(a[0][i]>eps)
                    e=i;
            if(!e)return 1;
            db k=1e18;
            for(int i=1;i<=m;i++)
                if(a[i][e]>eps && a[i][0]/a[i][e]<k)
                    k=a[i][0]/a[i][e],l=i;
            if(!l)return 0;
            pivot(l,e);
        }
    }

    inline int mina()
    {
        if(!init())return 0;
        if(!simplex())return 0;
        return (int)(-a[0][0]+0.5);
    }
}

int main()
{
    m=read();n=read();
    for(int i=1;i<=m;i++)
        a[i][0]=read();
    for(int i=1,l,r;i<=n;i++)
    {
        l=read();r=read();
        a[0][i]=read();
        for(int j=l;j<=r;j++)
            a[j][i]=1;
    }

    printf("%d\n",simplex::mina());
    return 0;
}