結果

問題 No.2009 Drunkers' Contest
ユーザー SSRS
提出日時 2022-07-15 22:32:39
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 950 bytes
コンパイル時間 1,817 ms
コンパイル使用メモリ 176,520 KB
実行使用メモリ 18,592 KB
最終ジャッジ日時 2024-06-27 19:17:14
合計ジャッジ時間 6,992 ms
ジャッジサーバーID
(参考情報)
judge1 / judge3
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 4
other AC * 17 WA * 3 TLE * 1 -- * 33
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <bits/stdc++.h>
using namespace std;
const long double INF = 1000000000;
int main(){
  cout << fixed << setprecision(20);
  int N;
  cin >> N;
  vector<int> A(N);
  for (int i = 0; i < N; i++){
    cin >> A[i];
  }
  vector<int> B(N);
  for (int i = 0; i < N; i++){
    cin >> B[i];
  }
  vector<long long> SA(N + 1), SB(N + 1);
  SA[N] = 0;
  SB[N] = 0;
  for (int i = N - 1; i >= 0; i--){
    SA[i] = SA[i + 1] + A[i];
    SB[i] = SB[i + 1] + B[i];
  }
  vector<long double> X;
  for (int i = 0; i < N; i++){
    long double x = sqrt((long double) SA[i] / SB[i]);
    if (x >= 1){
      X.push_back(x);
    }
  }
  X.push_back(1);
  sort(X.begin(), X.end());
  int M = X.size();
  vector<long double> dp(M, 0);
  for (int i = N - 1; i >= 0; i--){
    for (int j = 0; j < M; j++){
      dp[j] += A[i] / X[j] + B[i] * X[j];
    }
    for (int j = M - 2; j >= 0; j--){
      dp[j] = min(dp[j], dp[j + 1]);
    }
  }
  cout << dp[0] << endl;
}
0