結果

問題 No.2617 容量3のナップザック
ユーザー inksamuraiinksamurai
提出日時 2024-01-27 04:25:16
言語 C++17(gcc12)
(gcc 12.3.0 + boost 1.87.0)
結果
RE  
実行時間 -
コード長 2,037 bytes
コンパイル時間 2,417 ms
コンパイル使用メモリ 211,436 KB
実行使用メモリ 42,736 KB
最終ジャッジ日時 2024-09-28 09:32:48
合計ジャッジ時間 9,672 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 36 RE * 4
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#define int ll
using namespace std;
#define rep(i,n) for(int i=0;i<n;i++)
#define per(i,n) for(int i=n-1;i>=0;i--)
#define rng(i,c,n) for(int i=c;i<n;i++)
#define fi first
#define se second
#define pb push_back
#define all(a) a.begin(), a.end()
#define sz(a) ((int) a.size())
#define vec(...) vector<__VA_ARGS__>
#define _3PawRSC ios::sync_with_stdio(0),cin.tie(0)
typedef long long ll;
using vi=vector<int>;
using pii=pair<int,int>;
void print(){cout<<'\n';}
template<class h,class...t>
void print(const h&v,const t&...u){cout<<v<<' ',print(u...);}

void slv(){
	int n,k;
	cin>>n>>k;
	int seed,a,b,m;
	cin>>seed>>a>>b>>m;
	vi w(n),c(n);
	{
		vi f(2*n);
		f[0]=seed;
		rng(i,1,2*n){
			f[i]=(f[i-1]*a+b)%m;
		}
		rep(i,n){
			w[i]=f[i]%3+1;
		}
		rep(i,n){
			c[i]=w[i]*f[i+n];
		}
	}
	vec(vi) rbts(3);
	rep(i,n){
		rbts[w[i]-1].pb(c[i]);
		// print(w[i],c[i]);
	}
	rep(v,3){
		sort(rbts[v].begin(), rbts[v].end());
	}
	int si=sz(rbts[0]),si1=sz(rbts[1]),si2=sz(rbts[2]);
	vi ps;
	rep(i,si){
		ps.pb((!sz(ps)?0:ps.back())+rbts[0][i]);
	}
	vi ps1;
	rep(i,si1){
		ps1.pb((!sz(ps1)?0:ps1.back())+rbts[1][i]);
	}
	auto range_sum=[&](vi&ps,int l,int r)->int{
		if(l>r) return 0;
		// print(l,r);
		if(l<=0) l=0;
		return ps[r]-(!l?0:ps[l-1]);
	};
	// range_sum(ps,si-1)
	auto af=[&](int x,int take)->int{
		int now=range_sum(ps,si-1-take+1,si-1)+range_sum(ps1,si1-1-take+1,si1-1);
		int leftover=(x-take)*3;
		if(si-1-take>=0){
			now+=range_sum(ps,si-1-take-leftover+1,si-1-take);
		}
		return now;
	};
	// print(af(1,1));
	auto ask=[&](int x)->int{
		int l=0,r=min(x,max(si,si1));
		while(r-l>3){
			int ml=(2*l+r)/3;
			int mr=(l+2*r)/3;
			int fl=af(x,ml);
			int fr=af(x,mr);
			if(fl>fr){
				r=mr;
			}else{
				l=ml;
			}
		}
		int ret=0;
		rng(i,l,r+1){
			// print(x,i,af(x,i));
			ret=max(ret,af(x,i));
		}
		return ret;
	};
	// print(ask(1));
	int ans=0,ad=0;
	rep(i,min(k,si2)+1){
		ans=max(ans,ask(k-i)+ad);
		if(i<si2){
			ad+=rbts[2][si2-1-i];
		}
	}
	print(ans);
}

signed main(){
_3PawRSC;
	slv();
}
0