[BZOJ3265]志愿者招募加强版

原题地址

费用流神题，单纯形裸题(指原NOI题)…
这个题一不对偶就T&WA，也不知道是为什么OLZ…
对偶大法好…然而我不会证…

AC code:

#include <cstdio>
#include <cmath>
typedef double llf;
const int N=10010;
const int M=1010;
const int INF=1<<29;
const llf eps=1e-7;
int n,m;
int idb[N],idn[M];
llf v;
llf b[N],c[M],tc[M];
llf A[N][M];

void pivot(int l,int e,int type){
    b[l]/=A[l][e];
    for(int i=type;i<=m;i++) if(i!=e) A[l][i]/=A[l][e];
    A[l][e]=1/A[l][e];
    for(int i=type;i<=n;i++){
        if(i==l||fabs(A[i][e])<eps) continue;
        b[i]-=A[i][e]*b[l];
        for(int j=type;j<=m;j++){
            if(j==e) continue;
            A[i][j]-=A[i][e]*A[l][j];
        }
        A[i][e]*=-A[l][e];
    }
    v+=c[e]*b[l];
    for(int i=type;i<=m;i++){
        if(i==e) continue;
        c[i]-=c[e]*A[l][i];
    }
    c[e]*=-A[l][e];
    if(!type){
        int tmp=idb[l];idb[l]=idn[e];idn[e]=tmp;
    }
}

void opt(int type){
    while(1){
        int e=INF,l=INF;
        llf tmp=INF;
        for(int i=type;i<=m;i++){
            if(c[i]>eps){
                e=i;
                break;
            }
        }
        if(e==INF) return ;
        for(int i=type;i<=n;i++){
            if(A[i][e]>eps&&b[i]/A[i][e]<tmp){
                l=i;
                tmp=b[i]/A[i][e];
            }
        }
        if(l==INF){
            v=INF;
            return ;
        }
        pivot(l,e,type);
    }
}

bool init(){
    int l=0;
    for(int i=1;i<=n;i++) if((!l)||b[i]<b[l]) l=i;
    if(b[l]>=eps) return 1;
    for(int i=1;i<=n;i++) A[i][0]=-1;
    for(int i=1;i<=m;i++) tc[i]=c[i];
    c[0]=-1;
    for(int i=1;i<=m;i++) c[i]=0;
    pivot(l,0,0);
    opt(0);
    if(fabs(v)>eps) return 0;
    for(int i=1;i<=n;i++){
        if(!idb[i]){
            pivot(i,0,0);
            break;
        }
    }
    for(int i=1;i<=m;i++){
        if(!idn[i]){
            for(int j=i;j;j--) idn[j]=idn[j-1];
            for(int j=1;j<=n;j++)
                for(int k=i;k;k--)
                    A[j][k]=A[j][k-1];
            break;
        }
    }
    for(int i=1;i<=m;i++) c[i]=0;
    for(int i=1;i<=n;i++){
        if(idb[i]<=m){
            v+=tc[idb[i]]*b[i];
            for(int j=1;j<=m;j++) c[j]-=tc[idb[i]]*A[i][j];
        }
    }
    for(int i=1;i<=m;i++) if(idn[i]<=m) c[i]+=tc[idn[i]];
    return 1;
}

void simplex(){
    init();
    opt(1);
    printf("%d\n",(int)(v+0.5));
}

int main(){
    scanf("%d%d",&m,&n);
    for(int i=1;i<=n;i++) idb[i]=i+m;
    for(int i=1;i<=m;i++) idn[i]=i;
    for(int i=1;i<=m;i++) scanf("%lf",&c[i]);
    for(int i=1;i<=n;i++){
        int k,s,t;
        scanf("%d",&k);
        for(int j=1;j<=k;j++){
            scanf("%d%d",&s,&t);
            for(int jj=s;jj<=t;jj++) A[i][jj]=1;
        }
        scanf("%lf",&b[i]);
    }
    simplex();

    return 0;
}