再帰のループ展開その６

(B)に続いて、(C)の置き換えを行う。
(C)は単純な配列のコピーなので、ほとんどそのまま置き換えだけでよい。

(C)


    int l = 0;
    for(int j = 0; j < dim; ++j){
      if(j != i){
        for(int k = 1; k < dim; ++k)subarray[(l * (dim - 1)) + k - 1] = array[(j * dim) + k];
        ++l;
      }
    }

上記が、以下のようになる。

(C)


    int l = 0;
    for(int j = 0; j < d; ++j){
      if(j != i[d - 1]){
      for(int k = 1; k < d; ++k)array[d - 2][(l * (d - 1)) + k - 1] = array[d - 1][(j * d) + k];
        ++l;
      }
    }

array[d - 1] から、array[d - 2]へのコピーだ。
これを埋め込むと以下のようになる。

long[ ][ ] array = new long[dim][ ];
long[ ] return_value = new long[dim];
int[ ] i = new int[dim];
int d = dim;
int flg = 0;
do {
  if(d < 2){
    return_value[d - 1] = array[d - 1][0];
    d++;
    flg = 1;
  }else{
    if(flg == 0){
      if(i[d - 1] == 0){
        array[d - 2] = new long[(d - 1) * (d - 1)];
        long value = 0;
      }
      if(i[d - 1] < d){
        int l = 0;
        for(int j = 0; j < d; ++j){
          if(j != i[d - 1]){
          for(int k = 1; k < d; ++k)array[d - 2][(l * (d - 1)) + k - 1] = array[d - 1][(j * d) + k];
            ++l;
          }
        }
        d--;
        i[d - 1] = 0;
      }
    }else{
      if(i[d - 1] < d){
        (E)
        i[d - 1]++;
        flg = 0;
      }
      if(i[d - 1] == d){
        return_value[d - 1] = value;
        d++;
        flg = 1;
      }
    }
  }
}while(d <= dim);

今日はここまで。