#include <glpk.h>
#include <math.h>
#include <stdio.h>

#define int64 long long int

int main()
{
    int j;

    glp_prob *lp = glp_create_prob();

    glp_read_lp(lp, NULL, "test.cplex");

    glp_smcp parm;
    glp_init_smcp(&parm);
    parm.presolve = GLP_ON;

    glp_scale_prob(lp, GLP_SF_AUTO);
    glp_adv_basis(lp, 0);
    glp_simplex(lp, &parm);

    int lp_status = glp_get_status(lp);

    if (lp_status == GLP_NOFEAS) {
        glp_delete_prob(lp);
        return 1;
    }

    glp_intopt(lp, NULL);

    int ilp_status = glp_mip_status(lp);

    if (ilp_status == GLP_NOFEAS) {
        glp_delete_prob(lp);
        return 1;
    }

    double z = glp_mip_obj_val(lp);
    printf("z = %lf\n", z);

    for (j=0; j<glp_get_num_cols(lp); j++) {
        double x = glp_mip_col_val(lp, j+1);
        printf("c%d = %lld, ", j, (int64) round(x));
    }
    printf("\n");;
    glp_delete_prob(lp);

    return 0;
}