【BZOJ3112】[Zjoi2013]防守战线 单纯形法

【BZOJ3112】[Zjoi2013]防守战线

题解：依旧是转化成对偶问题，然后敲板子就行了~

建完表后发现跟志愿者招募的表正好是相反的，感觉很神奇~

#include <cstdio>
#include <cstring>
#include <iostream>
#include <cmath>
using namespace std;
const double inf=1e12;
const double eps=1e-6;
double A[1010][10010];
int p[10010];
int n,m;
int rd()
{
	int ret=0,f=1;	char gc=getchar();
	while(gc<'0'||gc>'9')	{if(gc=='-')f=-f;	gc=getchar();}
	while(gc>='0'&&gc<='9')	ret=ret*10+gc-'0',gc=getchar();
	return ret*f;
}
void pivot(int l,int e)
{
	double t=A[l][e];A[l][e]=1;
	int i,j;
	for(j=0;j<=n;j++)	A[l][j]/=t;
	for(p[0]=j=0;j<=n;j++)	if(fabs(A[l][j])>eps)	p[++p[0]]=j;
	for(i=0;i<=m;i++)	if(i!=l&&fabs(A[i][e])>eps)
	{
		t=A[i][e],A[i][e]=0;
		for(j=1;j<=p[0];j++)	A[i][p[j]]-=t*A[l][p[j]];
	}
}
void simplex()
{
	while(1)
	{
		double mn=inf;
		int i,l=0,e=0;
		for(i=1;i<=n;i++)	if(A[0][i]>eps)
		{
			e=i;
			break;
		}
		if(!e)	return;
		for(i=1;i<=m;i++)	if(A[i][e]>eps&&A[i][0]/A[i][e]<mn)	mn=A[i][0]/A[i][e],l=i;
		if(!l)	return;
		pivot(l,e);
	}
}
int main()
{
	m=rd(),n=rd();
	int i,j,l,r;
	for(i=1;i<=m;i++)	A[i][0]=rd();
	for(i=1;i<=n;i++)
	{
		l=rd(),r=rd(),A[0][i]=rd();
		for(j=l;j<=r;j++)	A[j][i]=1;
	}
	simplex();
	printf("%d",int(-A[0][0]+0.5));
	return 0;
}