結果

問題 No.1064 ∪∩∩ / Cup Cap Cap
コンテスト
ユーザー Shibuyap
提出日時 2020-05-29 22:08:19
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 4 ms / 2,000 ms
コード長 1,040 bytes
コンパイル時間 3,652 ms
コンパイル使用メモリ 192,856 KB
最終ジャッジ日時 2025-01-10 17:25:22
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 36
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <bits/stdc++.h>
#define rep(i,n) for(int i = 0; i < (n); ++i)
#define srep(i,s,t) for (int i = s; i < t; ++i)
#define drep(i,n) for(int i = (n)-1; i >= 0; --i)
using namespace std;
typedef long long int ll;
typedef pair<int,int> P;
#define yn {puts("YES");}else{puts("NO");}
#define MAX_N 200005

int main() {
    ll a, b, c, d;
    cin >> a >> b >> c >> d;
    if((a-c)*(a-c) == 8 * (b - d)){
        cout << "Yes" << endl;
    }else if((a-c)*(a-c) < 8 * (b - d)){
        cout << "No" << endl;
    }else{
        double alpha = (double)(a - c) / 2.0;
        double beta = (double)(b - d) / 2.0;
        double x1 = (-alpha + sqrt(alpha*alpha - 4.0 * beta)) / 2.0;
        double x2 = (-alpha - sqrt(alpha*alpha - 4.0 * beta)) / 2.0;
        double y1 = x1 * x1 + x1 * a + b;
        double y2 = x2 * x2 + x2 * a + b;
        double p = (y2 - y1) / (x2 - x1);
        double q = y1 - p * x1;
        cout << fixed << setprecision(10) << p << ' ';
        cout << fixed << setprecision(10) << q << endl;
    }
    return 0;
}
0