#include #include #include using namespace std; using ll = long long; using namespace __gnu_pbds; template > using nset = tree; const int mod = 998244353; int N, K; ll V[1010101]; int W[1010101]; int main(void){ ios::sync_with_stdio(false); cin.tie(nullptr); cin >> N >> K; { ll s, a, b, m; cin >> s >> a >> b >> m; for(int i = 0;i < N;i++){ W[i] = s % 3 + 1; s = (s * a + b) % m; } for(int i = 0;i < N;i++){ V[i] = s * W[i]; s = (s * a + b) % m; } } nset> st[3]; for(int i = 0;i < N;i++){ st[W[i] - 1].insert(V[i]); } for(int i = 0;i < N;i++){ st[0].insert(0); } ll ans = 0; for(int i = 0;i < K;i++){ vector a; a.push_back((st[2].empty() ? 0 : *st[2].find_by_order(0))); a.push_back((st[0].empty() || st[1].empty() ? 0 : *st[0].find_by_order(0) + *st[1].find_by_order(0))); a.push_back((st[0].size() < 3 ? 0 : *st[0].find_by_order(0) + *st[0].find_by_order(1) + *st[0].find_by_order(2))); int idx = max_element(a.begin(), a.end()) - a.begin(); if(a[idx] == 0)break; ans += a[idx]; if(idx == 2){ st[0].erase(st[0].begin()); st[0].erase(st[0].begin()); st[0].erase(st[0].begin()); }else if(idx == 1){ st[0].erase(st[0].begin()); st[1].erase(st[1].begin()); }else{ st[2].erase(st[2].begin()); } } cout << ans << endl; return 0; }