模板：线性规划

//数组开不下时勿用
const int maxn=110;
const double eps=1e-10;
double a[maxn][maxn];
double b[maxn],c[maxn],v;
int N[maxn],B[maxn],n,m;

int SGN(double x){
    return (x>eps)-(x<-eps);
}

void Init(){
    N[0]=B[0]=0;v=0.0;
    for(int i=1;i<=n;i++)N[++N[0]]=i;
    for(int i=1;i<=m;i++)B[++B[0]]=n+i;
}

void Pivot(int l,int e){
    b[e]=b[l]/a[l][e];
    a[e][l]=1.0/a[l][e];
    for(int i=1;i<=N[0];i++)
        if(N[i]!=e)a[e][N[i]]=a[l][N[i]]/a[l][e];
    
    for(int i=1;i<=B[0];i++)
        if(B[i]!=l){
            b[B[i]]-=a[B[i]][e]*b[e];
            a[B[i]][l]=-a[B[i]][e]*a[e][l];
            for(int j=1;j<=N[0];j++)
                if(N[j]!=e)a[B[i]][N[j]]-=a[B[i]][e]*a[e][N[j]];
        }    
        
    v+=b[e]*c[e];
    c[l]=-c[e]*a[e][l];
    for(int i=1;i<=N[0];i++)
        if(N[i]!=e)c[N[i]]-=c[e]*a[e][N[i]];
    for(int i=1;i<=N[0];i++)if(N[i]==e)N[i]=l;
    for(int i=1;i<=B[0];i++)if(B[i]==l)B[i]=e;
}

bool Simplex(){
    while(true){
        double lam=-1;
        int e=maxn,l=maxn;
        for(int i=1;i<=N[0];i++)
            if(SGN(c[N[i]])>0&&e>N[i])
                e=N[i];
        if(e==maxn)break;
        
        for(int i=1;i<=B[0];i++)
            if(SGN(a[B[i]][e])>0){
                double tmp=b[B[i]]/a[B[i]][e];
                if(lam==-1||SGN(tmp-lam)<0||SGN(tmp-lam)==0&&l>B[i])
                    lam=tmp,l=B[i];
            }
        if(l==maxn)return false;
        Pivot(l,e);    
    }
    return true;
}

 　　这个是整数规划的模板，总能对，但正确性不显然。

const int maxr=1010;
const int maxc=10010;

int n,m,nxt[maxc];
int a[maxr][maxc];

void Pivot(int l,int e){
    int pre=maxc-1;
    for(int i=0;i<=n;i++)
        if(a[l][i]!=0){nxt[pre]=i;pre=i;}
    nxt[pre]=-1;
            
    for(int i=0,t;i<=m;i++)
        if(i!=l&&(t=a[i][e])){
            a[i][e]=0;
            for(int j=nxt[maxc-1];j!=-1;j=nxt[j])
                a[i][j]+=t*a[l][j];
        }
}

void Simplex(){
    while(true){
        int e=0,l=0;
        for(int i=1;i<=n;i++)
            if(a[0][i]>0){e=i;break;}
        if(e==0)break;        
        for(int i=1;i<=m;i++)
            if(a[i][e]<0&&(!l||a[l][0]>a[i][0]))
                {l=i;}
        
        Pivot(l,e);    
    }
}