bzoj 3112: [Zjoi2013]防守战线 单纯形

        单纯形转对偶图。。

       由于是全幺模，直接用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