Codeforces Round #499 (Div. 2) C. Fly

問題ページ
Problem - C - Codeforces

考えたこと

• 燃料xを積んだときに旅行を達成できるか?を考えると単調性がある
• 判定O(N)でできそうなのでにぶたんできそう
• 燃料を10^9までしか使わないことが保証されている
• -1の判定はこれを使えばできそう
• 不等号に=を付け忘れてWA
• -1判定で10^9+1を確かめないといけないのに10^9を使ってWA
• 意味不明なWAを出したけど直したら通った

本番だったらシステスで死んでた、順位表見たらかなり虐殺っぽかった
二分探索しなくて後ろから見ていくだけでいいらしい

ソースコード

#include <bits/stdc++.h>

using namespace std;
using ll = long long;
#define int ll
using VI = vector<int>;
using VVI = vector<VI>;
using PII = pair<int, int>;

#define FOR(i, a, n) for (ll i = (ll)a; i < (ll)n; ++i)
#define REP(i, n) FOR(i, 0, n)
#define ALL(x) x.begin(), x.end()
#define PB push_back

const ll INF = (1LL<<60);
const int MOD = 1000000007;

template <typename T> T &chmin(T &a, const T &b) { return a = min(a, b); }
template <typename T> T &chmax(T &a, const T &b) { return a = max(a, b); }
template <typename T> bool IN(T a, T b, T x) { return a<=x&&x<b; }
template<typename T> T ceil(T a, T b) { return a/b + !!(a%b); }
template<class S,class T>
ostream &operator <<(ostream& out,const pair<S,T>& a){
  out<<'('<<a.first<<','<<a.second<<')';
  return out;
}
template<class T>
ostream &operator <<(ostream& out,const vector<T>& a){
  out<<'[';
  REP(i, a.size()) {out<<a[i];if(i!=a.size()-1)out<<',';}
  out<<']';
  return out;
}

int dx[] = {0, 1, 0, -1}, dy[] = {1, 0, -1, 0};

signed main(void)
{
  cin.tie(0);
  ios::sync_with_stdio(false);

  int n, m;
  cin >> n >> m;
  VI a(n), b(n);
  REP(i, n) cin >> a[i];
  REP(i, n) cin >> b[i];

  auto check = [&](long double mid) -> bool {
    long double x = mid + m;
    long double fuel = 0;
    fuel += x / a[0];
    x -= x / a[0];
    if(x <= 0) return false;
    FOR(i, 1, n) {
      fuel += x / b[i];
      x -= x / b[i];
      if(x <= 0) return false;
      fuel += x / a[i];
      x -= x / a[i];
      if(x <= 0) return false;
    }
    fuel += x / b[0];
    x -= x / b[0];
    if(x <= 0) return false;

    if(mid >= fuel) return true;
    return false;
  };

  if(!check(1e9+1)) {
    cout << -1 << endl;
    return 0;
  }

  long double lb = 0, ub = 1e9+10;
  REP(i, 100) {
    long double mid = (lb+ub)/2;
    // cout << lb << " " << mid << " " << ub << endl;
    if(check(mid)) ub = mid;
    else lb = mid;
  }

  cout << fixed << setprecision(9) << lb << endl;

  return 0;
}