解题思路:
如果不清楚单纯形法的可以看这里:
http://blog.sina.com.cn/s/blog_4b1046f80101m9ou.html
题目很像标准单纯形形式但求的是最大值,且所有限制条件都是大于等于。
根据对偶原理,我们只需要把单纯形法整个ABC矩阵转置再求单纯形即是答案。
这是1061的代码,3265的改一下输入即可。
#include<bits/stdc++.h>
using namespace std;
int getint()
{
int i=0,f=1;char c;
for(c=getchar();(c<'0'||c>'9')&&c!='-';c=getchar());
if(c=='-')f=-1,c=getchar();
for(;c>='0'&&c<='9';c=getchar())i=(i<<3)+(i<<1)+c-'0';
return i*f;
}
const int N=1005,M=10005;
const double eps=1e-8;
int n,m;
double ans,a[M][N],b[M],c[N];
void pivot(int l,int e)
{
b[l]/=a[l][e];
for(int i=1;i<=n;i++)
if(i!=e)a[l][i]/=a[l][e];
a[l][e]=1.0/a[l][e];
for(int i=1;i<=m;i++)
if(i!=l&&abs(a[i][e])>eps)
{
b[i]-=a[i][e]*b[l];
for(int j=1;j<=n;j++)
if(j!=e)a[i][j]-=a[i][e]*a[l][j];
a[i][e]=-a[i][e]*a[l][e];
}
ans+=c[e]*b[l];
for(int i=1;i<=n;i++)
if(i!=e)c[i]-=c[e]*a[l][i];
c[e]=-c[e]*a[l][e];
}
void simplex()
{
int i,e,l;
while(true)
{
for(i=1;i<=n;i++)
if(c[i]>eps){e=i;break;}
if(i==n+1)return;
double t=1e20;
for(i=1;i<=m;i++)
if(a[i][e]>eps&&b[i]/a[i][e]<t)
l=i,t=b[i]/a[i][e];
if(t==1e20){ans=t;return;}
pivot(l,e);
}
}
int main()
{
//freopen("lx.in","r",stdin);
n=getint(),m=getint();
for(int i=1;i<=n;i++)c[i]=getint();
for(int i=1;i<=m;i++)
{
int l=getint(),r=getint();b[i]=getint();
for(int j=l;j<=r;j++)a[i][j]++;
}
simplex();
printf("%d",(int)(ans+0.5));
return 0;
}