单纯形转对偶图。。
由于是全幺模,直接用int存就好了。
常数没有优化慢的飞起。。。
AC代码如下:
#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int n,m,a[1005][10005];
void pvt(int x,int y){
int i,j;
for (i=0; i<=n; i++) if (i!=x && a[i][y]){
for (j=0; j<=m; j++) if (a[x][j] && j!=y)
a[i][j]-=a[i][y]*a[x][j];
a[i][y]=-a[i][y];
}
}
int solve(){
int i,j,k; a[n+1][0]=1000000000;
while (1){
for (i=1; i<=m; i++) if (a[0][i]>0) break;
if (i>m) return -a[0][0];
for (j=1,k=n+1; j<=n; j++)
if (a[j][i]>0 && a[j][0]<a[k][0]) k=j;
pvt(k,i);
}
}
int main(){
scanf("%d%d",&n,&m); int i,j,x,y;
for (i=1; i<=n; i++) scanf("%d",&a[i][0]);
for (i=1; i<=m; i++){
scanf("%d%d%d",&x,&y,&a[0][i]);
for (j=x; j<=y; j++) a[j][i]=1;
}
printf("%d\n",solve());
return 0;
}
by lych
2016.5.27