сортировка изображений

Сравнил попарно все изображения - получается матрица чисел (для удобства сжаты до 0-254) N x N, как теперь упорядочить чтобы похожие были рядом т.е

int r[][] = {
    { 255, 100,  93 },
    { 100, 255, 113 },
    {  93, 113, 255 }};

в этом примере вроде бы правильный порядок of = 1,0,2. n-число изобр. будет меньше 100000, числа в массиве это насколько различны изображения т.е 1 строка 0/0(такие особый случай = 255), 0 и 1, 0 и 2..

1 вариант, функция поиска ищет максимально похожее место вставки; работает совсем не плохо, лучше чем полный перебор. На длинной дистанции результат вообще интересный получился, как preview или нарезка. Лучше бы конечно если 1 куском непрерывно было, а так надо соединять на границах получившиеся сцены (баг или сравнение не точное)

of.push_back(0);
for (int i = 1; i < n; i++)
{
    auto find_bm = [&](int j) {
        std::pair<int, int> m = { 0, 255 };
        for (int i = 0; i < of.size(); i++)
        {
            int _1 = of[i], x = r[j][_1],
                _2, y = x;
            if (i < of.size() - 1)
                _2 = of[i + 1], y = r[j][_2];
            if (m.second > x + y)
                m = { i, x + y };
        }
        return of.begin() + m.first + 1;
    };
    auto it = find_bm(i);
    of.emplace(it, i);
}

создание матрицы, код слишком интегрирован, надеюсь понятно так:

srcc = avs_vfbuf_t(new avs_vfbuf(_srcc, env));
cv::Ptr<cv::img_hash::ImgHashBase> func = cv::img_hash::BlockMeanHash::create();
std::vector<cv::Mat> img_hash(n);
for (int i = 0; i < n; i++)
{
    avs_vf_t srcf = srcc->get(i);
    s_vf_plane srcp = srcf->plane();
    cv::Mat img(srcp.h, srcp.w, CV_8U, srcp.p(), srcp.pitch);
    func->compute(img, img_hash[i]);
}
for (int i = 0; i < n - 1; i++)
{
    for (int j = i + 1; j < n; j++)
        r[i][j] = r[j][i] = std::min(int(func->compare(img_hash[i], img_hash[j])), 254);
}

такой вариант не работает больше чем с ~15 размером, и я не уверен что вообще правильно:

int of[] = { 0, 1, 2 };
std::pair<int, std::vector<int>> bm = { INT32_MAX, of };
do
{
    int _ = 0;
    for (int i = 1; i < n; i++)
        _ += r[of[i - 1]][of[i]];
    if (_ < bm.first)
        bm = { _, of };
} while (std::next_permutation(of.begin(), of.end()));
of = bm.second;