#489559 (C++14) No.1064 ∪∩∩ / Cup Cap Cap

• 問題: No.1064 ∪∩∩ / Cup Cap Cap
• ユーザー: chocono2230
• 提出日時: 2020-05-30 11:29:30
• 言語: C++14 (gcc 12.3.0 + boost 1.83.0)
• 結果: AC
• 実行時間: 3 ms / 2,000 ms
• コード長: 1,889 bytes
• コンパイル時間: 1,546 ms
• コンパイル使用メモリ: 166,564 KB
• 実行使用メモリ: 5,248 KB
• 最終ジャッジ日時: 2024-11-07 07:35:21
• 合計ジャッジ時間: 2,780 ms
• ジャッジサーバーID （参考情報）: judge4 / judge5

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	2 ms; 5,248 KB
testcase_01	AC	2 ms; 5,248 KB
testcase_02	AC	2 ms; 5,248 KB
testcase_03	AC	2 ms; 5,248 KB
testcase_04	AC	2 ms; 5,248 KB
testcase_05	AC	3 ms; 5,248 KB
testcase_06	AC	2 ms; 5,248 KB
testcase_07	AC	2 ms; 5,248 KB
testcase_08	AC	2 ms; 5,248 KB
testcase_09	AC	2 ms; 5,248 KB
testcase_10	AC	2 ms; 5,248 KB
testcase_11	AC	2 ms; 5,248 KB
testcase_12	AC	2 ms; 5,248 KB
testcase_13	AC	2 ms; 5,248 KB
testcase_14	AC	2 ms; 5,248 KB
testcase_15	AC	2 ms; 5,248 KB
testcase_16	AC	2 ms; 5,248 KB
testcase_17	AC	2 ms; 5,248 KB
testcase_18	AC	2 ms; 5,248 KB
testcase_19	AC	2 ms; 5,248 KB
testcase_20	AC	2 ms; 5,248 KB
testcase_21	AC	2 ms; 5,248 KB
testcase_22	AC	2 ms; 5,248 KB
testcase_23	AC	2 ms; 5,248 KB
testcase_24	AC	2 ms; 5,248 KB
testcase_25	AC	2 ms; 5,248 KB
testcase_26	AC	2 ms; 5,248 KB
testcase_27	AC	2 ms; 5,248 KB
testcase_28	AC	2 ms; 5,248 KB
testcase_29	AC	2 ms; 5,248 KB
testcase_30	AC	2 ms; 5,248 KB
testcase_31	AC	2 ms; 5,248 KB
testcase_32	AC	2 ms; 5,248 KB
testcase_33	AC	2 ms; 5,248 KB
testcase_34	AC	2 ms; 5,248 KB
testcase_35	AC	2 ms; 5,248 KB
testcase_36	AC	2 ms; 5,248 KB
testcase_37	AC	2 ms; 5,248 KB
testcase_38	AC	2 ms; 5,248 KB
testcase_39	AC	2 ms; 5,248 KB

ソースコード

#include <bits/stdc++.h>
#define rep(i,n) for(int i = 0; i < (int)(n); i++)
#define rrep(ri,n) for(int ri = (int)(n-1); ri >= 0; ri--)
#define rep2(i,x,n) for(int i = (int)(x); i < (int)(n); i++)
#define repit(itr,x) for(auto itr = x.begin(); itr != x.end(); itr++)
#define rrepit(ritr,x) for(auto ritr = x.rbegin(); ritr != x.rend(); ritr++)
#define ALL(x) x.begin(), x.end()
using ll = long long;
using namespace std;

const long double eps = 1e-12;

long double fc(long double n, long double a, long double b, long double c, long double d){
  long double ret = -2*n*n + (c - a)*n + (d-b);
  return ret;
}

int main(){
  int ia, ib, ic, id;
  long double  a, b, c, d;
  cin >> ia >> ib >> ic >> id;
  a = ia; b = ib; c = ic; d = id;
  long double D = (c-a)*(c-a) + 8*(d-b);
  if(D < 0){
    cout << "No" << endl;
    return 0;
  }else if(D == 0){
    cout << "Yes" << endl;
    return 0;
  }
  // long double r = (c-a) * (c-a) + 8*(d-b);
  // long double sq = mysqrtl(r);
  // long double xa = ((c - a) + sq)/4;
  // long double xb = ((c - a) - sq)/4;
  // long double ya = xa*xa + a*xa + b;
  // long double yb = xb*xb + a*xb + b;
  // long double p = (yb-ya)/(xb-xa);
  // long double q = -p * xa + ya;
  long double m = (c-a)/4;
  // cerr << m << endl;
  long double lim = -1e30, l = m;
  while(l-lim > eps){
    long double mid = (l+lim)/2;
    long double fa = fc(mid, a, b, c, d);
    if(fa < 0) lim = mid;
    else l = mid;
  }
  lim = 1e30;
  long double r = m;
  while(lim-r > eps){
    long double mid = (r+lim)/2;
    long double fa = fc(mid, a, b, c, d);
    if(fa >= 0) r = mid;
    else lim = mid;
  }
  // cout << fixed << setprecision(15) << l<< " " << r << endl;
  long double ya = l*l + a*l + b;
  long double yb = r*r + a*r + b;
  long double p = (yb-ya)/(r-l);
  long double q = -p * l + ya;
  cout << fixed << setprecision(15) << p << " " << q << endl;
  return 0;
}