int main() {
init();
// read input
cin >> n >> m;
for (int i=0;i<m;i++) {
for (int j=0;j<3;j++) cin >> e[i][j];
nextE[e[i][0]].push_back(i);
nextE[e[i][1]].push_back(i);
next[e[i][0]].push_back(e[i][1]);
next[e[i][1]].push_back(e[i][0]);
}
// set up linear program
nA=0;
for (int i=0;i<n;i++) cycles(i);
for (int i=0;i<m;i++) {
A[nA][2*i+1]=1; b[nA]=e[i][2];
nA++;
c[2*i]=-1; c[2*i+1]=1;
}
// solve
Simplex s = Simplex(A, b, c, nA, 2*m);
Fraction f = s.solve();
f = f.neg();
for (int i=0;i<m;i++) f = f.add(Fraction(e[i][2]));
cout << f.str() << endl;
return 0;
}
int main() {
/*
Simplex::pivot_struct test_struct;
Simplex::pivot_struct out_struct;
test_struct.v=0;
test_struct.N.insert(1);
test_struct.N.insert(2);
test_struct.N.insert(3);
test_struct.B.insert(4);
test_struct.B.insert(5);
test_struct.B.insert(6);
test_struct.c[1]=3;
test_struct.c[2]=1;
test_struct.c[3]=2;
test_struct.b[4]=30;
test_struct.b[5]=24;
test_struct.b[6]=36;
std::map<int,float> A_l1;
std::map<int,float> A_l2;
std::map<int,float> A_l3;
test_struct.A[4]=A_l1;
test_struct.A[5]=A_l2;
test_struct.A[6]=A_l3;
test_struct.A[4][1]=1;
test_struct.A[5][1]=2;
test_struct.A[6][1]=4;
test_struct.A[4][2]=1;
test_struct.A[5][2]=2;
test_struct.A[6][2]=1;
test_struct.A[4][3]=3;
test_struct.A[5][3]=5;
test_struct.A[6][3]=2;
int e = 1;
int l = 6;
//std::cout << "Before Pivot" <<std::endl;
//Simplex::print_pivot_struct(test_struct);
out_struct = Simplex::pivot(test_struct,e,l);
//std::cout << "After Pivot" <<std::endl;
//Simplex::print_pivot_struct(out_struct);
*/
/*
std::vector<float> b;
b.push_back(30);
b.push_back(24);
b.push_back(36);
std::vector<float> c;
c.push_back(3);
c.push_back(1);
c.push_back(2);
std::vector<std::vector<float>> A;
A.push_back(std::vector<float> ());
A.push_back(std::vector<float> ());
A.push_back(std::vector<float> ());
A[0].push_back(1);
A[0].push_back(1);
A[0].push_back(3);
A[1].push_back(2);
A[1].push_back(2);
A[1].push_back(5);
A[2].push_back(4);
A[2].push_back(1);
A[2].push_back(2);
std::cout << "Create Simplex" <<std::endl;
Simplex teste = Simplex(A,b,c);
std::cout << "Print" <<std::endl;
teste.print_internal_struct();
std::vector<float> x = teste.solve();
std::cout << "Solution" <<std::endl;
for(auto i : x){
std::cout << "x = " << i << std::endl;
}
std::cout << "Without initial feasible solution: " <<std::endl;
std::vector<float> b1;
b1.push_back(2);
b1.push_back(-4);
std::vector<float> c1;
c1.push_back(2);
c1.push_back(-1);
std::vector<std::vector<float>> A1;
A1.push_back(std::vector<float> ());
A1.push_back(std::vector<float> ());
A1[0].push_back(2);
A1[0].push_back(-1);
A1[1].push_back(1);
A1[1].push_back(-5);
Simplex teste1 = Simplex(A1,b1,c1);
std::cout << "Print" <<std::endl;
teste1.print_internal_struct();
std::vector<float> x1 = teste1.solve();
std::cout << "Solution" <<std::endl;
std::cout << "Size: " << x1.size() << std::endl;
for(auto i : x1){
std::cout << "x = " << i << std::endl;
}
std::cout << std::endl<< std::endl;
*/
// TEST with Arduino LP
/*
# Created by Octave 3.8.1, Thu Dec 10 09:46:19 2015 WET <jmirandadealme@JMdeA-Zenbook>
# name: A
# type: matrix
# rows: 18
# columns: 9
-10.21830350832508 -0.7271861952428118 -0.09748519653650334 -0.7271861952428118 -0.2769996701683155 -0.06890464691960219 -0.09748519653650334 -0.06890464691960219 -0.03280392761204531
-0.7271861952428118 -10.21830350832508 -0.7271861952428118 -0.2769996701683155 -0.7271861952428118 -0.2769996701683155 -0.06890464691960219 -0.09748519653650334 -0.06890464691960219
-0.09748519653650334 -0.7271861952428118 -10.21830350832508 -0.06890464691960219 -0.2769996701683155 -0.7271861952428118 -0.03280392761204531 -0.06890464691960219 -0.09748519653650334
-0.7271861952428118 -0.2769996701683155 -0.06890464691960219 -10.21830350832508 -0.7271861952428118 -0.09748519653650334 -0.7271861952428118 -0.2769996701683155 -0.06890464691960219
-0.2769996701683155 -0.7271861952428118 -0.2769996701683155 -0.7271861952428118 -10.21830350832508 -0.7271861952428118 -0.2769996701683155 -0.7271861952428118 -0.2769996701683155
-0.06890464691960219 -0.2769996701683155 -0.7271861952428118 -0.09748519653650334 -0.7271861952428118 -10.21830350832508 -0.06890464691960219 -0.2769996701683155 -0.7271861952428118
-0.09748519653650334 -0.06890464691960219 -0.03280392761204531 -0.7271861952428118 -0.2769996701683155 -0.06890464691960219 -10.21830350832508 -0.7271861952428118 -0.09748519653650334
-0.06890464691960219 -0.09748519653650334 -0.06890464691960219 -0.2769996701683155 -0.7271861952428118 -0.2769996701683155 -0.7271861952428118 -10.21830350832508 -0.7271861952428118
-0.03280392761204531 -0.06890464691960219 -0.09748519653650334 -0.06890464691960219 -0.2769996701683155 -0.7271861952428118 -0.09748519653650334 -0.7271861952428118 -10.21830350832508
1 0 0 0 0 0 0 0 0
0 1 0 0 0 0 0 0 0
0 0 1 0 0 0 0 0 0
0 0 0 1 0 0 0 0 0
0 0 0 0 1 0 0 0 0
0 0 0 0 0 1 0 0 0
0 0 0 0 0 0 1 0 0
0 0 0 0 0 0 0 1 0
0 0 0 0 0 0 0 0 1
*/
/*
# name: b
# type: matrix
# rows: 18
# columns: 1
-2.234087037518483
-6.370072142271482
-2.314137495496334
-4.195003918678363
-2.89337934873546
-7.706424278239382
-7.357001197283779
-5.483351507948527
-7.286096857845781
1
1
1
1
1
1
1
1
1
*/
/*
c = -1 -1 -1 -1 -1 -1 -1 -1 -1
*/
/*
# name: random_pwm
# type: matrix
# rows: 9
# columns: 1
0.1374661759009127
0.5587605221930771
0.1228113038292256
0.3096195642559854
0.104062656935522
0.659452543924277
0.6504349725253338
0.4053939889252992
0.6207688892129481
Solution
x = 0.137466
x = 0.55876
x = 0.122811
x = 0.309627
x = 0.104059
x = 0.659453
x = 0.650444
x = 0.405394
x = 0.620769
*/
/*
std::vector<std::vector<float>> A2;
A2.push_back(std::vector<float> ());
A2.push_back(std::vector<float> ());
A2.push_back(std::vector<float> ());
A2.push_back(std::vector<float> ());
A2.push_back(std::vector<float> ());
A2.push_back(std::vector<float> ());
A2.push_back(std::vector<float> ());
A2.push_back(std::vector<float> ());
A2.push_back(std::vector<float> ());
A2.push_back(std::vector<float> ());
A2.push_back(std::vector<float> ());
A2.push_back(std::vector<float> ());
A2.push_back(std::vector<float> ());
A2.push_back(std::vector<float> ());
A2.push_back(std::vector<float> ());
A2.push_back(std::vector<float> ());
A2.push_back(std::vector<float> ());
A2.push_back(std::vector<float> ());
A2[0].push_back(-10.2183035083);
A2[0].push_back(-0.7271861952);
A2[0].push_back(-0.0974851965);
A2[0].push_back(-0.7271861952);
A2[0].push_back(-0.2769996702);
A2[0].push_back(-0.0689046469);
A2[0].push_back(-0.0974851965);
A2[0].push_back(-0.0689046469);
A2[0].push_back(-0.0328039276);
A2[1].push_back(-0.7271861952428118);
A2[1].push_back(-10.21830350832508);
A2[1].push_back(-0.7271861952428118);
A2[1].push_back(-0.2769996701683155);
A2[1].push_back(-0.7271861952428118);
A2[1].push_back(-0.2769996701683155);
A2[1].push_back(-0.06890464691960219);
A2[1].push_back(-0.09748519653650334);
A2[1].push_back(-0.06890464691960219);
A2[2].push_back(-0.09748519653650334);
A2[2].push_back(-0.7271861952428118);
A2[2].push_back(-10.21830350832508);
A2[2].push_back(-0.06890464691960219);
A2[2].push_back(-0.2769996701683155);
A2[2].push_back(-0.7271861952428118);
A2[2].push_back(-0.03280392761204531);
A2[2].push_back(-0.06890464691960219);
A2[2].push_back(-0.09748519653650334);
A2[3].push_back(-0.7271861952428118);
A2[3].push_back(-0.2769996701683155 );
A2[3].push_back(-0.06890464691960219);
A2[3].push_back(-10.21830350832508);
A2[3].push_back(-0.7271861952428118);
A2[3].push_back(-0.09748519653650334);
A2[3].push_back(-0.7271861952428118);
A2[3].push_back(-0.2769996701683155);
A2[3].push_back(-0.06890464691960219);
A2[4].push_back(-0.2769996701683155);
A2[4].push_back(-0.7271861952428118);
A2[4].push_back(-0.2769996701683155);
A2[4].push_back(-0.7271861952428118);
A2[4].push_back(-10.21830350832508);
A2[4].push_back(-0.7271861952428118);
A2[4].push_back(-0.2769996701683155);
A2[4].push_back(-0.7271861952428118);
A2[4].push_back(-0.2769996701683155);
A2[5].push_back(-0.06890464691960219);
A2[5].push_back(-0.2769996701683155);
A2[5].push_back(-0.7271861952428118);
A2[5].push_back(-0.09748519653650334);
A2[5].push_back(-0.7271861952428118);
A2[5].push_back(-10.21830350832508);
A2[5].push_back(-0.06890464691960219);
A2[5].push_back(-0.2769996701683155);
A2[5].push_back(-0.7271861952428118);
A2[6].push_back(-0.09748519653650334);
A2[6].push_back(-0.06890464691960219);
A2[6].push_back(-0.03280392761204531);
A2[6].push_back(-0.7271861952428118);
A2[6].push_back(-0.2769996701683155);
A2[6].push_back(-0.06890464691960219);
A2[6].push_back(-10.21830350832508);
A2[6].push_back(-0.7271861952428118);
A2[6].push_back(-0.09748519653650334);
A2[7].push_back(-0.06890464691960219);
A2[7].push_back(-0.09748519653650334);
A2[7].push_back(-0.06890464691960219);
A2[7].push_back(-0.2769996701683155);
A2[7].push_back(-0.7271861952428118);
A2[7].push_back(-0.2769996701683155);
A2[7].push_back(-0.7271861952428118);
A2[7].push_back(-10.21830350832508);
A2[7].push_back(-0.7271861952428118);
A2[8].push_back(-0.03280392761204531);
A2[8].push_back(-0.06890464691960219);
A2[8].push_back(-0.09748519653650334);
A2[8].push_back(-0.06890464691960219);
A2[8].push_back(-0.2769996701683155);
A2[8].push_back(-0.7271861952428118);
A2[8].push_back(-0.09748519653650334);
A2[8].push_back(-0.7271861952428118);
A2[8].push_back(-10.21830350832508);
A2[9].push_back(1);
A2[9].push_back(0);
A2[9].push_back(0);
A2[9].push_back(0);
A2[9].push_back(0);
A2[9].push_back(0);
A2[9].push_back(0);
A2[9].push_back(0);
A2[9].push_back(0);
A2[10].push_back(0);
A2[10].push_back(1);
A2[10].push_back(0);
A2[10].push_back(0);
A2[10].push_back(0);
A2[10].push_back(0);
A2[10].push_back(0);
A2[10].push_back(0);
A2[10].push_back(0);
A2[11].push_back(0);
A2[11].push_back(0);
A2[11].push_back(1);
A2[11].push_back(0);
A2[11].push_back(0);
A2[11].push_back(0);
A2[11].push_back(0);
A2[11].push_back(0);
A2[11].push_back(0);
A2[12].push_back(0);
A2[12].push_back(0);
A2[12].push_back(0);
A2[12].push_back(1);
A2[12].push_back(0);
A2[12].push_back(0);
A2[12].push_back(0);
A2[12].push_back(0);
A2[12].push_back(0);
A2[13].push_back(0);
A2[13].push_back(0);
A2[13].push_back(0);
A2[13].push_back(0);
A2[13].push_back(1);
A2[13].push_back(0);
A2[13].push_back(0);
A2[13].push_back(0);
A2[13].push_back(0);
A2[14].push_back(0);
A2[14].push_back(0);
A2[14].push_back(0);
A2[14].push_back(0);
A2[14].push_back(0);
A2[14].push_back(1);
A2[14].push_back(0);
A2[14].push_back(0);
A2[14].push_back(0);
A2[15].push_back(0);
A2[15].push_back(0);
A2[15].push_back(0);
A2[15].push_back(0);
A2[15].push_back(0);
A2[15].push_back(0);
A2[15].push_back(1);
A2[15].push_back(0);
A2[15].push_back(0);
A2[16].push_back(0);
A2[16].push_back(0);
A2[16].push_back(0);
A2[16].push_back(0);
A2[16].push_back(0);
A2[16].push_back(0);
A2[16].push_back(0);
A2[16].push_back(1);
A2[16].push_back(0);
A2[17].push_back(0);
A2[17].push_back(0);
A2[17].push_back(0);
A2[17].push_back(0);
A2[17].push_back(0);
A2[17].push_back(0);
A2[17].push_back(0);
A2[17].push_back(0);
A2[17].push_back(1);
std::vector<float> b2;
b2.push_back(-2.234087037518483);
b2.push_back(-6.370072142271482);
b2.push_back(-2.314137495496334);
b2.push_back(-4.195003918678363);
b2.push_back(-2.89337934873546);
b2.push_back(-7.706424278239382);
b2.push_back(-7.357001197283779);
b2.push_back(-5.483351507948527);
b2.push_back(-7.286096857845781);
b2.push_back(1);
b2.push_back(1);
b2.push_back(1);
b2.push_back(1);
b2.push_back(1);
b2.push_back(1);
b2.push_back(1);
b2.push_back(1);
b2.push_back(1);
std::vector<float> c2;
c2.push_back(-1);
c2.push_back(-1);
c2.push_back(-1);
c2.push_back(-1);
c2.push_back(-1);
c2.push_back(-1);
c2.push_back(-1);
c2.push_back(-1);
c2.push_back(-1);
Simplex teste2 = Simplex(A2,b2,c2);
std::cout << "Print" <<std::endl;
teste2.print_internal_struct();
std::vector<float> x2 = teste2.solve();
std::cout << "Solution" <<std::endl;
for(auto i : x2){
std::cout << "x = " << i << std::endl;
}
*/
/* // Mr Silva test, lab 12:
std::vector<std::vector<float>> A3;
A3.push_back(std::vector<float> ());
A3.push_back(std::vector<float> ());
A3.push_back(std::vector<float> ());
A3[0].push_back(2);
A3[0].push_back(1);
A3[1].push_back(4);
A3[1].push_back(3);
std::vector<float> b3;
b3.push_back(100);
b3.push_back(240);
std::vector<float> c3;
c3.push_back(7);
c3.push_back(5);
Simplex teste3 = Simplex(A3,b3,c3);
std::cout << "Print" <<std::endl;
teste3.print_internal_struct();
std::vector<float> x3 = teste3.solve();
std::cout << "Solution" <<std::endl;
for(auto i : x3){
std::cout << "x = " << i << std::endl;
}
*/
// Mr Silva with unnecessary constrain
/*
std::vector<std::vector<float>> A4;
A4.push_back(std::vector<float> ());
A4.push_back(std::vector<float> ());
A4.push_back(std::vector<float> ());
A4[0].push_back(-1);
A4[0].push_back(0);
A4[1].push_back(2);
A4[1].push_back(1);
A4[2].push_back(4);
A4[2].push_back(3);
std::vector<float> b4;
b4.push_back(-20);
b4.push_back(100);
b4.push_back(240);
std::vector<float> c4;
c4.push_back(7);
c4.push_back(5);
Simplex teste4 = Simplex(A4,b4,c4);
std::cout << "Print" <<std::endl;
teste4.print_internal_struct();
std::vector<float> x4 = teste4.solve();
std::cout << "Solution" <<std::endl;
for(auto i : x4){
std::cout << "x = " << i << std::endl;
}
*/
/*
N = {2,5,}
B = {1,3,4,}
A = {
0,-1,
1,2,
3,4,
}
b = {20,60,160,}
c = {5,7,}
v = 140
*/
/*
std::vector<std::vector<float>> A5;
A5.push_back(std::vector<float> ());
A5.push_back(std::vector<float> ());
A5.push_back(std::vector<float> ());
A5[0].push_back(0);
A5[0].push_back(-1);
A5[1].push_back(1);
A5[1].push_back(2);
A5[2].push_back(3);
A5[2].push_back(4);
std::vector<float> b5;
b5.push_back(20);
b5.push_back(60);
b5.push_back(160);
std::vector<float> c5;
c5.push_back(5);
c5.push_back(7);
Simplex teste5 = Simplex(A5,b5,c5);
std::vector<float> x5 = teste5.solve();
std::cout << "Solution" <<std::endl;
for(auto i : x5){
std::cout << "x = " << i << std::endl;
}
*/
/*isto é para aberto
O_vals: 198 236 260
E_vals:
708 275 269
579 640 284
380 560 581
*/
/*
std::vector<std::vector<float>> E;
std::vector<float> O;
std::vector<float> L;
E.push_back(std::vector<float>());
E.push_back(std::vector<float>());
E.push_back(std::vector<float>());
O.push_back(198.0/1024);
O.push_back(236.0/1024);
O.push_back(260.0/1024);
L.push_back(600.0/1024);
L.push_back(600.0/1024);
L.push_back(600.0/1024);
E[0].push_back(708.0/1024);
E[0].push_back(275.0/1024);
E[0].push_back(269.0/1024);
E[1].push_back(579.0/1024);
E[1].push_back(640.0/1024);
E[1].push_back(284.0/1024);
E[2].push_back(380.0/1024);
E[2].push_back(560.0/1024);
E[2].push_back(581.0/1024);
lp_struct lp = generate_linear_programm(E,O,L);
Simplex teste = Simplex(lp.A,lp.b,lp.c);
std::vector<float> x = teste.solve();
for(auto i : x){
std::cout << "x = " << i*255 << std::endl;
}
*/
/*fechado
O_vals: 5 3 3
E_vals:
745 137 67
598 704 113
317 586 674*/
std::vector<std::vector<float>> E;
std::vector<float> O;
std::vector<float> L;
E.push_back(std::vector<float>());
E.push_back(std::vector<float>());
E.push_back(std::vector<float>());
O.push_back(5.0/1024);
O.push_back(3.0/1024);
O.push_back(3.0/1024);
L.push_back(600.0/1024);
L.push_back(600.0/1024);
L.push_back(600.0/1024);
E[0].push_back(745.0/1024);
E[0].push_back(137.0/1024);
E[0].push_back(67.0/1024);
E[1].push_back(598.0/1024);
E[1].push_back(704.0/1024);
E[1].push_back(113.0/1024);
E[2].push_back(317.0/1024);
E[2].push_back(586.0/1024);
E[2].push_back(674.0/1024);
lp_struct lp = generate_linear_programm(E,O,L);
Simplex teste = Simplex(lp.A,lp.b,lp.c);
std::vector<float> x = teste.solve();
for(auto i : x){
std::cout << "x = " << i*255 << std::endl;
}
}
bool ReductionToLP::solve()
{
Simplex simplex (var,res, C, B, A,min, EQ);
simplex.solve();
return 0;
}
