bzoj1061题解

【解题思路】

　　设类型i的志愿者，即第S_i天~第T_i天工作的志愿者，共招募x_i个，于是有不等式组Σx_j≥A_i(Sj≤i≤Tj)。

　　这样，题目就变成了求一组满足一次不等式组的x_i，使ΣC_ix_i最小，即标准的线性规划形式。

　　本人比较懒。。并不想建图跑费用流之类的。。于是写了单纯形。。复杂度O(松)。

【参考程序】

 #pragma GCC optimize(2)
 #include <cmath>
 #include <cstdio>
 #include <cstring>
 #define REP(i,low,high) for(register int i=(low);i<=(high);i++)
 #define INF 1e10
 #define eps 1e-7
 using namespace std;
 inline bool getmin(double &tar,const double &pat) {return pat+eps<tar?tar=pat,:;}
 double a[][]; int m,n;
 inline int check() {REP(i,,m) if(a[i][]>eps) return i; return ;}
 inline double Simplex()
 {
     while(int x=check())
     {
         int y=; double lim=INF; REP(i,,n) if(a[x][i]>eps&&getmin(lim,a[][i]/a[x][i])) y=i;
         if(!y) return a[][]=INF; double p=1.0/a[x][y]; a[x][y]=1.0; REP(i,,m) a[i][y]*=p;
         REP(i,,n) if(i!=y)
         {
             double now=a[x][i];
             if(fabs(now)>eps) {a[][i]-=now*a[][y]; REP(j,,m) a[j][i]-=now*a[j][y]; a[x][i]=-now*p;}
         }
         double now=a[x][]; a[][]+=now*a[][y]; REP(i,,m) a[i][]-=now*a[i][y]; a[x][]=-now*p;
     }
     return a[][];
 }
 int main()
 {
     scanf("%d%d",&m,&n); REP(i,,m) scanf("%lf",&a[i][]);
     REP(i,,n) {int S,T; scanf("%d%d",&S,&T); REP(j,S,T) a[j][i]=1.0; scanf("%lf",&a[][i]);}
     return printf("%.0lf\n",Simplex()+eps),;
 }