输入初始单纯形表后的单纯形程序，线性规划 未写完

虽然学了单纯形，但学得很差，终于看着书把最简单的例子给写出来了，对不起白老师啊（教我们运筹学的老师）。

该程序执行的初试条件是自己把初试单纯形表搞出来，本程序只负责进行迭代，（偷笑）。

唉，希望有人能得到帮助，本来做这个题是为了去解一个线性规划的题的，但感觉有点麻烦，还是先歇着吧，看以后变强了再来A这个实数幻方的题

该题链接：http://122.207.68.93/OnlineJudge/problem.php?cid=2031&pid=9

Problem J: Real Numbers Magic Square

咳咳，以后再说吧

贴代码：

//这里在找进基变量和出基变量时有时会出现死循环，尽管概率很低
//可以使用勃兰特规则，该程序没有写，勃兰特规则自己百度，或参考运筹学教材,清华大学出版社出版，36页,退化那一块
#include <cstdio>
#define MAXN 100
#define INF 0x7fffffff
int ind ,outd;
double a[MAXN][MAXN];
double cb[MAXN],b[MAXN];
int xb[MAXN],nxb[MAXN];
double theta[MAXN],c[MAXN];
int n,m;
void init()
{
    for(int i=1; i<=m; ++i)
        for(int j=1; j<=n; ++j)
            scanf("%lf",&a[i][j]);
    for(int i=1; i<=n; ++i)
        scanf("%lf",&c[i]);
    for(int i=1; i<=m; ++i)
        scanf("%lf%d%lf",&cb[i],&xb[i],&b[i]);
    for(int i=1; i<=n-m; ++i)
        scanf("%d",&nxb[i]);
}
void print()
{
    for(int i=1; i<=m; ++i)
    {
        for(int j=1; j<=n; ++j)
        {
            printf("a[%d][%d] = %lf ",i,j,a[i][j]);
        }
        puts("");
    }
//    for(int i=1; i<=n; ++i)
//    {
//        printf("c[%d] = %lf ",i,c[i]);
//    }
//    puts("");
    for(int i=1; i<=m; ++i)
    {
        printf("cb[%d] = %lf xb[%d] =%d b[%d] = %lf\n",i,cb[i],i,xb[i],i,b[i]);
    }
    for(int i=1; i<=n-m; ++i)
    {
        printf("nxb[%d] = %d\n",i,nxb[i]);
    }
}
int findInVarible()//找进基变量
{
    double max = 0;
    ind = -1;
    for(int i=1; i<=n-m; ++i)
    {
        int j = nxb[i];
        double sig = c[j];
        for(int k =1; k<=m; ++k)
        {
            sig -= cb[k]*a[k][j];
        }
        bool flag = true;
        if(sig > 0)
        {
            for(int k=1; k<=m; ++k)
            {
                if(a[k][j] > 0)
                {
                    flag = false;
                    break;
                }
            }
            if(flag) return -2;//出现无界解
            if(sig > max)
            {
                printf("sig = %lf\n",sig);
                max = sig;
                ind  = i;
            }
        }
    }
    return ind;
}
int findOutVarible()//找出基变量
{
    double min = INF;
    outd = -1;
    int j = nxb[ind];
    for(int k=1; k<=m; ++k)
    {
        if(a[k][j] > 0 && b[k]/a[k][j] < min)
        {
            min = b[k]/a[k][j] ;
            outd = k;
            printf("theta = %lf\n",min);
        }
    }
    return outd;
}
void Guass(int x,int y)//高斯变换化进基变量所在列为单位向量
{
    double de = a[y][x];
    for(int j = 1; j<=n; ++j)
        a[y][j] = a[y][j]/de;
    b[y] = b[y]/de;
    for(int i=1; i<=m; ++i)
    {
        if(i == y) continue;
        de = a[i][x];
        for(int j=1; j<=n; ++j)
            a[i][j] =a[i][j]-a[y][j]*de;
        b[i] =b[i]-b[y]*de;
    }
}
int main()
{
//    freopen("in.cpp","r",stdin);
//    freopen("out.cpp","w",stdout);
    while(~scanf("%d%d",&m,&n))
    {
        init();//建立初始单纯形表，这里为了简化，直接输入的就是已经初试化好了的单纯形表
        print();
        int ind,outd;//进基变量和出基变量的序号
        while(1)
        {
            ind = findInVarible();
            if(ind == -2 || ind == -1) break;
            outd = findOutVarible();
            printf("ind = %d outd = %d\n",nxb[ind],xb[outd]);
            Guass(nxb[ind],outd);
            int tmp = nxb[ind];
            nxb[ind] = xb[outd];
            xb[outd] = tmp;
            cb[outd] = c[tmp];
            print();
        }
        if(ind == -2) printf("No solution\n");
        else
        {
            double sum = 0;
            for(int i=1; i<=m; ++i)
                sum += cb[i]*b[i];
            printf("sum = %lf\n",sum);//输出最后的解
        }
    }
    return 0;
}