Lettiery Procopio / Simpex

main.cpp 3.76 KB
edit raw blame history 
   1
   2
   3
   4
   5
   6
   7
   8
   9
  10
  11
  12
  13
  14
  15
  16
  17
  18
  19
  20
  21
  22
  23
  24
  25
  26
  27
  28
  29
  30
  31
  32
  33
  34
  35
  36
  37
  38
  39
  40
  41
  42
  43
  44
  45
  46
  47
  48
  49
  50
  51
  52
  53
  54
  55
  56
  57
  58
  59
  60
  61
  62
  63
  64
  65
  66
  67
  68
  69
  70
  71
  72
  73
  74
  75
  76
  77
  78
  79
  80
  81
  82
  83
  84
  85
  86
  87
  88
  89
  90
  91
  92
  93
  94
  95
  96
  97
  98
  99
 100
 101
 102
 103
 104
 105
 106
 107
 108
 109
 110
 111
 112
 113
 114
 115
 116
 117
 118
 119
 120
 121
 122
 123
 124
 125
 126
 127
 128
 129
 130
 131
 132
 133
 134
 135
 136
//
//  main.cpp
//  Simplex
//
//  Created by Lettiery on 02/05/16.
//  Copyright © 2016 Lettiery. All rights reserved.
//

#include <iostream>
#include <vector>
#include "Simplex.h"

using namespace std;

void imprimirModelo(double model, int n_variaveis, int n_rest, vector<vector<double> >);
vector<vector<double> > lerModelo(int * n_restricoes, int * n_variaveis, double * model);
vector<vector<double> > criarModelo(int * n_restricoes, int * n_variaveis, double * model);

int main(int argc, const char * argv[]) {
    int n_restricoes, n_variaveis;
    double model;
    vector<vector<double> > tablo;
    
    char c;
    
    cout << "Deseja inserir um modelo? (Y/N): ";
    cin >> c;
    if(c == 'Y' || c == 'y'){
        tablo = lerModelo(&n_restricoes, &n_variaveis, &model);
    }else{
        tablo = criarModelo(&n_restricoes, &n_variaveis, &model);
    }
    
    //vector<vector<double> > tablo = lerModelo(&n_restricoes, &n_variaveis, &model);
    //vector<vector<double> > tablo = criarModelo(&n_restricoes, &n_variaveis, &model);
    
    
    
    imprimirModelo(model, n_variaveis, n_restricoes, tablo);
    
    Simplex alg (model , n_variaveis, n_restricoes, tablo);
    alg.imprimirTablo();
    
    alg.primal();
    
    return 0;
}

vector<vector<double> > criarModelo(int * n_restricoes, int * n_variaveis, double * model){
    *model = -1.0;// minimizar
    *n_variaveis = 2;
    *n_restricoes = 3;
    
    vector<vector<double> > tablo((*n_restricoes) + 1, vector<double>( (*n_variaveis) + (*n_restricoes) + 1 , 0));
    
    
    tablo[0][0] = 0;// inicia z = 0
    
    tablo[0][1] = 1.12; tablo[0][2] = 1.23;//1.12x1 + 1.23x2
    
    tablo[1][1] = 1; tablo[1][2] = 0; tablo[1][0] = 300;// 3x1 + 1x2+ 0x2 <= 15
    tablo[2][1] = 0; tablo[2][2] = 1; tablo[2][0] = 270;// 2x1 + 2x2 + 2x3 <= 20
    tablo[3][1] = -1; tablo[3][2] = -1; tablo[3][0] = -500;
    
    for(int c = 1; c < tablo[0].size();c++){
        tablo[0][c] *= *model;
    }
    return tablo;
}

vector<vector<double> > lerModelo(int * n_restricoes, int * n_variaveis, double * model){
    
    char max;
    
    printf("Entrada do Modelo\n\n");
    printf("MAXIMIZAR? (Y/N): "); scanf("%c", &max);
    printf("Quantidade do variáveis: "); scanf("%d", n_variaveis);
    printf("Quantidade de restrições: "); scanf("%d", n_restricoes);
    
    vector<vector<double> > tablo((*n_restricoes) + 1, vector<double>( (*n_variaveis) + (*n_restricoes) + 1 , 0));
    
    if (max == 'Y' || max=='y'){
        *model = 1.0;
    }else{
        *model = -1.0;
    }
    tablo[0][0] = 0;// inicia z = 0
    
    
    printf("Valores da FO:\n");
    for (int c = 1; c <= *n_variaveis; c++) {// lendo FO
        printf("X%d: ", c);
        cin >> tablo[0][c];
        tablo[0][c] *= *model;// zj-cj
    }
    
    for (int l = 1; l <= *n_restricoes; l++) {// lendo restricoes
        printf("\nRestrição %d:\n", l);
        for (int c = 1; c <= *n_variaveis; c++) {
            printf("X%d: ", c);
            cin >> tablo[l][c];
        }
        printf("Restrição %d <= ", l);
        cin >> tablo[l][0];
    }
    return tablo;
}

void imprimirModelo(double model, int n_variaveis, int n_restricoes, vector<vector<double> > tablo){
    cout << "\nModelo\n";
    
    if(model == 1.0){
        cout << "Max Z = ";
    }else{
        cout << "Min Z = ";
    }
    
    for(int c = 1; c <= n_variaveis; c++){
        printf("%0.fX%d", tablo[0][c] * model , c);
        if(c < n_variaveis){
            cout << " + ";
        }
    }
    cout << endl;
    
    for(int l = 1; l <= n_restricoes; l++){
        for(int c = 1; c <= n_variaveis; c++){
            printf("%0.fX%d", tablo[l][c] , c);
            if(c < n_variaveis){
                cout << " + ";
            }
        }
        printf(" <= %0.f\n", tablo[l][0] );
    }
    cout << endl;
}