题面:http://img.d1kt.cn/programming/1347_1181.png
题解:对偶单纯形裸题。
代码
#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#define ll long long
#define mod 10000007
using namespace std;
inline int read()
{
int x=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if (ch=='-') f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
return x*f;
}
#define inf 10000000000.0
#define eps (1e-9)
double ans,a[1005][10005],c[10005],b[10005];
int n,m;
inline 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/a[l][e];
for (int i=1;i<=m;i++)if (i!=l&&fabs(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];
}
inline void simplex()
{
int l,e;
while (1)
{
for (e=1;e<=n;e++) if (c[e]>eps) break;
if (e>n) break;
double t=inf;
for (int i=1;i<=m;i++)if (a[i][e]>eps&&b[i]/a[i][e]<t)
{
t=b[i]/a[i][e];
l=i;
}
pivot(l,e);
}
}
int main()
{
n=read();m=read();swap(n,m);
for (int i=1;i<=m;i++) b[i]=read();
for (int i=1;i<=n;i++)
{
int l=read(),r=read();c[i]=read();
for (int j=l;j<=r;j++) a[j][i]++;
}
simplex();
printf("%.0lf",ans+eps);
return 0;
}