数学(线性规划): ZJOI2013 防守战线

时间:2023-03-10 02:38:08
数学(线性规划): ZJOI2013 防守战线

数学(线性规划): ZJOI2013 防守战线

  偷懒用的线性规划。

 #include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;
const int maxr=;
const int maxc=; int n,m,nxt[maxc];
int a[maxr][maxc]; void Pivot(int l,int e){
int pre=maxc-;
for(int i=;i<=n;i++)
if(a[l][i]!=){nxt[pre]=i;pre=i;}
nxt[pre]=-; for(int i=,t;i<=m;i++)
if(i!=l&&(t=a[i][e])){
a[i][e]=;
for(int j=nxt[maxc-];j!=-;j=nxt[j])
a[i][j]+=t*a[l][j];
}
} void Simplex(){
while(true){
int e=,l=;
for(int i=;i<=n;i++)
if(a[][i]>){e=i;break;}
if(e==)break;
for(int i=;i<=m;i++)
if(a[i][e]<&&(!l||a[l][]>a[i][]))
{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=;i<=m;i++)
scanf("%d",&a[i][]);
for(int i=,l,r,t;i<=n;i++){
scanf("%d%d%d",&l,&r,&t);
for(int j=l;j<=r;j++)
a[j][i]=-;
a[][i]=t;
}
Simplex();
printf("%d\n",a[][]);
return ;
}