【BZOJ3112】【ZJOI2013】—防守战线（线性规划+对偶）

传送门

---

Solution：

题意即求：
$$Min\sum _{i=1}^n{C}_i{x}_i$$
满足约束$$\sum _{l_j\le i\le r_j}{x}_i\ge D_j,j\in [1,m]$$
$$x_i\ge 0$$

对偶得：

$$Max\sum _{j=1}^m{D}_j{y}_j$$
满足约束$$\sum _{l_j\le i\le r_j}{y}_j\le C_i.i\in [1,n]$$
$$y_j\ge 0$$

即对于一个$C_i$来说，所有包含它的区间的权值和小于它
对于每一个$[l,r]$竖着填就是了

#include<bits/stdc++.h>
using namespace std;
const int RLEN=1<<20|1;
inline char gc(){
    static char ibuf[RLEN],*ib,*ob;
    (ib==ob)&&(ob=(ib=ibuf)+fread(ibuf,1,RLEN,stdin));
    return ib==ob?EOF:*ib++;
}
#define gc getchar
#define pb push_back
inline int read(){
    char ch=gc();
    int res=0,f=1;
    while(!isdigit(ch))f^=ch=='-',ch=gc();
    while(isdigit(ch))res=(res+(res<<2)<<1)+(ch^48),ch=gc();
    return f?res:-res;
}
const int N=10005,M=1005;
int n,m;
double f[M][N];
const double eps=1e-6,inf=1e18;
inline void povit(int l,int e){
	double t=f[l][e];
	for(int j=0;j<=n;j++)f[l][j]/=t;
	for(int i=0;i<=m;i++)if(l!=i&&fabs(f[i][e])>0){
		t=f[i][e],f[i][e]=0;
		for(int j=0;j<=n;j++)f[i][j]-=t*f[l][j];
	}
}
inline void simplex(){
	while(1){
		int l=0,e=0;double mn=inf;
		for(int i=1;i<=n;i++)if(f[0][i]>eps){e=i;break;}
		if(!e)break;
		for(int i=1;i<=m;i++)if(f[i][e]>eps&&f[i][0]/f[i][e]<mn)
			mn=f[i][0]/f[i][e],l=i;
		povit(l,e);
	}
}
int main(){
	m=read(),n=read();
	for(int i=1;i<=m;i++)f[i][0]=read();
	for(int i=1;i<=n;i++){
		int l=read(),r=read();f[0][i]=read();
		for(int j=l;j<=r;j++)f[j][i]=1;
	}
	simplex();
	cout<<(int)-f[0][0];
}