#pragma GCC optimize("Ofast") #include using namespace std; typedef long long int ll; typedef unsigned long long int ull; mt19937_64 rng(chrono::steady_clock::now().time_since_epoch().count()); ll myRand(ll B) { return (ull)rng() % B; } int main() { ios::sync_with_stdio(false); cin.tie(nullptr); int n, k; cin >> n >> k; ll seed, a, b, m; cin >> seed >> a >> b >> m; vector f(2 * n); f[0] = seed; for (int i = 1; i < 2 * n; i++) { f[i] = (f[i - 1] * a + b) % m; } vector v[3]; vector rs[3]; for (int i = 0; i < n; i++) { int w = f[i] % 3 + 1; v[w - 1].push_back(w * f[n + i]); } for (int i = 0; i < 3; i++) { sort(v[i].rbegin(), v[i].rend()); rs[i].resize(v[i].size() + 1); for (int j = 0; j < v[i].size(); j++) { rs[i][j + 1] = rs[i][j] + v[i][j]; } } vector dp(n + 1); int r = 0; for (int i = 1; i <= n; i++) { auto cal = [&](int t) -> ll { ll res = rs[1][min((int)v[1].size(), t)] + rs[0][min((int)v[0].size(), 3 * i - t * 2)]; return res; }; while (r < v[1].size() and r < i) { ll pre = cal(r); ll nxt = cal(r + 1); if (pre <= nxt) { r += 1; } else { break; } } dp[i] = cal(r); } ll res = 0; ll sum = 0; for (int i = 0; i <= v[2].size(); i++) { if (i > k) break; res = max(res, sum + dp[k - i]); if (i == v[2].size()) break; sum += v[2][i]; } cout << res << endl; }