|
1 //#include <stdio.h> |
|
2 //#include <stdlib.h> |
|
3 #include <iostream> |
|
4 #include <lemon/lp_glpk.h> |
|
5 using namespace lemon; |
|
6 |
|
7 int main() |
|
8 { |
|
9 LpGlpk lp; |
|
10 typedef LpGlpk::Row Row; |
|
11 typedef LpGlpk::Col Col; |
|
12 //std::cout << "Hello" << std::endl; |
|
13 |
|
14 // LPX *lp; |
|
15 // int ia[1+1000], ja[1+1000]; |
|
16 // double ar[1+1000], Z, x1, x2, x3; |
|
17 // s1: lp = lpx_create_prob(); |
|
18 // s2: lpx_set_prob_name(lp, "sample"); |
|
19 // s3: lpx_set_obj_dir(lp, LPX_MAX); |
|
20 lp.max(); |
|
21 // s4: lpx_add_rows(lp, 3); |
|
22 // s5: lpx_set_row_name(lp, 1, "p"); |
|
23 // s6: lpx_set_row_bnds(lp, 1, LPX_UP, 0.0, 100.0); |
|
24 |
|
25 // Row p = lp.addRow(); |
|
26 |
|
27 // s7: lpx_set_row_name(lp, 2, "q"); |
|
28 // s8: lpx_set_row_bnds(lp, 2, LPX_UP, 0.0, 600.0); |
|
29 |
|
30 // Row q = lp.addRow(); |
|
31 |
|
32 // s9: lpx_set_row_name(lp, 3, "r"); |
|
33 // s10: lpx_set_row_bnds(lp, 3, LPX_UP, 0.0, 300.0); |
|
34 |
|
35 // Row r = lp.addRow(); |
|
36 |
|
37 // s11: lpx_add_cols(lp, 3); |
|
38 // s12: lpx_set_col_name(lp, 1, "x1"); |
|
39 // s13: lpx_set_col_bnds(lp, 1, LPX_LO, 0.0, 0.0); |
|
40 |
|
41 // Col x1 = lp.addCol(); |
|
42 |
|
43 // s14: lpx_set_obj_coef(lp, 1, 10.0); |
|
44 // s15: lpx_set_col_name(lp, 2, "x2"); |
|
45 // s16: lpx_set_col_bnds(lp, 2, LPX_LO, 0.0, 0.0); |
|
46 // s17: lpx_set_obj_coef(lp, 2, 6.0); |
|
47 // s18: lpx_set_col_name(lp, 3, "x3"); |
|
48 // s19: lpx_set_col_bnds(lp, 3, LPX_LO, 0.0, 0.0); |
|
49 // s20: lpx_set_obj_coef(lp, 3, 4.0); |
|
50 // s21: ia[1] = 1, ja[1] = 1, ar[1] = 1.0; /* a[1,1] = 1 */ |
|
51 // s22: ia[2] = 1, ja[2] = 2, ar[2] = 1.0; /* a[1,2] = 1 */ |
|
52 // s23: ia[3] = 1, ja[3] = 3, ar[3] = 1.0; /* a[1,3] = 1 */ |
|
53 // s24: ia[4] = 2, ja[4] = 1, ar[4] = 10.0; /* a[2,1] = 10 */ |
|
54 // s25: ia[5] = 3, ja[5] = 1, ar[5] = 2.0; /* a[3,1] = 2 */ |
|
55 // s26: ia[6] = 2, ja[6] = 2, ar[6] = 4.0; /* a[2,2] = 4 */ |
|
56 // s27: ia[7] = 3, ja[7] = 2, ar[7] = 2.0; /* a[3,2] = 2 */ |
|
57 // s28: ia[8] = 2, ja[8] = 3, ar[8] = 5.0; /* a[2,3] = 5 */ |
|
58 // s29: ia[9] = 3, ja[9] = 3, ar[9] = 6.0; /* a[3,3] = 6 */ |
|
59 // s30: lpx_load_matrix(lp, 9, ia, ja, ar); |
|
60 // s31: lpx_simplex(lp); |
|
61 // s32: Z = lpx_get_obj_val(lp); |
|
62 // s33: x1 = lpx_get_col_prim(lp, 1); |
|
63 // s34: x2 = lpx_get_col_prim(lp, 2); |
|
64 // s35: x3 = lpx_get_col_prim(lp, 3); |
|
65 // s36: printf("\nZ = %g; x1 = %g; x2 = %g; x3 = %g\n", Z, x1, x2, x3); |
|
66 // s37: lpx_delete_prob(lp); |
|
67 // return 0; |
|
68 return 0; |
|
69 } |