单纯形法模板

#include <bits/stdc++.h>
#define eps 1e-8
#define INF 1e100
using namespace std;
long double a[105][105],x[205],ans,c[105];
int n,m,xx,yy,t,B[105],N[105],d[105],tt;
void pivot(int l,int e)
{
    swap(B[l],N[e]);
    for(int i=0;i<=n;i++)if(i!=e)a[l][i]/=a[l][e];
    a[l][e]=(long double)1/a[l][e];
    for(int i=0;i<=m;i++)if(i!=l&&fabs(a[i][e])>eps)
    {
        for(int j=0;j<=n;j++)if(j!=e)a[i][j]-=a[i][e]*a[l][j];
        a[i][e]=-a[i][e]*a[l][e];
    }
}
int simplex()
{
    int l,e;
    while(1)
    {
        e=n+1; tt=0;
        for(int i=1;i<=n;i++)if(a[0][i]>eps)d[++tt]=i;
        if(tt==0)break; e=d[rand()%tt+1]; long double temp=INF;
        for(int i=1;i<=m;i++)if((a[i][e]>eps)and(a[i][0]/a[i][e]<temp)){ l=i; temp=a[i][0]/a[i][e]; }
        if(temp==INF){ printf("Unbounded
"); return 0; }
        pivot(l,e);
    }
    ans=-a[0][0]; return 1;
}
int main()
{
    scanf("%d%d%d
",&n,&m,&t);
    for(int i=1;i<=n;i++)scanf("%Lf",&c[i]); a[0][n+1]=-1;
    for(int i=1;i<=n;i++)N[i]=i; N[n+1]=0;
    for(int i=1;i<=m;i++)B[i]=n+i;
    for(int i=1;i<=m;i++)
    {
        for(int j=1;j<=n;j++)scanf("%Lf",&a[i][j]); a[i][n+1]=-1;
        scanf("%Lf",&a[i][0]);
    }
    n++; int mm=1; for(int i=2;i<=m;i++)if(a[i][0]<a[mm][0])mm=i; 
    if(a[mm][0]<-eps)pivot(mm,n);
    if(!simplex())return 0;
    if(ans<-eps){ printf("Infeasible
"); return 0; }
    xx=0; for(int i=1;i<=n;i++)if(N[i]==0){ xx=i; break; }
    if(xx==0)
    {
        yy=0;
        for(int i=1;i<=m;i++)if(B[i]==0){ yy=i; break; }
        for(xx=0;xx<=n;xx++)if(abs(a[yy][xx])>eps)break; 
        if(xx<=n)pivot(yy,xx);
    }
    if(xx<=n)
    {
        for(int i=xx;i<=n-1;i++)N[i]=N[i+1]; N[n]=0;
        for(int i=1;i<=m;i++)for(int j=xx;j<=n-1;j++)a[i][j]=a[i][j+1];
        n--; for(int i=0;i<=n;i++)a[0][i]=0;
        for(int i=1;i<=n;i++)
        {
            xx=0; for(int j=1;j<=n;j++)if(N[j]==i){ xx=j; break; }
            if(xx>0)a[0][xx]+=c[i];else
            {
                for(int j=1;j<=m;j++)if(B[j]==i){ xx=j; break; }
                for(int j=0;j<=n;j++)a[0][j]=a[0][j]-c[i]*a[xx][j];
            }
        }
    }else
    {
        for(int i=yy;i<=m-1;i++)for(int j=0;j<=n;j++)a[i][j]=a[i+1][j];
        m--;
    }
    if(!simplex())return 0;
    printf("%.10Lf
",ans);
    if(t==1)
    {
        for(int i=1;i<=m;i++)x[B[i]]=a[i][0];
        for(int i=1;i<=n;i++)printf("%.10Lf ",x[i]);
        printf("
");
    }
    return 0;
}