BZOJ 3265 志愿者招募(BZOJ 1061)加强版 && BZOJ 3112 [Zjoi2013]防守战线 单纯形

BZOJ 3265
题目大意：现要进行n天活动，第i天需要ai个志愿者。有m种志愿者，每种志愿者工作时间为若干个区间[li,ri]（BZOJ 1061：一个区间），招募一个费用为ci，求最小花费

m个变量n个限制条件求函数最值，单纯形裸上就好…
http://wenku.baidu.com/view/ce5784754a7302768f99391d

应该需要用一下对偶原理

#include <cstdio>
#include <algorithm>
#include <cmath>
#define N 1005
#define M 10005
#define eps 1e-7
#define INF 1e10
using namespace std;
int n,m;
double a[M][N],b[M],c[N],v;
void Pivot(int l,int e) {
    double tmp=a[l][e];
    a[l][e]=1;
    for(int i=1;i<=n;++i) a[l][i]/=tmp;
    b[l]/=tmp;
    for(int i=1;i<=m;++i) {
        if(i==l || fabs(a[i][e])<=eps) continue;
        for(int j=1;j<=n;++j)
            if(j!=e) a[i][j]-=a[i][e]*a[l][j];
        b[i]-=a[i][e]*b[l];
        a[i][e]*=-a[l][e];
    }
    for(int i=1;i<=n;++i)
        if(i!=e) c[i]-=c[e]*a[l][i];
    v+=c[e]*b[l];
    c[e]*=-a[l][e];
    return ;
}
double Simplex() {
    while(1) {
        int l,e;
        for(e=1;e<=n;++e)
            if(c[e]>eps) break;
        if(e==n+1) return v;
        double tmp=INF;
        for(int i=1;i<=m;++i)
            if(a[i][e]>eps && b[i]/a[i][e]<tmp)
                tmp=b[i]/a[i][e], l=i;
        if(tmp==INF) return INF;
        Pivot(l,e);
    }
    return -1;
}
int main() {
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;++i) scanf("%lf",c+i);
    for(int i=1,tmp,x,y,z;i<=m;++i) {
        scanf("%d",&tmp);
        while(tmp--) {
            scanf("%d%d",&x,&y);
            for(int j=x;j<=y;++j) a[i][j]=1;
        }
        scanf("%d",&z);
        b[i]=z;
    }
    printf("%d\n",int(Simplex()+0.5));
    return 0;
}

BZOJ 3112
题目大意：现有战场可以看做使一个长度为n的序列，现在需要在这个序列上建塔来防守敌兵，在序列第i 号位置上建一座塔有Ci 的花费，且一个位置可以建任意多的塔，费用累加计算。有m 个区间[L1, R1], [L2, R2], …, [Lm, Rm]，在第i 个区间的范围内要建至少Di 座塔。求最少花费。

和上面那道差不多…

#include <cstdio>
#include <cmath>
#define N 1005
#define eps 1e-7
#define INF 1e10
using namespace std;
int n,m;
double a[N][N*10],b[N],c[N*10],v;
void Pivot(int l,int e) {
    double tmp=a[l][e];
    a[l][e]=1;
    for(int i=1;i<=m;++i) a[l][i]/=tmp;
    b[l]/=tmp;
    for(int i=1;i<=n;++i) {
        if(i==l || fabs(a[i][e])<eps) continue;
        for(int j=1;j<=m;++j)
            if(j!=e) a[i][j]-=a[i][e]*a[l][j];
        b[i]-=b[l]*a[i][e];
        a[i][e]*=-a[l][e];
    }
    for(int i=1;i<=m;++i)
        if(i!=e) c[i]-=c[e]*a[l][i];
    v+=b[l]*c[e];
    c[e]*=-a[l][e];
    return ;
}
double Simplex() {
    while(1) {
        int l,e;
        for(e=1;e<=m;++e)
            if(c[e]>eps) break;
        if(e==m+1) return v;
        double tmp=INF;
        for(int i=1;i<=n;++i)
            if(a[i][e]>eps && b[i]/a[i][e]<tmp)
                tmp=b[i]/a[i][e], l=i;
        if(tmp==INF) return INF;
        Pivot(l,e);
    }
    return -1;
}
int main() {
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;++i) scanf("%lf",b+i);
    for(int i=1,x,y,z;i<=m;++i) {
        scanf("%d%d%d",&x,&y,&z);
        for(int j=x;j<=y;++j) a[j][i]=1;
        c[i]=z;
    }
    printf("%d\n",int(Simplex()+0.5));
    return 0;
}