C++实现对偶单纯形算法

背景

这学期修了卜东波老师的算法课，要编程实现对偶单纯形算法解决线性规划问题，笔者在网上找了好久，没有看到编程实现对偶算法的。下面借助一个实例，编程实现对偶算法

问题描述

For the problem

代码实现

/*
对偶单纯形算法
输入格式：
A：单纯形表矩阵，各变量的系数和常数 M * N+1 维
a_11 a_12 ··· a_1N b_1
a_21 a_22 ··· a_2N b_2
 ·     ·        ·   ·
 ·     ·        ·   ·
 ·     ·        ·   ·
a_M1 a_M2 ··· a_MN b_M

C：目标函数系数向量 N 维
Lambda：检验数向量 N 维
base_V：基向量下标 M 维
B：基向量在目标函数中对应的系数 M 维
goal：目标值

flygeda 2017.12.16
*/
#include <stdio.h>
#include <iostream>

#define N 6  
#define M 3 //N个变量，M个约束条件

using namespace std;

double A[M][N + 1] = { { 3, -1, 1, -2, 0, 0, -3},
                    { 2, 1, 0, 1, 1, 0, 4 },
                    { -1, 3, 0, -3,