#1174684 (C++23) No.3604 Min of Max of Div of Sum

• 問題: No.3604 Min of Max of Div of Sum
• ユーザー: kwm_t
• 提出日時: 2026-08-01 11:42:19
• 言語: C++23 (gcc 15.2.0 + boost 1.90.0)
• 結果: AC
• 実行時間: 103 ms / 2,000 ms
• コード長: 3,733 bytes
• コンパイル時間: 2,076 ms
• コンパイル使用メモリ: 336,604 KB
• 実行使用メモリ: 6,272 KB
• 最終ジャッジ日時: 2026-08-05 23:36:36
• 合計ジャッジ時間: 4,513 ms
• ジャッジサーバーID （参考情報）: judge1_0 / judge3_0

ファイルパターン	結果
sample	AC * 3
other	AC * 25

ソースコード

#include <bits/stdc++.h>
//#include <atcoder/all>
using namespace std;
// using namespace atcoder;
// using mint = modint1000000007;
// const int mod = 1000000007;
// using mint = modint998244353;
// const int mod = 998244353;
// const int INF = 1e9;
// const long long LINF = 1e18;
#define rep(i, n) for (int i = 0; i < (n); ++i)
#define rep2(i, l, r) for (int i = (l); i < (r); ++i)
#define rrep(i, n) for (int i = (n)-1; i >= 0; --i)
#define rrep2(i, l, r) for (int i = (r)-1; i >= (l); --i)
#define all(x) (x).begin(), (x).end()
#define allR(x) (x).rbegin(), (x).rend()
#define P pair<int, int>
template<typename A, typename B> inline bool chmax(A& a, const B& b) { if (a < b) { a = b; return true; } return false; }
template<typename A, typename B> inline bool chmin(A& a, const B& b) { if (a > b) { a = b; return true; } return false; }
#ifndef KWM_T_ALGORITHM_SEARCH_BINARYSEARCH_HPP
#define KWM_T_ALGORITHM_SEARCH_BINARYSEARCH_HPP

#include <cstdlib>

namespace kwm_t::algorithm::search {
/**
 * @brief 二分探索 (整数)
 *
 * @details
 *  単調性を持つ述語 f(x) に対して、f(x)=true となる境界を探索する。
 *  初期状態で ok は true、ng は false を満たしている必要がある。
 *
 * Typical Use:
 *  最大/最小の条件を満たす整数を探す
 *
 * Complexity:
 *  Time: O(log |ok-ng|)
 *
 * Template Parameters:
 *  T : 整数型 (int / long long など)
 *  F : bool f(T) を返す関数オブジェクト
 *
 * Parameters:
 *  ok : 条件を満たす値
 *  ng : 条件を満たさない値
 *  f  : 判定関数
 *
 * Return:
 *  条件を満たす境界値
 *
 * Requirements:
 *  f(ok) == true
 *  f(ng) == false
 *
 * Notes:
 *  ok と ng の大小関係は問わない
 *
 * Example:
 *  ```cpp
 *  auto f = [&](long long x){
 *      return x * x <= 100;
 *  };
 *
 *  long long ans = binarySearch(0LL, 11LL, f);
 *  ```
 *
 * Verified:
 *
 */

template <class T, class F>
T binarySearch(T ok, T ng, const F& f) {
	while (std::abs(ok - ng) > 1) {
		T mid = (ok + ng) / 2;
		(f(mid) ? ok : ng) = mid;
	}
	return ok;
}

/**
 * @brief 二分探索 (実数)
 *
 * @details
 *  実数値に対する二分探索。
 *  回数固定で探索を行う。
 *
 * Typical Use:
 *  実数解の近似探索
 *
 * Complexity:
 *  Time: O(counter)
 *
 * Template Parameters:
 *  T : 浮動小数点型 (double / long double)
 *  F : bool f(T) を返す関数オブジェクト
 *
 * Parameters:
 *  ok      : 条件を満たす値
 *  ng      : 条件を満たさない値
 *  f       : 判定関数
 *  counter : 反復回数
 *
 * Return:
 *  条件を満たす近似値
 *
 * Requirements:
 *  f(ok) == true
 *  f(ng) == false
 *
 * Example:
 *  ```cpp
 *  auto f = [&](double x){
 *      return x * x <= 2.0;
 *  };
 *
 *  double ans = binarySearchDouble(0.0, 2.0, f);
 *  ```
 *
 * Verified:
 *
 */

template <class T, class F>
T binarySearchDouble(T ok, T ng, const F& f, int counter = 100) {
	while (counter--) {
		T mid = (ok + ng) / 2;
		(f(mid) ? ok : ng) = mid;
	}
	return ok;
}

} // namespace kwm_t::algorithm::search

#endif // KWM_T_ALGORITHM_SEARCH_BINARYSEARCH_HPP
int main() {
	std::ios::sync_with_stdio(false);
	std::cin.tie(nullptr);
	int n; cin >> n;
	vector<int>a(n), b(n);
	rep(i, n)cin >> a[i];
	rep(i, n)cin >> b[i];
	auto f = [&](double target) {
		vector<double>x(n);
		rep(i, n)x[i] = (double)a[i] - target * b[i];
		vector<double>u(n + 1), v(n + 1);
		rep(i, n)chmax(u[i + 1], u[i] + x[i]);
		rrep(i, n)chmax(v[i], v[i + 1] + x[i]);
		rep(i, n)if (u[i] + x[i] + v[i + 1] < 0)return false;
		return true;
	};
	auto ans = kwm_t::algorithm::search::binarySearchDouble(0.0, (double)1e15, f);
	cout << fixed << setprecision(16) << ans << endl;
	return 0;
}