結果
| 問題 |
No.2617 容量3のナップザック
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2024-01-31 13:58:03 |
| 言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 167 ms / 2,000 ms |
| コード長 | 1,662 bytes |
| コンパイル時間 | 2,808 ms |
| コンパイル使用メモリ | 222,200 KB |
| 最終ジャッジ日時 | 2025-02-19 00:53:21 |
|
ジャッジサーバーID (参考情報) |
judge2 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 2 |
| other | AC * 40 |
ソースコード
#pragma GCC optimize("Ofast")
#include <bits/stdc++.h>
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<ll> f(2 * n);
f[0] = seed;
for (int i = 1; i < 2 * n; i++) {
f[i] = (f[i - 1] * a + b) % m;
}
vector<ll> v[3];
vector<ll> 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<ll> 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;
}