結果

問題 No.1064 ∪∩∩ / Cup Cap Cap
コンテスト
ユーザー mkmkmkmk
提出日時 2020-05-29 22:10:54
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 2 ms / 2,000 ms
コード長 2,154 bytes
コンパイル時間 1,703 ms
コンパイル使用メモリ 168,932 KB
実行使用メモリ 6,820 KB
最終ジャッジ日時 2024-11-06 05:04:44
合計ジャッジ時間 2,741 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 36
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#include <fstream>
using namespace std;
typedef long long lint;
#define rep(i,n) for(lint (i)=0;(i)<(n);(i)++)
#define repp(i,m,n) for(lint (i)=(m);(i)<(n);(i)++)
#define repm(i,n) for(lint (i)=(n-1);(i)>=0;(i)--)
#define INF (1ll<<60)
#define all(x) (x).begin(),(x).end()
const lint MOD =1000000007;
const lint MAX = 1000000;
using Graph =vector<vector<lint>>;
typedef pair<lint,lint> P;
typedef map<lint,lint> M;
 
lint fac[MAX], finv[MAX], inv[MAX];
 
void COMinit() 
{
    fac[0] = fac[1] = 1;
    finv[0] = finv[1] = 1;
    inv[1] = 1;
    for (lint i = 2; i < MAX; i++)
    {
        fac[i] = fac[i - 1] * i % MOD;
        inv[i] = MOD - inv[MOD % i] * (MOD / i) % MOD;
        finv[i] = finv[i - 1] * inv[i] % MOD;
    }
}
 
long long COM(lint n, lint k)
{
    if (n < k)
        return 0;
    if (n < 0 || k < 0)
        return 0;
    return fac[n] * (finv[k] * finv[n - k] % MOD) % MOD;
}
 
lint primary(lint num)
{
    if (num < 2) return 0;
    else if (num == 2) return 1;
    else if (num % 2 == 0) return 0;
 
    double sqrtNum = sqrt(num);
    for (int i = 3; i <= sqrtNum; i += 2)
    {
        if (num % i == 0)
        {
            return 0;
        }
    }
 
    return 1;
}
   long long modpow(long long a, long long n, long long mod) {
    long long res = 1;
    while (n > 0) {
        if (n & 1) res = res * a % mod;
        a = a * a % mod;
        n >>= 1;
    }
    return res;
}
    lint lcm(lint a,lint b){
        return a/__gcd(a,b)*b;
    }
     lint gcd(lint a,lint b){
        return __gcd(a,b);
    } 
    
 
  int main(){
     lint a,b,c,d;
     cin>>a>>b>>c>>d;
     lint s=a-c;
     lint z=b-d;
     if((s*s-8*z)<0)cout<<"No"<<endl;
     else if((s*s-8*z)==0)cout<<"Yes";
     else{
         long double x,x1,y,y1;
         x=(long double)(-s+sqrtl(s*s-8*z))/4;
         x1=(long double)(-s-sqrtl(s*s-8*z))/4;
         y=x*x+a*x+b;
         y1=x1*x1+a*x1+b;
        long double p=(y1-y)/(x1-x);
        long double q=y-(y1-y)/(x1-x)*x;
        if(p==0)p=0;
        if(q==0)q=0;
        printf("%.12Lf %.12Lf\n",p,q);
        
     }
          }
0