Dev - C++ - por favor me pueden ayudar con este ejercicio

 
Vista:

por favor me pueden ayudar con este ejercicio

Publicado por isak (10 intervenciones) el 18/07/2020 02:54:25
realizar un código de programación lineal que resuelva ejercicios con el método simplex dual
aqui ya tenia algo pero tiene algunos errores

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
137
138
139
140
141
142
143
144
145
146
147
148
#include <iostream>
#include <stdlib.h>
#include <stdio.h>
using namespace std;
double funcion[100], restricciones[100][100];
double restricc[100][100];
double independiente[100];
 
int main(){
 
    int rest, var, ni ;
    int i,j, k, columna, renglon=0, count, cuenta;
    double maximo = 0, menor, aux,  temp;
    int c;
	cout<<" Funcion a Optimizar \nNumero de Variables: ";
	cin>>var;
	cout<<"\nFUNCION OBJETIVO"<<endl;
	for(c=0;c<var;c++){
		cin>>funcion[c];
		}
		cout<<"\nNumero de Restricciones: ";
		cin>>rest;
 
		for(i=0;i<rest;i++){
			cout<<"Restriccion "<<i+1<<endl;
		for(c=0;c<var;c++){
			cin>>restricc[i][c];
			}
 
    cout<<"Termino independiente ";
    cin>>independiente[i];
   }
   	   for( i=0;i<var;i++)
		   for( j=0;j<rest;j++)
		   if(restricc[j][i]!=0)
			   restricciones[j][i]=restricc[j][i]*-1;
		   else
               restricciones[j][i]=restricc[j][i];
    //Crear tabla inicial simplex
    for(i=0;i<rest;i++){
        for(j=0; j<rest; j++)
           if (i==j)
           restricciones[i][var+j]=1;
           else
           restricciones[i][var+j]=0;
   }
 
       for(i=0; i<rest;i++)
           restricciones[i][var+rest] = independiente[i]*-1;
 
        for(j=0; j<var; j++)
			restricciones[rest][j] = funcion[j]*-1;
        for(j=var; j<rest+var; j++)
           restricciones[rest][j] = 0;
 
    ni=1;
    cout<<"\n-----------------------------------------------";
    cout<<"\n ITERACION "<<ni;
    cout<<"\n-----------------------------------------------\n";
    ni++;
    //rest cantidad de restricciones para la tabla
    for(i=0; i<=rest; i++) {
		for(j=0; j<=rest+var; j++)
		cout<<"  \t"<<restricciones[i][j];
		cout<<endl;
		}
		//Encontrar la variable de decision que entra a la base
		//y la de holgura que saldra de la base
		do{
			///Encontrar la variable de decision que entrara a la base
            maximo = abs(restricciones[rest][0]);
			columna = 0;
            count = 0;
		    aux = 120 ;
			menor =120;
			for(i=0; i<rest; i++) {
				if(restricciones[i][columna] < 0)
				aux = restricciones[i][rest+var] / restricciones[i][columna];
				else
				count++;
				if( menor < aux ) {
					menor = aux;
					renglon = i;
				}
				else
			continue;
			}
			for(j=0; j<=rest+var; j++)
				if( restricciones[rest][j] < 0 ){
				temp = -1*restricciones[rest][j];
				if ( maximo < temp){
					maximo = temp;
					columna = j;
				}
		     else
			continue;
			}
			for(i=0; i<rest; i++) {
				if(restricciones[i][columna] < 0)
				aux = restricciones[i][rest+var] / restricciones[i][columna];
				else
				count++;
				if( menor > aux ) {
					menor = aux;
					renglon = i;
				}
			}
			if(count==var) {
			cout<<"Solucion no acotada";
			exit(1);
			}
			cout<<"\nPivote:"<<restricciones[renglon][columna];
     //Encontrar los coeficientes de la nueva tabla
     aux = restricciones[renglon][columna];
     for(j=0; j<=(rest+var); j++)
     restricciones[renglon][j] = restricciones[renglon][j] / aux;
     for(i=0; i<=rest; i++) {
		 if (i == renglon ){
		 continue;}
		 else {
			 aux = restricciones[i][columna]*-1;
 
			 for(k=0; k <= rest+var; k++) {
 
				 restricciones[i][k] = (aux *restricciones[renglon][k]) + restricciones[i][k];
				 }
		      }
	}
	cout<<"\n-----------------------------------------------";
	cout<<"\n ITERACION "<<ni;
	cout<<"\n-----------------------------------------------\n";
	ni++;
	///la nueva tabla
	for(i=0; i<=rest; i++) {
		for(j=0; j<=rest+var; j++)
		cout<<"   \t"<<restricciones[i][j];
		cout<<"\n";
	}
	cuenta=0;
	for(j=0;j<rest+var;j++)
	if(restricciones[var][j] >= 0)
	cuenta++;
	if (cuenta == rest+var)
	break;
	}while(cuenta);
 
	return 0;
}
Valora esta pregunta
Me gusta: Está pregunta es útil y esta claraNo me gusta: Está pregunta no esta clara o no es útil
0
Responder
Imágen de perfil de Rodrigo
Val: 1.755
Plata
Ha mantenido su posición en Dev - C++ (en relación al último mes)
Gráfica de Dev - C++

por favor me pueden ayudar con este ejercicio

Publicado por Rodrigo (539 intervenciones) el 18/07/2020 03:30:48
SI quieres ayuda no digas solamente que "tiene algunos errores". Explicitalos.
Valora esta respuesta
Me gusta: Está respuesta es útil y esta claraNo me gusta: Está respuesta no esta clara o no es útil
0
Comentar

por favor me pueden ayudar con este ejercicio

Publicado por Isak (10 intervenciones) el 18/07/2020 05:19:25
El problema esta que cuando compilo y introdusco los datos solo realiza 2 iteraciones y va copiando lo mismo sin detenerse
Tiene que realizar las iteraciones hasta que los independientes que son multiplicado por -1 sean mayores o iguales a cero
Por ejemplo con
2 variables
2000 2000
3 restricciones
Restriccion1: 1 2 independiente 80

Restriccion2: 3 2 independiente 160

Restriccion3: 5 2 independiente 200


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
137
138
139
140
141
142
143
144
145
146
147
148
#include <iostream>
#include <stdlib.h>
#include <stdio.h>
using namespace std;
double funcion[100], restricciones[100][100];
double restricc[100][100];
double independiente[100];
 
int main(){
 
int rest, var, ni ;
int i,j, k, columna, renglon=0, count, cuenta;
double maximo = 0, menor, aux, temp;
int c;
cout<<" Funcion a Optimizar \nNumero de Variables: ";
cin>>var;
cout<<"\nFUNCION OBJETIVO"<<endl;
for(c=0;c<var;c++){
cin>>funcion[c];
}
cout<<"\nNumero de Restricciones: ";
cin>>rest;
 
for(i=0;i<rest;i++){
cout<<"Restriccion "<<i+1<<endl;
for(c=0;c<var;c++){
cin>>restricc[i][c];
}
 
cout<<"Termino independiente ";
cin>>independiente[i];
}
for( i=0;i<var;i++)
for( j=0;j<rest;j++)
if(restricc[j][i]!=0)
restricciones[j][i]=restricc[j][i]*-1;
else
restricciones[j][i]=restricc[j][i];
//Crear tabla inicial simplex
for(i=0;i<rest;i++){
for(j=0; j<rest; j++)
if (i==j)
restricciones[i][var+j]=1;
else
restricciones[i][var+j]=0;
}
 
for(i=0; i<rest;i++)
restricciones[i][var+rest] = independiente[i]*-1;
 
for(j=0; j<var; j++)
restricciones[rest][j] = funcion[j]*-1;
for(j=var; j<rest+var; j++)
restricciones[rest][j] = 0;
 
ni=1;
cout<<"\n-----------------------------------------------";
cout<<"\n ITERACION "<<ni;
cout<<"\n-----------------------------------------------\n";
ni++;
//rest cantidad de restricciones para la tabla
for(i=0; i<=rest; i++) {
for(j=0; j<=rest+var; j++)
cout<<" \t"<<restricciones[i][j];
cout<<endl;
}
//Encontrar la variable de decision que entra a la base
//y la de holgura que saldra de la base
do{
///Encontrar la variable de decision que entrara a la base
maximo = abs(restricciones[rest][0]);
columna = 0;
count = 0;
aux = 120 ;
menor =120;
for(i=0; i<rest; i++) {
if(restricciones[i][columna] < 0)
aux = restricciones[i][rest+var] / restricciones[i][columna];
else
count++;
if( menor < aux ) {
menor = aux;
renglon = i;
}
else
continue;
}
for(j=0; j<=rest+var; j++)
if( restricciones[rest][j] < 0 ){
temp = -1*restricciones[rest][j];
if ( maximo < temp){
maximo = temp;
columna = j;
}
else
continue;
}
for(i=0; i<rest; i++) {
if(restricciones[i][columna] < 0)
aux = restricciones[i][rest+var] / restricciones[i][columna];
else
count++;
if( menor > aux ) {
menor = aux;
renglon = i;
}
}
if(count==var) {
cout<<"Solucion no acotada";
exit(1);
}
cout<<"\nPivote:"<<restricciones[renglon][columna];
//Encontrar los coeficientes de la nueva tabla
aux = restricciones[renglon][columna];
for(j=0; j<=(rest+var); j++)
restricciones[renglon][j] = restricciones[renglon][j] / aux;
for(i=0; i<=rest; i++) {
if (i == renglon ){
continue;}
else {
aux = restricciones[i][columna]*-1;
 
for(k=0; k <= rest+var; k++) {
 
restricciones[i][k] = (aux *restricciones[renglon][k]) + restricciones[i][k];
}
}
}
cout<<"\n-----------------------------------------------";
cout<<"\n ITERACION "<<ni;
cout<<"\n-----------------------------------------------\n";
ni++;
///la nueva tabla
for(i=0; i<=rest; i++) {
for(j=0; j<=rest+var; j++)
cout<<" \t"<<restricciones[i][j];
cout<<"\n";
}
cuenta=0;
for(j=0;j<rest+var;j++)
if(restricciones[var][j] >= 0)
cuenta++;
if (cuenta == rest+var)
break;
}while(cuenta);
 
return 0;
}
Valora esta respuesta
Me gusta: Está respuesta es útil y esta claraNo me gusta: Está respuesta no esta clara o no es útil
0
Comentar