結果
| 問題 |
No.67 よくある棒を切る問題 (1)
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2022-05-06 17:52:15 |
| 言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
TLE
|
| 実行時間 | - |
| コード長 | 1,887 bytes |
| コンパイル時間 | 1,927 ms |
| コンパイル使用メモリ | 194,388 KB |
| 実行使用メモリ | 11,680 KB |
| 最終ジャッジ日時 | 2025-03-03 11:58:19 |
| 合計ジャッジ時間 | 14,607 ms |
|
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| other | TLE * 1 -- * 29 |
ソースコード
// 相対誤差に応じて二分探索を終了できるように改良
// とりあえず、doulbe 型のみの対応
#include <bits/stdc++.h>
using namespace std;
using ll = long long;
const int INF = (int)1e9 + 1001010;
const ll llINF = (ll)4e18 + 11000010;
#define ALL(x) x.begin(),x.end()
#define RALL(x) x.rbegin(),x.rend()
ll ceil(ll a, ll b){return (a+b-1) / b;};
// ================================== ここまでテンプレ ==================================
// とある条件を満たす区間の境界を見つける
// 探索区間は [ok, ng) または (ng, ok]
// ok はつねに「とある条件」を満たす
// ng はつねに「とある条件」を満たさない
// 「とある条件」を満たすかどうかは、judge関数によって求められる
template<class Judgement>
double Binary_Search(double ok, double ng, Judgement judge, double tolerance = 1){
// assert(!is_integral_v<double> || tolerance == 1);
// while((fabs(ok - ng) > tolerance) && (fabs(ok - ng) > tolerance * fabs((ok + ng) * 0.5))){
while((fabs(ok - ng) > tolerance)){
double mid;
// if(is_integral_v<double>){
// if(ok > 0 ^ ng > 0) mid = ok + (ng - ok)/2;
// else mid = (ok + ng)/2;
// }
// else mid = (ok + ng) * 0.5;
mid = (ok + ng) * 0.5;
if(judge(mid)) ok = mid;
else ng = mid;
}
return ok;
};
int main(){
int n; cin >> n;
vector<double> l(n);
for(int i = 0; i < n; i++) cin >> l[i];
ll k; cin >> k;
auto judge = [&](double mid) -> bool{
if(mid <= 0) return true;
ll count = 0;
for(int i = 0; i < n; i++){
count += floor(l[i] / mid);
if(count >= k) return true;
}
return false;
};
double ans = Binary_Search(-1.0, 1e12, judge, 1e-10);
printf("%.12lf\n", ans);
return 0;
}