結果

問題 No.2617 容量3のナップザック
ユーザー shobonvipshobonvip
提出日時 2024-01-26 22:21:47
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 199 ms / 2,000 ms
コード長 3,780 bytes
コンパイル時間 4,724 ms
コンパイル使用メモリ 268,056 KB
実行使用メモリ 52,648 KB
最終ジャッジ日時 2024-01-26 22:21:57
合計ジャッジ時間 9,966 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,676 KB
testcase_01 AC 2 ms
6,676 KB
testcase_02 AC 2 ms
6,676 KB
testcase_03 AC 2 ms
6,676 KB
testcase_04 AC 2 ms
6,676 KB
testcase_05 AC 2 ms
6,676 KB
testcase_06 AC 2 ms
6,676 KB
testcase_07 AC 2 ms
6,676 KB
testcase_08 AC 2 ms
6,676 KB
testcase_09 AC 2 ms
6,676 KB
testcase_10 AC 2 ms
6,676 KB
testcase_11 AC 2 ms
6,676 KB
testcase_12 AC 138 ms
38,936 KB
testcase_13 AC 143 ms
45,700 KB
testcase_14 AC 132 ms
37,480 KB
testcase_15 AC 97 ms
31,124 KB
testcase_16 AC 107 ms
28,244 KB
testcase_17 AC 122 ms
36,060 KB
testcase_18 AC 119 ms
36,424 KB
testcase_19 AC 21 ms
9,200 KB
testcase_20 AC 132 ms
42,796 KB
testcase_21 AC 69 ms
20,900 KB
testcase_22 AC 147 ms
46,436 KB
testcase_23 AC 156 ms
46,812 KB
testcase_24 AC 9 ms
6,676 KB
testcase_25 AC 92 ms
30,376 KB
testcase_26 AC 112 ms
29,964 KB
testcase_27 AC 56 ms
17,640 KB
testcase_28 AC 7 ms
6,676 KB
testcase_29 AC 56 ms
19,500 KB
testcase_30 AC 141 ms
43,836 KB
testcase_31 AC 123 ms
35,340 KB
testcase_32 AC 149 ms
51,732 KB
testcase_33 AC 175 ms
52,496 KB
testcase_34 AC 180 ms
51,540 KB
testcase_35 AC 174 ms
51,660 KB
testcase_36 AC 199 ms
52,648 KB
testcase_37 AC 193 ms
51,732 KB
testcase_38 AC 161 ms
52,492 KB
testcase_39 AC 164 ms
50,468 KB
testcase_40 AC 171 ms
52,184 KB
testcase_41 AC 82 ms
50,468 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
using namespace std;

//* ATCODER
#include<atcoder/all>
using namespace atcoder;
typedef modint998244353 mint;
//*/

/* BOOST MULTIPRECISION
#include<boost/multiprecision/cpp_int.hpp>
using namespace boost::multiprecision;
//*/

typedef long long ll;

#define rep(i, s, n) for (int i = (int)(s); i < (int)(n); i++)
#define rrep(i, s, n) for (int i = (int)(n)-1; i >= (int)(s); i--)

template <typename T> bool chmin(T &a, const T &b) {
	if (a <= b) return false;
	a = b;
	return true;
}

template <typename T> bool chmax(T &a, const T &b) {
	if (a >= b) return false;
	a = b;
	return true;
}

template <typename T> T max(vector<T> &a){
	assert(!a.empty());
	T ret = a[0];
	for (int i=0; i<(int)a.size(); i++) chmax(ret, a[i]);
	return ret;
}

template <typename T> T min(vector<T> &a){
	assert(!a.empty());
	T ret = a[0];
	for (int i=0; i<(int)a.size(); i++) chmin(ret, a[i]);
	return ret;
}

template <typename T> T sum(vector<T> &a){
	T ret = 0;
	for (int i=0; i<(int)a.size(); i++) ret += a[i];
	return ret;
}

int main(){
	int n, k; cin >> n >> k;
	ll seed; cin >> seed;
	ll a, b, m; cin >> a >> b >> m;

	vector<ll> v(n), w(n);
	vector<ll> f(n * 2);

	f[0] = seed;
	for (int i=1; i<n * 2; i++){
		f[i] = (a * f[i-1] + b) % m;
	}
	for (int i=0; i<n; i++){
		w[i] = f[i] % 3 + 1;
		v[i] = w[i] * f[n + i];
	}

	vector<ll> w1(0);
	vector<ll> w2(0);
	vector<ll> w3(0);
	rep(i,0,n){
		if (w[i] == 1){
			w1.push_back(v[i]);
		}else if(w[i] == 2){
			w2.push_back(v[i]);
		}else if(w[i] == 3){
			w3.push_back(v[i]);
		}
	}
	sort(w1.begin(), w1.end());
	sort(w2.begin(), w2.end());
	sort(w3.begin(), w3.end());
	reverse(w1.begin(), w1.end());
	reverse(w2.begin(), w2.end());
	reverse(w3.begin(), w3.end());

	int n1 = w1.size();
	int n2 = w2.size();
	int n3 = w3.size();

	vector<ll> rui1(n1 + 1);
	vector<ll> rui2(n2 + 1);
	vector<ll> rui3(n3 + 1);
	rep(i,0,n1) rui1[i+1] = w1[i];
	rep(i,0,n2) rui2[i+1] = w2[i];
	rep(i,0,n3) rui3[i+1] = w3[i];
	
	rep(i,0,n1) rui1[i+1] += rui1[i];
	rep(i,0,n2) rui2[i+1] += rui2[i];
	rep(i,0,n3) rui3[i+1] += rui3[i];

	vector<int> threes(3);
	vector<int> ones(3);
	vector<ll> var(3);

	ll ans = 0;
	rep(l,0,k+1){
		ll tmp = 0;
		tmp += rui2[min(l, n2)];
		while (ones[l % 3] < min(l, n1)){
			var[l % 3] += w1[ones[l % 3]];
			ones[l % 3]++;
		}
		while((ones[l % 3] - min(l, n1)) % 3 != 0){
			ones[l % 3]--;
			var[l % 3] -= w1[ones[l % 3]];
		}
		if(ones[l % 3] > l){
			ones[l % 3]--;
			var[l % 3] -= w1[ones[l % 3]];
			ones[l % 3]--;
			var[l % 3] -= w1[ones[l % 3]];
			ones[l % 3]--;
			var[l % 3] -= w1[ones[l % 3]];
		}
		if(ones[l % 3] > l){
			ones[l % 3]--;
			var[l % 3] -= w1[ones[l % 3]];
			ones[l % 3]--;
			var[l % 3] -= w1[ones[l % 3]];
			ones[l % 3]--;
			var[l % 3] -= w1[ones[l % 3]];
		}
		if(ones[l % 3] > l){
			ones[l % 3]--;
			var[l % 3] -= w1[ones[l % 3]];
			ones[l % 3]--;
			var[l % 3] -= w1[ones[l % 3]];
			ones[l % 3]--;
			var[l % 3] -= w1[ones[l % 3]];
		}
		if (threes[l % 3] >= 1){
			threes[l % 3]--;
			var[l % 3] -= w3[threes[l % 3]];
		}
		if (threes[l % 3] >= 1){
			threes[l % 3]--;
			var[l % 3] -= w3[threes[l % 3]];
		}
		if (threes[l % 3] >= 1){
			threes[l % 3]--;
			var[l % 3] -= w3[threes[l % 3]];
		}
		while ((ones[l % 3] < n1 || threes[l%3] < n3) && (ones[l % 3] - min(l, n1) + 2) / 3 + threes[l % 3] < k - l){
			if (rui1[min(ones[l%3] + 3, n1)] - rui1[min(ones[l%3], n1)]
			<= rui3[min(threes[l%3] + 1, n3)] - rui3[min(threes[l%3], n3)]){
				var[l%3] += rui3[min(threes[l%3] + 1, n3)] - rui3[min(threes[l%3], n3)];
				threes[l%3]++;
				chmin(threes[l%3], n3);
			}else{
				var[l%3] += rui1[min(ones[l%3] + 3, n1)] - rui1[min(ones[l%3], n1)];
				ones[l%3] += 3;
				chmin(ones[l%3], n1);
			}
		}
		chmax(ans, var[l%3] + tmp);
	}
	cout << ans << '\n';
}
0