#include #include #define rep(i,b) for(int i=0;i=0;i--) #define rep1(i,b) for(int i=1;i=x;i--) #define fore(i,a) for(auto& i:a) #define rng(x) (x).begin(), (x).end() #define rrng(x) (x).rbegin(), (x).rend() #define sz(x) ((int)(x).size()) #define pb push_back #define fi first #define se second #define pcnt __builtin_popcountll using namespace std; using namespace atcoder; using ll = long long; using ld = long double; template using mpq = priority_queue, greater>; template bool chmax(T &a, const T &b) { if (a bool chmin(T &a, const T &b) { if (b ll sumv(const vector&a){ll res(0);for(auto&&x:a)res+=x;return res;} bool yn(bool a) { if(a) {cout << "Yes" << endl; return true;} else {cout << "No" << endl; return false;}} #define retval(x) {cout << #x << endl; return;} #define cout2(x,y) cout << x << " " << y << endl; #define coutp(p) cout << p.fi << " " << p.se << endl; #define out cout << ans << endl; #define outd cout << fixed << setprecision(20) << ans << endl; #define outm cout << ans.val() << endl; #define outv fore(yans , ans) cout << yans << "\n"; #define outdv fore(yans , ans) cout << yans.val() << "\n"; #define assertmle(x) if (!(x)) {vi v(3e8);} #define asserttle(x) if (!(x)) {while(1){}} #define coutv(v) {fore(vy , v) {cout << vy << " ";} cout << endl;} #define coutv2(v) fore(vy , v) cout << vy << "\n"; #define coutvm(v) {fore(vy , v) {cout << vy.val() << " ";} cout << endl;} #define coutvm2(v) fore(vy , v) cout << vy.val() << "\n"; using pll = pair;using pil = pair;using pli = pair;using pii = pair;using pdd = pair; using vi = vector;using vd = vector;using vl = vector;using vs = vector;using vb = vector; using vpii = vector;using vpli = vector;using vpll = vector;using vpil = vector; using vvi = vector>;using vvl = vector>;using vvs = vector>;using vvb = vector>; using vvpii = vector>;using vvpli = vector>;using vvpll = vector;using vvpil = vector; using mint = modint998244353; //using mint = modint1000000007; //using mint = dynamic_modint<0>; using vm = vector; using vvm = vector>; vector dx={1,0,-1,0,1,1,-1,-1},dy={0,1,0,-1,1,-1,1,-1}; ll gcd(ll a, ll b) { return a?gcd(b%a,a):b;} ll lcm(ll a, ll b) { return a/gcd(a,b)*b;} #define yes {cout <<"Yes"<>n>>k>>s>>a>>b>>m; vl f(n*2); f[0] = s; rep1(i,n*2) f[i] = (f[i-1] * a + b) % m; vl w(n),v(n); rep(i,n) w[i] = f[i] % 3 + 1; rep(i,n) v[i] = w[i] * f[i+n]; vvl x(3); rep(i,n) x[w[i]-1].pb(v[i]); while(sz(x[0]) < 3*n) x[0].pb(0); while(sz(x[1]) < n) x[1].pb(0); while(sz(x[2]) < n) x[2].pb(0); rep(i,3) sort(rrng(x[i])); vvl sum(3); rep(i,3){ sum[i] = vl(sz(x[i])); rep(j,sz(x[i])){ sum[i][j] += x[i][j]; if (j) sum[i][j] += sum[i][j-1]; } } ll ans = sum[2][k-1]; rep(i,k){ ll tmp = 0; if (i) tmp = sum[2][i-1]; int l = 0, u = k-i-1; auto check = [&](int j)->bool{ int cnt1 = j + ((k-i)-j) * 3; int id1 = cnt1 - 1; int id2 = cnt1 - 2; int id3 = j; ll dif = x[1][id3] - x[0][id1] - x[0][id2]; return (dif>=0); }; if (check(u)){ tmp += sum[1][k-i-1]; tmp += sum[0][k-i-1]; chmax(ans, tmp); continue; }else if(!check(l)){ tmp += sum[0][(k-i)*3-1]; chmax(ans, tmp); continue; } while(u-l>1){ int mid = (l+u) / 2; if (check(mid)) l = mid; else u = mid; } tmp += sum[1][l] + sum[0][(k-i) + (k-i-(l+1))*2 - 1]; chmax(ans, tmp); } out; return; } int main(){ ios::sync_with_stdio(false); cin.tie(0); int t = 1; //cin>>t; rep(i,t){ solve(); } return 0; }