void error(void)
{
	int ia[] = {0, 1, 2, 2, 2, 3, 3, 4, 4, 5, 5, 6, 6, 7, 7, 8, 8, 9, 9, 10, 10, 10, 11, 11, 12, 12, 12, 13, 13, 14,
		14, 14, 15, 15, 16, 16, 16, 17, 17, 18, 18, 19, 19, 20, 20, 21, 21, 22, 22, 23, 23, 24, 24, 25, 25, 26, 26, 27, 27,
		28, 28, 29, 29, 30, 30, 31, 31, 32, 32, 33, 33, 34, 34, 35, 35, 36, 36, 37, 37, 38, 38, 39, 39, 40, 40, 40, 41, 41,
		42, 42, 43, 43, 44, 44, 45 };
	int ja[] = {0, 1, 1, 2, 7, 1, 2, 2, 3, 2, 3, 3, 4, 3, 4, 4, 5, 4, 5, 5, 6, 7, 5, 6, 6, 7, 21, 6, 7, 7, 8, 9, 7,
		8, 8, 9, 24, 8, 9, 9, 10, 9, 10, 10, 11, 10, 11, 11, 12, 11, 12, 12, 13, 12, 13, 13, 14, 13, 14, 14, 15, 14, 15, 15,
		16, 15, 16, 16, 17, 16, 17, 17, 18, 17, 18, 18, 19, 18, 19, 19, 20, 19, 20, 20, 22, 23, 20, 21, 21, 22, 22, 23, 23,
		24, 24 };
	double ar[] = {0.0, 1.0, 1.0, -1.0, -1.0, 1.0, -1.0, 1.0, -2.0, 1.0, -1.0, 1.0, -2.0, 1.0, -1.0, 1.0, -2.0, 1.0,
		-1.0, 1.0, -1.0, -1.0, 1.0, -1.0, 1.0, -1.0, -1.0, 1.0, -1.0, 1.0, -1.0, -1.0, 1.0, -1.0, 1.0, -1.0, -1.0, 1.0, -1.0,
		1.0, -2.0, 1.0, -1.0, 1.0, -2.0, 1.0, -1.0, 1.0, -2.0, 1.0, -1.0, 1.0, -2.0, 1.0, -1.0, 1.0, -2.0, 1.0, -1.0, 1.0,
		-2.0, 1.0, -1.0, 1.0, -2.0, 1.0, -1.0, 1.0, -2.0, 1.0, -1.0, 1.0, -2.0, 1.0, -1.0, 1.0, -2.0, 1.0, -1.0, 1.0, -2.0,
		1.0, -1.0, 1.0, -1.0, -1.0, 1.0, -1.0, 1.0, -1.0, 1.0, -1.0, 1.0, -1.0, 1.0};

	glp_prob *lp = glp_create_prob();
	glp_set_prob_name(lp, "Addition Chain Bounds");
	glp_set_obj_dir(lp, GLP_MIN);
	glp_add_rows(lp, 45);
	glp_add_cols(lp, 24);
	glp_set_row_bnds(lp, 1, GLP_FX, 1000001.0, 1000001.000000);
	glp_set_row_bnds(lp, 2, GLP_FX, 0.0, 0.0);
	glp_set_row_bnds(lp, 3, GLP_LO, 1.0, 0.0);
	glp_set_row_bnds(lp, 4, GLP_FX, 0.0, 0.0);
	glp_set_row_bnds(lp, 5, GLP_LO, 1.0, 0.0);
	glp_set_row_bnds(lp, 6, GLP_FX, 0.0, 0.0);
	glp_set_row_bnds(lp, 7, GLP_LO, 1.0, 0.0);
	glp_set_row_bnds(lp, 8, GLP_FX, 0.0, 0.0);
	glp_set_row_bnds(lp, 9, GLP_LO, 1.0, 0.0);
	glp_set_row_bnds(lp, 10, GLP_FX, 0.0, 0.0);
	glp_set_row_bnds(lp, 11, GLP_LO, 1.0, 0.0);
	glp_set_row_bnds(lp, 12, GLP_FX, 0.0, 0.0);
	glp_set_row_bnds(lp, 13, GLP_LO, 1.0, 0.0);
	glp_set_row_bnds(lp, 14, GLP_FX, 0.0, 0.0);
	glp_set_row_bnds(lp, 15, GLP_LO, 1.0, 0.0);
	glp_set_row_bnds(lp, 16, GLP_FX, 0.0, 0.0);
	glp_set_row_bnds(lp, 17, GLP_LO, 1.0, 0.0);
	glp_set_row_bnds(lp, 18, GLP_FX, 0.0, 0.0);
	glp_set_row_bnds(lp, 19, GLP_LO, 1.0, 0.0);
	glp_set_row_bnds(lp, 20, GLP_FX, 0.0, 0.0);
	glp_set_row_bnds(lp, 21, GLP_LO, 1.0, 0.0);
	glp_set_row_bnds(lp, 22, GLP_FX, 0.0, 0.0);
	glp_set_row_bnds(lp, 23, GLP_LO, 1.0, 0.0);
	glp_set_row_bnds(lp, 24, GLP_FX, 0.0, 0.0);
	glp_set_row_bnds(lp, 25, GLP_LO, 1.0, 0.0);
	glp_set_row_bnds(lp, 26, GLP_FX, 0.0, 0.0);
	glp_set_row_bnds(lp, 27, GLP_LO, 1.0, 0.0);
	glp_set_row_bnds(lp, 28, GLP_FX, 0.0, 0.0);
	glp_set_row_bnds(lp, 29, GLP_LO, 1.0, 0.0);
	glp_set_row_bnds(lp, 30, GLP_FX, 0.0, 0.0);
	glp_set_row_bnds(lp, 31, GLP_LO, 1.0, 0.0);
	glp_set_row_bnds(lp, 32, GLP_FX, 0.0, 0.0);
	glp_set_row_bnds(lp, 33, GLP_LO, 1.0, 0.0);
	glp_set_row_bnds(lp, 34, GLP_FX, 0.0, 0.0);
	glp_set_row_bnds(lp, 35, GLP_LO, 1.0, 0.0);
	glp_set_row_bnds(lp, 36, GLP_FX, 0.0, 0.0);
	glp_set_row_bnds(lp, 37, GLP_LO, 1.0, 0.0);
	glp_set_row_bnds(lp, 38, GLP_FX, 0.0, 0.0);
	glp_set_row_bnds(lp, 39, GLP_LO, 1.0, 0.0);
	glp_set_row_bnds(lp, 40, GLP_FX, 0.0, 0.0);
	glp_set_row_bnds(lp, 41, GLP_LO, 1.0, 0.0);
	glp_set_row_bnds(lp, 42, GLP_LO, 1.0, 0.0);
	glp_set_row_bnds(lp, 43, GLP_LO, 1.0, 0.0);
	glp_set_row_bnds(lp, 44, GLP_LO, 1.0, 0.0);
	glp_set_row_bnds(lp, 45, GLP_LO, 1.0, 0.0);
	glp_set_col_bnds(lp, 1, GLP_LO, 1.0, 0.0);
	glp_set_col_kind(lp, 1, GLP_IV);
	glp_set_col_bnds(lp, 2, GLP_LO, 1.0, 0.0);
	glp_set_col_kind(lp, 2, GLP_IV);
	glp_set_col_bnds(lp, 3, GLP_LO, 1.0, 0.0);
	glp_set_col_kind(lp, 3, GLP_IV);
	glp_set_col_bnds(lp, 4, GLP_LO, 1.0, 0.0);
	glp_set_col_kind(lp, 4, GLP_IV);
	glp_set_col_bnds(lp, 5, GLP_LO, 1.0, 0.0);
	glp_set_col_kind(lp, 5, GLP_IV);
	glp_set_col_bnds(lp, 6, GLP_LO, 1.0, 0.0);
	glp_set_col_kind(lp, 6, GLP_IV);
	glp_set_col_bnds(lp, 7, GLP_LO, 1.0, 0.0);
	glp_set_col_kind(lp, 7, GLP_IV);
	glp_set_col_bnds(lp, 8, GLP_LO, 1.0, 0.0);
	glp_set_col_kind(lp, 8, GLP_IV);
	glp_set_col_bnds(lp, 9, GLP_LO, 1.0, 0.0);
	glp_set_col_kind(lp, 9, GLP_IV);
	glp_set_col_bnds(lp, 10, GLP_LO, 1.0, 0.0);
	glp_set_col_kind(lp, 10, GLP_IV);
	glp_set_col_bnds(lp, 11, GLP_LO, 1.0, 0.0);
	glp_set_col_kind(lp, 11, GLP_IV);
	glp_set_col_bnds(lp, 12, GLP_LO, 1.0, 0.0);
	glp_set_col_kind(lp, 12, GLP_IV);
	glp_set_col_bnds(lp, 13, GLP_LO, 1.0, 0.0);
	glp_set_col_kind(lp, 13, GLP_IV);
	glp_set_col_bnds(lp, 14, GLP_LO, 1.0, 0.0);
	glp_set_col_kind(lp, 14, GLP_IV);
	glp_set_col_bnds(lp, 15, GLP_LO, 1.0, 0.0);
	glp_set_col_kind(lp, 15, GLP_IV);
	glp_set_col_bnds(lp, 16, GLP_LO, 1.0, 0.0);
	glp_set_col_kind(lp, 16, GLP_IV);
	glp_set_col_bnds(lp, 17, GLP_LO, 1.0, 0.0);
	glp_set_col_kind(lp, 17, GLP_IV);
	glp_set_col_bnds(lp, 18, GLP_LO, 1.0, 0.0);
	glp_set_col_kind(lp, 18, GLP_IV);
	glp_set_col_bnds(lp, 19, GLP_LO, 1.0, 0.0);
	glp_set_col_kind(lp, 19, GLP_IV);
	glp_set_col_bnds(lp, 20, GLP_LO, 1.0, 0.0);
	glp_set_col_kind(lp, 20, GLP_IV);
	glp_set_col_bnds(lp, 21, GLP_LO, 1.0, 0.0);
	glp_set_col_kind(lp, 21, GLP_IV);
	glp_set_col_bnds(lp, 22, GLP_LO, 1.0, 0.0);
	glp_set_col_kind(lp, 22, GLP_IV);
	glp_set_col_bnds(lp, 23, GLP_LO, 1.0, 0.0);
	glp_set_col_kind(lp, 23, GLP_IV);
	glp_set_col_bnds(lp, 24, GLP_LO, 1.0, 0.0);
	glp_set_col_kind(lp, 24, GLP_IV);
	glp_set_obj_coef(lp, 21, 1.0);
	glp_term_out(GLP_OFF);
	int iai = sizeof(ia) / sizeof(ia[0]) - 1;
	glp_load_matrix(lp, iai, ia, ja, ar);
	glp_iocp param;
	glp_init_iocp(&param);
	param.presolve = GLP_ON;
	if (glp_intopt(lp, &param) == 0 && glp_mip_status(lp) == GLP_OPT) {
		double z = glp_mip_obj_val(lp);
	}
}