UOJ#179. 线性规划[模板]

传送门 http://uoj.ac/problem/179

震惊，博主竟然还不会线性规划！

单纯形实在学不会啊……背个板子当黑盒用……

学（chao）了NanoApe dalao的板子

 #include<iostream>
 #include<algorithm>
 #include<cstring>
 #include<cstdio>
 #include<cmath>
 using namespace std;
 const double eps=1e-;
 const int mxn=;
 int read(){
     int x=,f=;char ch=getchar();
     while(ch<'' || ch>''){if(ch=='-')f=-;ch=getchar();}
     while(ch>='' && ch<=''){x=x*+ch-'';ch=getchar();}
     return x*f;
 }

 int n,m,t;
 double c[mxn];
 double a[mxn][mxn];
 int idx[mxn],idy[mxn];
 int st[mxn],top=;
 void Pivot(int x,int y){//用idy代换idx
     swap(idy[x],idx[y]);
     double tmp=a[x][y];a[x][y]=/a[x][y];
     int i,j;top=;
     for(i=;i<=n;i++)if(y!=i)a[x][i]/=tmp;
     for(i=;i<=n;i++)if((y!=i) && fabs(a[x][i])>eps)st[++top]=i;
     for(i=;i<=m;i++){
         if((i==x)||(fabs(a[i][y])<eps))continue;
         for(j=;j<=top;j++)a[i][st[j]]-=a[x][st[j]]*a[i][y];
         a[i][y]=-a[i][y]/tmp;
     }
     return;
 }
 void Simplex(){
     int i,j;
     while(){
         int x=,y=;double mn=1e15;
         for(i=;i<=n;i++)if(a[][i]>eps){y=i;break;}
         if(!y)break;
         for(i=;i<=m;i++){
             if(a[i][y]>eps && (a[i][]/a[i][y]<mn)){
                 mn=a[i][]/a[i][y];
                 x=i;
             }
         }
         if(!x){printf("Unbounded\n");exit();}
         Pivot(x,y);
     }
     return;
 }
 int main(){
     int i,j;
     n=read();m=read();t=read();
     for(i=;i<=n;i++)a[][i]=read();//目标函数
     for(i=;i<=m;i++){
         for(j=;j<=n;j++)a[i][j]=read();
         a[i][]=read();
     }
     for(i=;i<=n;i++)idx[i]=i;//基变量
     for(i=;i<=m;i++)idy[i]=i+n;//非基变量
     while(){
         int x=,y=;
         for(i=;i<=m;i++)if(a[i][]<-eps && ((!x)||(rand()&)))x=i;
         if(!x)break;
         for(i=;i<=n;i++)if(a[x][i]<-eps && ((!y)||(rand()&)))y=i;
         if(!y){printf("Infeasible\n");return ;}
         Pivot(x,y);
     }
     Simplex();
     printf("%.8f\n",-a[][]);
     if(!t)return ;
     for(i=;i<=n;i++)a[][i]=;
     for(i=;i<=m;i++)if(idy[i]<=n)a[][idy[i]]=a[i][];
     for(i=;i<=n;i++)printf("%.8f ",a[][i]);puts("");
     return ;
 }