Симплекс метод

#include <iostream.h>
#include <fstream.h>
#include <iomanip.h>
#include <string.h>
using namespace std;
 
fstream in( "lptask.txt", ios_base::in );
 
double tbl[4][8]={
{3, 11, 0, 0, 0, 0, 0, 0},
{-3, 8, 1, 0, 0, 30, 0, 2},
{2, 7, 0, 1, 0, 62, 0, 3},
{1, 6, 0, 0, 1, 98, 0, 4}
};
int m=8,n=4, col, row;
 
 
void input_tbl() { // çàïîëíåíèå ñèìïëåêñ-òàáëèöû ÷òåíèåì ââîäà
int i=0, j=0;
while(true) {
string x;
in >> x;
if( x == "/" ) {
if( j==0 ) break;
if( i==0 ) m = j + 1;
else {
int j1; in >> j1;
tbl[i][m+1] = j1;
}
tbl[i][m] = 0;
j = 0;
++i;
continue;
}
tbl[i][j] = atof(x.c_str());
++j;
}
n = i;
tbl[0][m-1] = 0;
tbl[0][m+1] = 0;
}
 
void output_tbl() { // âûâîäèò òàáëèöó ñèìïëåêñ-ìåòîäà â óäîáíîì äëÿ ÷òåíèÿ âèäå
for( int i=0; i<n; ++i ) {
for( int j=0; j<=m+1; ++j ) {
cout << setw(9) << tbl[i][j];
}
cout << "\n";
}
}
 
void choose_col_row() { // îïðåäåëÿåò âåäóùèé ñòîëáåö è âåäóùóþ ñòðîêó
 
col = 0;
bool flag = (tbl[0][col] < 0);
for( int j=0; j<m; ++j ) {
double ci = tbl[0][j];
if( ci<0 && ci<tbl[0][col] ) {
flag = true;
col = j;
}
}
if( !flag ) {
col = -1;
row = -1;
return;
}
 
for( int i=1; i<n; ++i ) {
if( tbl[i][m] != 0 ) continue;
double den = tbl[i][col];
if( den == 0 )
tbl[i][m] = 1.e100;
else
tbl[i][m] = tbl[i][m-1]/den;
}
row = 1;
for( int i=1; i<n; ++i ) {
if( tbl[i][m] < tbl[row][m] ) row = i;
}
if( tbl[row][m] >= 1.e98 ) row = -1;
}
 
void transform() {
double den = tbl[row][col];
for( int j=0; j<m; ++j ) {
tbl[row][j] /= den;
}
for( int i=0; i<n; ++i ) {
if( i == row ) continue;
double k = tbl[i][col];
for( int j=0; j<m; ++j ) {
tbl[i][j] -= k*tbl[row][j];
}
}
tbl[row][m+1] = col;
 
for( int i=1; i<=n; ++i ) {
tbl[i][m] = 0;
}
tbl[row][m] = 1e100;
}
 
bool solve() { // ðåøàåò çàäà÷ó - âûïîëíÿåò íóæíî êîëè÷åñòâî èòåðàöèé
while(1) {
cout << "\n";
 
choose_col_row();
 
output_tbl();
cout << "col = " << col << "\n";
cout << "row = " << row << "\n";
 
if( col == -1 ) {
cout << "\nnajdeno optimalnoe reshenie\n";
return true;
}
if( row == -1 ) {
cout << "\nluchevoe reshenie\n";
return false;
}
 
transform();
}
}
 
int main(void) { // ãëàâíàÿ òî÷êà âõîäà
 
input_tbl();
 
bool flag = solve();
 
if( flag ) {
for( int i=1; i<n; ++i ) {
cout << tbl[i][m+1] << " = " << tbl[i][m-1] << "\n";
}
}
}

tbl[i][j] = atof(x.c_str());

//#include "stdafx.h"
#include <iostream.h>
#include <stdlib.h>
#include <windows.h> 
#include <fstream.h>
 
using namespace std;
 
const int str=3;
const int stb=5;
float mas[str][stb];
float R_mas[str][stb];
bool flag=false;
int minI, minJ=1;
 
//Ôóíêöèè
int min_J(float temp[str][stb], int minI)
{
int num;
        if      ((temp[str-1][minJ]/temp[str-1][1]) <= (temp[str][minJ]/temp[str][1]) )
                num=1;
        else
                num=0;
        return num;
}//??
 
int min_I(float temp[str][stb])
{
        int min[2]={};
        for(int i=0; i<str; i++)
        if (temp[i][stb-1]<min[0])
        {
                min[0]=int(temp[0][i]);
                min[1]=i;
        }
        return min[1];
}
 
 
int main()
{
//×òåíèå èç ôàéëà
        //FILE *myFile;
       // myFile=fopen("3.txt", "rt"); 
        char temp[4];
        for (int i=0; i<str; i++)
                for (int j=0; j<stb; j++)
                {
                        //fgets(temp, 4, myFile);
                        mas[i][j]=float(atoi(temp));
                }
 
//îòðèöàòåëüíûå 1,2 è 5 ñòîëáöû
for (int i=0; i<str; i++)
        for (int j=0; j<stb; j++)
                if ((j != 2) && (j != 3) )
                        mas[i][j]*=-1;
 
//Âûâîä 
        for (int i=0; i<str; i++)
        {
                cout<<endl;
                for (int j=0; j<stb; j++)
                        cout<<mas[i][j]<<"\t";
        }
 
        cout<<endl<<endl;
        int min[2]={};
        for (int j=1; j<stb; j++)
                
                
 
        while(!flag)
                        
//Ïîèñê îòðèöàòåëüíûõ ÷èñåë
                for (int j=1; j<stb; j++)
                        for(int i=0; i<stb-1;i++)
                {                       
                        if(mas[i][j]<0)
                        {
                                flag=true; //ôëàã
                                
//Ïîèñê ìèíèìàëüíîé ñòðîêè                              
                                minI=min_I(mas);
 
//Ïîèñê ìèíèìàëüíîãî ñòîëáöà
                                minJ=min_J(mas, minI);
 
//Îñíîâíûå âû÷èñëåíèÿ
                                for(int i=0; i<stb;i++)
                                        R_mas[minI][i]=mas[minI][i]/mas[minI][minJ];            
                                for (int i=0; i<str;i++)
                                        for(int j=0; j<stb;j++)
                                                if (i != minI)
                                                        R_mas[i][j]=mas[i][j]-(R_mas[minI][j] * mas[i][minJ]);
                    }
                }
       
                        
 
//Âûâîä ìàññèâà
                for (int i=0; i<str; i++)
                {
                        for (int j=0; j<stb; j++)
                cout<<R_mas[i][j]<<"\t";
                cout<<endl;
                        
                }
 
        }