uoj179 线性规划（单纯形模板题）

今天的一天：《学习单纯形解线性规划—从入门到放弃》
放弃.jpg

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
#define eps 1e-7
#define N 25
#define inf 0x3f3f3f3f
inline char gc(){
    static char buf[1<<16],*S,*T;
    if(S==T){T=(S=buf)+fread(buf,1,1<<16,stdin);if(T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=gc();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
    return x*f;
}
int n,m,type;
double a[N][N],ans[N];
int id[N<<1];
int q[N];
inline void pivot(int l,int e){
    swap(id[n+l],id[e]);
    double t=a[l][e];a[l][e]=1;
    for(int j=0;j<=n;++j) a[l][j]/=t;
    for(int i=0;i<=m;++i){
        if(i==l||abs(a[i][e])<eps) continue;
        t=a[i][e];a[i][e]=0;
        for(int j=0;j<=n;++j) a[i][j]-=a[l][j]*t;
    }
}
inline bool init(){
    while(1){
        int l=0,e=0;
        for(int i=1;i<=m;++i) if(a[i][0]<-eps&&(!l||(rand()&1))) l=i;
        if(!l) return 1;
        for(int j=1;j<=n;++j) if(a[l][j]<-eps&&(!e||(rand()&1))) e=j;
        if(!e){puts("Infeasible");return 0;}
        pivot(l,e);
    }
}
inline bool simplex(){
    while(1){
        int l=0,e=0;double mn=inf;
        for(int j=1;j<=n;++j) if(a[0][j]>eps){e=j;break;}
        if(!e) return 1;
        for(int i=1;i<=m;++i)
            if(a[i][e]>eps&&a[i][0]/a[i][e]<mn) mn=a[i][0]/a[i][e],l=i;
        if(!l){puts("Unbounded");return 0;}
        pivot(l,e);
    } 
}
int main(){
//  freopen("a.in","r",stdin);
    srand(20000712);
    n=read();m=read();type=read();
    for(int i=1;i<=n;++i) a[0][i]=read();
    for(int i=1;i<=m;++i){
        for(int j=1;j<=n;++j) a[i][j]=read();
        a[i][0]=read();
    }for(int i=1;i<=n;++i) id[i]=i;
    if(init()&&simplex()){
        printf("%.8lf\n",-a[0][0]);
        if(type){
            for(int i=1;i<=m;++i) ans[id[n+i]]=a[i][0];
            for(int i=1;i<=n;++i) printf("%.8lf ",ans[i]);
        }return 0;
    }
}