数学（线性规划）： ZJOI2013 防守战线

#include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;
const int maxr=1010;
const int maxc=10010;

int n,m,nxt[maxc];
int a[maxr][maxc];

void Pivot(int l,int e){
int pre=maxc-1;
for(int i=0;i<=n;i++)
if(a[l][i]!=0){nxt[pre]=i;pre=i;}
nxt[pre]=-1;

for(int i=0,t;i<=m;i++)
if(i!=l&&(t=a[i][e])){
a[i][e]=0;
for(int j=nxt[maxc-1];j!=-1;j=nxt[j])
a[i][j]+=t*a[l][j];
}
}

void Simplex(){
while(true){
int e=0,l=0;
for(int i=1;i<=n;i++)
if(a[0][i]>0){e=i;break;}
if(e==0)break;
for(int i=1;i<=m;i++)
if(a[i][e]<0&&(!l||a[l][0]>a[i][0]))
{l=i;}

Pivot(l,e);
}
}

int main(){
#ifndef ONLINE_JUDGE
freopen("zjoi13_defend.in","r",stdin);
freopen("zjoi13_defend.out","w",stdout);
#endif
scanf("%d%d",&m,&n);
for(int i=1;i<=m;i++)
scanf("%d",&a[i][0]);
for(int i=1,l,r,t;i<=n;i++){
scanf("%d%d%d",&l,&r,&t);
for(int j=l;j<=r;j++)
a[j][i]=-1;
a[0][i]=t;
}
Simplex();
printf("%d\n",a[0][0]);
return 0;
}