結果

問題 No.475 最終日 - Writerの怠慢
コンテスト
ユーザー tanzaku
提出日時 2016-12-25 01:02:05
言語 Java
(openjdk 23)
結果
WA  
実行時間 -
コード長 1,520 bytes
コンパイル時間 3,792 ms
コンパイル使用メモリ 79,376 KB
実行使用メモリ 61,692 KB
最終ジャッジ日時 2024-12-14 18:53:21
合計ジャッジ時間 11,874 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 3 WA * 12
権限があれば一括ダウンロードができます

ソースコード

diff # 

import java.io.BufferedInputStream;
import java.io.BufferedReader;
import java.io.IOException;
import java.io.InputStreamReader;
import java.util.*;

public class _475 {
	public static void main(String[] args) throws IOException {
		new _475().solve();
	}

	void solve() throws IOException {
		try (final Scanner in = new Scanner(System.in)) {
			int n = in.nextInt();
			long S = in.nextInt();
			int wid = in.nextInt();
			long[] a = new long[n];
			for (int i = 0; i < n; i++) {
				a[i] = in.nextInt();
			}
			a[wid] += 100 * S;
			int[] cnt = new int[n];
			for (int i = 0; i < n; i++) if (i != wid) {
				int low = 0, high = n;
				while (high - low > 1) {
					int mid = (low + high) / 2;
					long x = a[i] + (long)(50 * S + 50 * S / (0.8 + 0.2 * mid));
					if (x > a[wid]) {
						low = mid;
					} else {
						high = mid;
					}
				}
				cnt[low]++;
//				dump(i, low, p);
			}

			for (int i = 1; i < n; i++) {
				logfact[i] = logfact[i-1] + Math.log(i);
			}
			
			double ans = 1;
			int s = 0;
			for (int i = 0; i < n; i++) {
				// sCi/n-1Ci
				double p = s < i ? 0 : prob(n-1, s, i);
				s += cnt[i];
				ans *= p;
//				dump(p, s);
			}
			System.out.printf("%.10f\n", ans);
		}
	}

	double[] logfact = new double[100000];
	double logcomb(int n, int r) {
		return logfact[n] - logfact[r] - logfact[n-r];
	}
	double prob(int n, int s, int i) {
		return Math.exp(logcomb(s, i) - logcomb(n, i));
	}
	
	// for debug
	static void dump(Object... o) {
		System.err.println(Arrays.deepToString(o));
	}
}
0