Question

我需要使用visual c ++中的可调用库编写一个cplex程序。

我需要以这种方式使用numrows和numcols。

由于我的问题刚开始，我只把我的程序开头。我的程序崩溃，我发现在哪里。它发生在循环之后增加数量。看起来它不能走出循环并且std :: cout＆lt;＆lt; NUMROWS＆lt;＆lt; “和”;永远不会出现如果我在循环中写它，我将看到值，但不会在之后。我找不到原因。你知道为什么吗？

由于

//subfunction pk, xik, yjk
int P(int k){
return k-1;
}

int X(int i, int k){
int p, n;
return p + (i-1)*(n-1) + (k-1);
}

int Y(int j, int k){
int p, n, m;
return p + n*p + (j-1)*(m-1) + (k-1);
}

int main (int argc, char **argv){

srand(time(0));
int status = 0;

int project = 4; 
int employee = 5; 
int time = 5; 

int empl [] = {2, 2, 2, 3};
int meet [] = {2, 4, 3, 3};

int n, m, k, p;
n=5;
m=5;
p=4;
int NUMCOLS=n+m;

CPXENVptr env = NULL;
CPXLPptr  lp  = NULL;

double *obj = new double [NUMCOLS]; 

//Objective function
int profit [] = {10, 20, 5, 15};
for (int i=0; i<project; i++){
        obj[i]=profit[i];
}

int      solstat;
double   objval;

double *lb = new double [NUMCOLS];
double *ub = new double [NUMCOLS];
double *x  = new double [NUMCOLS];
int    *matbeg = new int [NUMCOLS];
int    *matcnt = new int [NUMCOLS];
char   *ctype = new char [NUMCOLS];

int **F = new int*[employee]; 

for(int a = 0; a < employee; a++){      
    F[a] = new int [time];
    for(int b = 0; b < time; b++){
        F[a][b]=rand()%2;
        std::cout << F[a][b] << ", ";
    }
}

int NUMROWS=0;  

for(int i=1; i<=n; i++){ //Each xi      
    for(int j=1; j<=m; j++){ //Each yj
        if(F[i][j] ==0){
            NUMROWS++;
        } 
    }
}std::cout << NUMROWS << "and ";

double *rhs = new double [NUMROWS+1];
char   *sense = new char [NUMROWS+1];

for(int i=0; i < NUMROWS; i++) { //Each row
    rhs[i]=1;
    sense[i]='L';
}

int num_entries=-1;
for(int i=0; i < n; i++) {
    for(int j=0; j < m; j++) {
        if(F[i+1][j+1] ==0) { //1st constraint (pk, xik, yjk)
            num_entries++;
            num_entries++;
            num_entries++;
            std::cout << "try ";

        }
    }
}

int NUMNZ = NUMROWS*NUMCOLS;
int    *matind = new int [2*NUMROWS+m];
double *matval  = new double [2*NUMROWS+m];
matrix_entry *M = new matrix_entry  [num_entries+1+m];

num_entries=-1;
int row=-1;
for(int i=0; i < n; i++) {//1st constraint
    for(int j=0; j < m; j++) {
        if(F[i+1][j+1] ==0) {
            row++;
            num_entries++;
            M[num_entries].col=P(project); //pk
            M[num_entries].row=row;
            M[num_entries].val=-1;
            num_entries++;
            M[num_entries].col=X(i,project); //xik
            M[num_entries].row=row;
            M[num_entries].val=1;
            num_entries++;
            M[num_entries].col=Y(j,project); //yik
            M[num_entries].row=row;
            M[num_entries].val=1;
        }
    }
}

Answer 1

您有以下声明。

int employee = 5; 
int time = 5; 
n=5;
m=5;

然后使用 -

初始化F数组

for(int a = 0; a < employee; a++){      
    F[a] = new int [time];
    for(int b = 0; b < time; b++){
         F[a][b]=rand()%2;

然后使用 -

访问它

for(int i=1; i<=n; i++){ //Each xi      
    for(int j=1; j<=m; j++){ //Each yj
        if(F[i][j] ==0){

所以你不受1的限制。你的循环应该是

for(int i=0 i<n; i++){ //Each xi      
    for(int j=0; j<m; j++){ //Each yj

建议：

请使用std::vector或其他一些容器（由Basile Starynkevitch建议）
请更好地命名您的变量（易于理解和记住）
请在循环中使用的变量保持一致，你会忘记边界。

希望这有帮助。

Cplex c ++循环

1 个答案: