Наибольшая возрастающая подпоследовательность за O(n*log(n) с восстановлением ответа

#include <iostream>
#include <vector>
#include <algorithm>
using namespace std;
int main()      
{                                       //В этой задаче мы используем массив, в котором на i-ом месте хранится минимальное число, стоящее на i-ом месте в некоторой возрастающей подпоследовательности. В дальнейшем, под подпоследовательностью мы будем подразумевать возрастающую подпоследовательность
    long long n, a1, k, b, m;
    cin >> n >> a1 >> k >> b >> m;
    vector<long long> a(n), from(m, -1), dp, ans;       //Массив from служит для восстановления ответа. Так как мы делим с остатком на m, то это значит, что числа у нас могут быть от 0 до m - 1, на i-ой позиции стоит число, стоящее перед числом i в некоторой возрастающей подпоследовательности, -1 значит, что это число стоит в начале некоторой подпоследовательности
    a[0] = a1;
    for (int i = 1; i < n; i++)
    {
        a[i] = (k * a[i - 1] + b) % m;
    }
    dp.push_back(a[0]);
    for (int i = 1; i < n; i++)
    {
        auto it = lower_bound(dp.begin(), dp.end(), a[i]);
        if (it == dp.begin())
        {
            dp[it - dp.begin()] = a[i];
        }
        else if (it == dp.end())
        {
            from[a[i]] = dp[dp.size() - 1];
            dp.push_back(a[i]);
        }
        else
        {
            if (a[i] != dp[it - dp.begin()])           //Если мы встречаем число, которое до этого уже встречалось, то мы не обновляем from, потому что множество подпоследовательностей от второго числа - подмножество множества подпоследовательностей от первого числа
            {
                from[a[i]] = dp[it - dp.begin() - 1];
            }
            dp[it - dp.begin()] = a[i];
        }
    }
    long long x = dp[dp.size() - 1];
    while (from[x] != -1)
    {
        ans.push_back(x);
        x = from[x];
    }
    ans.push_back(x);
    for (auto it = ans.rbegin(); it != ans.rend(); it++)
    {
        cout << *it << " ";
    }
}