單純形模板

板子來自神犇zzq

 

  1 #include <ctime>
  2 #include <cmath>
  3 #include <cstdio>
  4 #include <cstdlib>
  5 #include <cstring>
  6 #include <algorithm>
  7 #define N1 111
  8 #define dd double
  9 #define ld long double
 10 #define ll long long 
 11 #define ull unsigned long long 
 12 using namespace std;
 13 
 14 const ld eps=1e-9;
 15 const int inf=0x3f3f3f3f;
 16 int gint()
 17 {
 18     int ret=0,fh=1;char c=getchar();
 19     while(c<'0'||c>'9'){if(c=='-')fh=-1;c=getchar();}
 20     while(c>='0'&&c<='9'){ret=ret*10+c-'0';c=getchar();}
 21     return ret*fh;
 22 }
 23 
 24 int n,m,T;
 25 int id[N1];
 26 ld a[N1][N1],val[N1];
 27 
 28 int dcmp(ld x)
 29 {
 30     if(x<-eps) return -1;
 31     if(x>eps) return 1;
 32     return 0;
 33 }
 34 void pivot(int r,int c)
 35 {
 36     swap(id[r+n],id[c]);
 37     int i,j;
 38     ld x=-a[r][c]; a[r][c]=-1;
 39     for(i=0;i<=n;i++) a[r][i]/=x; 
 40     for(i=0;i<=m;i++)
 41     {
 42         if(!dcmp(a[i][c])||i==r) continue;
 43         x=a[i][c]; a[i][c]=0;
 44         for(j=0;j<=n;j++) a[i][j]+=x*a[r][j];
 45     }
 46 }
 47 void solve()
 48 {
 49     int i,j,x,y,fl;
 50     for(i=1;i<=n;i++) id[i]=i;
 51     while(1)
 52     {
 53         x=0,y=0;
 54         for(i=1;i<=m;i++)
 55             if( dcmp(a[i][0])<0 && (!x||rand()&1) ) x=i;
 56         if(!x) break;
 57         for(i=1;i<=n;i++)
 58             if( dcmp(a[x][i])>0 && (!y||rand()&1) ) y=i;
 59         if(!y){ puts("Infeasible"); return; } 
 60         pivot(x,y);
 61     }
 62     ld w,t;
 63     while(1)
 64     {
 65         x=0,y=0;
 66         for(i=1;i<=n;i++)
 67             if( dcmp(a[0][i])>0 && (!x||rand()&1) ) x=i;
 68         if(!x) break;
 69         w=1e18; t=0; 
 70         for(i=1;i<=m;i++)
 71         {
 72             if( dcmp(a[i][x])>=0 ) continue;
 73             t=-a[i][0]/a[i][x];
 74             if(t<w||!y) y=i, w=t;
 75         }
 76         if(!y){ puts("Unbounded"); return; } 
 77         pivot(y,x);
 78     }
 79     printf("%.8Lf\n",a[0][0]);
 80     if(!T) return;
 81     for(i=1;i<=n;i++) val[id[i]]=0;
 82     for(i=n+1;i<=n+m;i++) val[id[i]]=a[i-n][0];
 83     for(i=1;i<=n;i++) printf("%.8Lf ",val[i]);
 84     puts("");
 85 }
 86 
 87 int main()
 88 {
 89     srand(time(NULL));
 90     int i,j;
 91     scanf("%d%d%d",&n,&m,&T);
 92     for(i=1;i<=n;i++) scanf("%Lf",&a[0][i]);
 93     for(i=1;i<=m;i++)
 94     {
 95         for(j=1;j<=n;j++) scanf("%Lf",&a[i][j]), a[i][j]*=-1;
 96         scanf("%Lf",&a[i][0]);
 97     }
 98     solve();
 99     return 0;
100 }