志愿者招募
申奥成功后,布布经过不懈努力,终于成为奥组委下属公司人力资源部门的主管。布布刚上任就遇到了一个难
题:为即将启动的奥运新项目招募一批短期志愿者。经过估算,这个项目需要N 天才能完成,其中第i 天至少需要
Ai 个人。 布布通过了解得知,一共有M 类志愿者可以招募。其中第i 类可以从第Si 天工作到第Ti 天,招募费用
是每人Ci 元。新官上任三把火,为了出色地完成自己的工作,布布希望用尽量少的费用招募足够的志愿者,但这
并不是他的特长!于是布布找到了你,希望你帮他设计一种最优的招募方案。
方案一:单纯形
思路:设第i类人招募了xi个,aij表示第i个人第j天是否工作,工作为1,不工作为0,则应满足对于所有j,
∑
i
n
a
i
j
x
i
>
=
A
j
\sum_i^na_{ij}x_i >= A_j
∑inaijxi>=Aj,需要最小化
z
=
∑
1
n
C
i
x
i
z=\sum_1^nC_ix_i
z=∑1nCixi,根据对偶原理将矩阵转置后变为最大化即可
#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int MAXN=10005,MAXM=1005,INF=1e9;
const double eps=1e-6;
int n,m;
double a[MAXN][MAXM],b[MAXN],c[MAXM],v;
//a为约束矩阵,有n行m列,c最大化z的系数,b为限制
inline void pivot(int l,int e)
{
b[l]/=a[l][e];
for(int j=1;j<=m;j++) if(j!=e) a[l][j]/=a[l][e];
a[l][e]=1/a[l][e];
for(int i=1;i<=n;i++) if(i!=l&&fabs(a[i][e])>0)
{
b[i]-=a[i][e]*b[l];
for(int j=1;j<=m;j++) if(j!=e) a[i][j]-=a[i][e]*a[l][j];
a[i][e]=-a[i][e]*a[l][e];
}
v+=c[e]*b[l];
for(int j=1;j<=m;j++) if(j!=e) c[j]-=c[e]*a[l][j];
c[e]=-c[e]*a[l][e];
}
inline double simplex()
{
while(true)
{
int e=0,l=0;
for(e=1;e<=m;e++) if(c[e]>eps) break;
if(e==m+1) return v;
double mn=INF;
for(int i=1;i<=n;i++)
if(a[i][e]>eps&&mn>b[i]/a[i][e]) mn=b[i]/a[i][e],l=i;
if(mn==INF) return INF;//无穷大
pivot(l,e);
}
}
int main()
{
while(~scanf("%d%d",&m,&n))
{
for(int i=1;i<=m;++i) scanf("%lf",&c[i]);
int s,t;
for(int i=1;i<=n;++i)
{
scanf("%d%d",&s,&t);
for(int j = 1;j < s;++j) a[i][j] = 0;
for(int j = s;j <= t;++j) a[i][j] = 1;
for(int j = t+1;j <= m;++j) a[i][j] = 0;
scanf("%lf",&b[i]);
}
printf("%.0f",simplex());
}
}
方案二:最小费用流
https://blog.csdn.net/xing_mo/article/details/103986110