結果
| 問題 |
No.2370 He ate many cakes
|
| ユーザー |
startcpp
|
| 提出日時 | 2023-07-03 21:45:53 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 10 ms / 2,000 ms |
| コード長 | 1,176 bytes |
| コンパイル時間 | 1,229 ms |
| コンパイル使用メモリ | 85,812 KB |
| 実行使用メモリ | 5,376 KB |
| 最終ジャッジ日時 | 2024-07-17 19:38:42 |
| 合計ジャッジ時間 | 2,054 ms |
|
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 4 |
| other | AC * 16 |
ソースコード
//二分探索 + 半分全列挙 O(2^(N/2) * log(MAX)) + QCFirm法
#pragma GCC target("avx2")
#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")
#include <iostream>
#include <algorithm>
#include <vector>
#define rep(i, n) for(i = 0; i < n; i++)
#define int long long
using namespace std;
int n, K, a[30];
void dfs(int l, int r, int sumA, vector<int> &vec) {
if (l == r) {
vec.push_back(sumA);
return;
}
dfs(l + 1, r, sumA, vec);
dfs(l + 1, r, sumA + a[l], vec);
}
//総和がX以上になる選び方がK通り以上存在するか
bool check(int X, vector<int> &vec1, vector<int> &vec2) {
int i, j = (int)vec2.size() - 1;
int ret = 0;
rep(i, vec1.size()) {
for (; j >= 0 && vec1[i] + vec2[j] >= X; j--);
ret += (int)vec2.size() - 1 - j;
}
return ret >= K;
}
signed main() {
int i;
cin >> n >> K;
rep(i, n) cin >> a[i];
vector<int> vec1, vec2;
dfs(0, n / 2, 0, vec1);
dfs(n / 2, n, 0, vec2);
sort(vec1.begin(), vec1.end());
sort(vec2.begin(), vec2.end());
int ok = -1e+11, ng = 1e+11, mid;
while (ng - ok >= 2) {
mid = (ok + ng) / 2;
if (check(mid, vec1, vec2)) ok = mid;
else ng = mid;
}
cout << ok << endl;
return 0;
}
startcpp