Сортировка матрицы методом Шелла

#include <iostream>
using namespace std;
 
// спёрто с Википедии
int increment(long inc[], long size) {
    // inc[] массив, в который заносятся инкременты
    // size размерность этого массива
    int p1, p2, p3, s;
 
    p1 = p2 = p3 = 1;
    s = -1;
    do {// заполняем массив элементов по формуле Роберта Седжвика
        if (++s % 2) {
            inc[s] = 8 * p1 - 6 * p2 + 1;
        }
        else {
            inc[s] = 9 * p1 - 9 * p3 + 1;
            p2 *= 2;
            p3 *= 2;
        }
        p1 *= 2;
        // заполняем массив, пока текущая инкремента хотя бы в 3 раза меньше количества элементов в массиве
    } while (3 * inc[s] < size);
 
    return s > 0 ? --s : 0;// возвращаем количество элементов в массиве
}
 
template<class T>
 
int shellSort(T a[], long size) {
    // inc инкремент, расстояние между элементами сравнения
    // i и j стандартные переменные цикла
    // seq[40] массив, в котором хранятся инкременты
    long inc, i, j, seq[40];
    int s;//количество элементов в массиве seq[40]
 
    int k = 0;
          // вычисление последовательности приращений
    s = increment(seq, size);
    while (s >= 0) {
        //извлекаем из массива очередную инкременту
        inc = seq[s--];
        // сортировка вставками с инкрементами inc
        for (i = inc; i < size; i++) {
            T temp = a[i];
            // сдвигаем элементы до тех пор, пока не дойдем до конца или не упорядочим в нужном порядке
            for (j = i - inc; (j >= 0) && (a[j] > temp); j -= inc) {
                a[j + inc] = a[j];
                k++;
            }
            // после всех сдвигов ставим на место j+inc элемент, который находился на i месте
            a[j + inc] = temp;
        }
    }
    return k;
}
// спёрто с Википедии
 
void my_array() {
    
    int n, i, j, k, a = 0;
 
    cout << "Enter the size of the array:" << endl;
    cin >> n;
    if (n < 1) {
        cout << "Invalid value entered" << endl;
        return;
    }
 
    int **arr = new int*[n];
    int *arr_sort = new int[n*n];
 
    for (i = 0; i < n; i++) {
        arr[i] = new int[n];
        for (j = 0; j < n; j++) {
            arr[i][j] = rand() % 100;
            arr_sort[a] = arr[i][j];
            a++;
            cout << arr[i][j] << '\t';
        }
        cout << endl;
    }
    cout << endl;
 
    k = shellSort(arr_sort, n*n);
 
    a = 0;
    for (i = 0; i < n; i++) {
        for (j = 0; j < n; j++) {
            arr[i][j] = arr_sort[a];
            cout << arr[i][j] << '\t';
            a++;
        }
        cout << endl;
    }
 
    for (int i = 0; i < n; i++) {
        delete [] arr[i];
    }
    delete[] arr;
    delete[] arr_sort;
    
    cout << "K=" << k << endl;
 
}
 
int main() {
        
    my_array();
 
    return 0;
 
}