-
Notifications
You must be signed in to change notification settings - Fork 248
/
LIS.cpp
80 lines (66 loc) · 1.88 KB
/
LIS.cpp
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
/****************************************************************************************************
Finding Longest Increasing Sequence in O(NlogN)
About it: http://e-maxx.ru/algo/longest_increasing_subseq_log
Based on problem http://informatics.mccme.ru/mod/statements/view3.php?id=766&chapterid=1794
****************************************************************************************************/
#include <iostream>
#include <fstream>
#include <cmath>
#include <algorithm>
#include <vector>
#include <set>
#include <map>
#include <stack>
#include <queue>
#include <cstdlib>
#include <cstdio>
#include <string>
#include <cstring>
#include <cassert>
#include <utility>
#include <iomanip>
using namespace std;
const int MAXN = 105000;
const int INF = 1000 * 1000 * 1000;
int n;
int k, b, m;
int a[MAXN];
int d[MAXN];
int ind[MAXN], pr[MAXN];
vector <int> ansv;
int ans = 1;
int main() {
//assert(freopen("input.txt","r",stdin));
//assert(freopen("output.txt","w",stdout));
scanf("%d", &n);
scanf("%d %d %d %d", &a[1], &k, &b, &m);
for (int i = 2; i <= n; i++)
a[i] = (k * a[i - 1] + b) % m;
d[0] = -INF;
for (int i = 1; i <= n; i++)
d[i] = INF;
for (int i = 1; i <= n; i++) {
int pos = upper_bound(d + 1, d + n + 1, a[i]) - d;
if (d[pos - 1] < a[i] && a[i] < d[pos]) {
d[pos] = a[i];
ind[pos] = i;
pr[i] = ind[pos - 1];
if (pos > ans) {
ans = pos;
}
}
}
if (ans == 1) {
printf("%d", a[1]);
}
else {
int cur = ind[ans];
while (cur != 0) {
ansv.push_back(a[cur]);
cur = pr[cur];
}
for (int i = (int) ansv.size() - 1; i >= 0; i--)
printf("%d ", ansv[i]);
}
return 0;
}