結果
| 問題 | 
                            No.306 さいたま2008
                             | 
                    
| コンテスト | |
| ユーザー | 
                             | 
                    
| 提出日時 | 2017-06-22 23:40:33 | 
| 言語 | D  (dmd 2.109.1)  | 
                    
| 結果 | 
                             
                                AC
                                 
                             
                            
                         | 
                    
| 実行時間 | 2 ms / 2,000 ms | 
| コード長 | 1,144 bytes | 
| コンパイル時間 | 6,115 ms | 
| コンパイル使用メモリ | 139,360 KB | 
| 実行使用メモリ | 6,820 KB | 
| 最終ジャッジ日時 | 2025-02-14 12:14:12 | 
| 合計ジャッジ時間 | 5,820 ms | 
| 
                            ジャッジサーバーID (参考情報)  | 
                        judge5 / judge4 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| other | AC * 26 | 
ソースコード
// tested by Hightail - https://github.com/dj3500/hightail
import std.stdio, std.string, std.conv, std.algorithm;
import std.range, std.array, std.math, std.typecons, std.container, core.bitop;
immutable inf = 10^^9 + 7;
int xa, ya;
int xb, yb;
void main() {
    scan(xa, ya);
    scan(xb, yb);
    if (ya == yb) {
        writeln(ya);
        return;
    }
    real top, btm;
    if (ya > yb) {
        swap(xa, xb);
        swap(ya, yb);
    }
    top = yb;
    btm = ya;
    foreach (i ; 0 .. 1000) {
        real m1 = btm + (top - btm) / 3;
        real m2 = btm + (top - btm) * 2 / 3;
        debug {
            writeln(m1, " ", m2);
        }
        real d1 = sqrt(xa^^2 + (m1 - ya)^^2) + sqrt(xb^^2 + (yb - m1)^^2);
        real d2 = sqrt(xa^^2 + (m2 - ya)^^2) + sqrt(xb^^2 + (yb - m2)^^2);
        if (d1 < d2) {
            top = m2;
        }
        else {
            btm = m1;
        }
    }
    writefln("%.10f", btm);
}
void scan(T...)(ref T args) {
    string[] line = readln.split;
    foreach (ref arg; args) {
        arg = line.front.to!(typeof(arg));
        line.popFront();
    }
    assert(line.empty);
}