結果

問題 No.306 さいたま2008
ユーザー YamaKasaYamaKasa
提出日時 2018-09-17 00:30:02
言語 Java21
(openjdk 21)
結果
AC  
実行時間 139 ms / 2,000 ms
コード長 1,726 bytes
コンパイル時間 2,033 ms
コンパイル使用メモリ 77,176 KB
実行使用メモリ 41,892 KB
最終ジャッジ日時 2024-07-18 07:34:30
合計ジャッジ時間 5,446 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 130 ms
41,604 KB
testcase_01 AC 118 ms
41,636 KB
testcase_02 AC 128 ms
41,868 KB
testcase_03 AC 122 ms
40,812 KB
testcase_04 AC 111 ms
41,708 KB
testcase_05 AC 114 ms
41,788 KB
testcase_06 AC 136 ms
41,648 KB
testcase_07 AC 120 ms
41,836 KB
testcase_08 AC 124 ms
41,632 KB
testcase_09 AC 110 ms
41,624 KB
testcase_10 AC 121 ms
41,892 KB
testcase_11 AC 109 ms
41,728 KB
testcase_12 AC 121 ms
40,936 KB
testcase_13 AC 133 ms
41,768 KB
testcase_14 AC 132 ms
41,696 KB
testcase_15 AC 134 ms
41,356 KB
testcase_16 AC 131 ms
41,796 KB
testcase_17 AC 131 ms
41,844 KB
testcase_18 AC 126 ms
41,788 KB
testcase_19 AC 130 ms
41,444 KB
testcase_20 AC 127 ms
41,644 KB
testcase_21 AC 139 ms
41,456 KB
testcase_22 AC 125 ms
41,744 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import java.util.Scanner;

public class Main {
	static double xa;
	static double ya;
	static double xb;
	static double yb;
	public static void main(String[] args) {
		Scanner scan = new Scanner(System.in);
		xa = scan.nextDouble();
		ya = scan.nextDouble();
		xb = scan.nextDouble();
		yb = scan.nextDouble();
		scan.close();
		double yp = 9;
		double r = Math.max(ya, yb);
		double l = Math.min(ya, yb);
		while(Math.abs(r - l) > 0.0000001) {
			double y1 = (2.0 * l + r) / 3.0;
			double y2 = (l + 2.0 * r) / 3.0;
//			System.out.println(l + " " + r);
//			System.out.println(y1 + " " + y2);
//			System.out.println(diff1_func(y1) + " " + diff1_func(y2));
//			System.out.println();
//			if(Math.abs(diff1_func(y1)) < Math.abs(diff1_func(y2))) {
//				r = y2;
//			}else {
//				l = y1;
//			}
			if(func(y1) < func(y2)) {
				r = y2;
			}else {
				l = y1;
			}
		}

		yp = (l + r) / 2.0;
		//System.out.println(yp);
		//System.exit(0);
		for(int i = 0; i < 100; i++) {
			double yp_new = yp - diff1_func(yp) / diff2_func(yp);
			yp = yp_new;
		}
		System.out.println(yp);
	}
	static double func(double yp) {
		return Math.pow(func1(yp), 0.5) +  Math.pow(func2(yp), 0.5);
	}
	static double diff1_func(double yp) {
		return (yp - ya) * Math.pow(xa * xa + (yp - ya) * (yp - ya), -0.5)
				+ (yp - yb) * Math.pow(xb * xb + (yp - yb) * (yp - yb), -0.5);
	}
	static double diff2_func(double yp) {
		return Math.pow(func1(yp), -0.5) - Math.pow(yp - ya, 2) * Math.pow(func1(yp), -1.5) +
				Math.pow(func2(yp), -0.5) - Math.pow(yp - yb, 2)* Math.pow(func2(yp), -1.5);
	}
	static double func1(double yp) {
		return xa * xa + (yp - ya) * (yp - ya);
	}
	static double func2(double yp) {
		return xb * xb + (yp - yb) * (yp - yb);
	}
}
0