結果

問題 No.425 ジャンケンの必勝法
ユーザー te-shte-sh
提出日時 2017-10-27 16:44:31
言語 D
(dmd 2.109.1)
結果
AC  
実行時間 5 ms / 2,000 ms
コード長 1,120 bytes
コンパイル時間 3,039 ms
コンパイル使用メモリ 152,864 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-06-12 22:09:55
合計ジャッジ時間 3,968 ms
ジャッジサーバーID
(参考情報)
judge5 / judge1
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 5
other AC * 18
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

import std.algorithm, std.conv, std.range, std.stdio, std.string;
void main()
{
auto rd = readln.split.to!(int[]), p = rd[0], q = rd[1];
auto n = 101, a = new real[][](n, n+1);
auto r2 = real(1)/2, r3 = real(1)/3;
foreach (i; 0..n) {
a[i][] = 0;
a[i][i] = 1;
a[i][max(i-q, 0)] -= r2 * i / 100;
a[i][min(i+q, 100)] -= r3 * (100 - i) / 100;
a[i][n] = (r2 * i + r3 * (100 - i)) / 100;
}
gaussElimination(a, n);
writefln("%.7f", (a[p][n] + 1) * r3);
}
auto gaussElimination(T)(ref T[][] a, size_t n)
{
import std.math;
size_t p;
T pmax;
foreach (k; 0..n-1) {
p = k;
pmax = a[k][k].abs;
foreach (i; k+1..n)
if (a[i][k].abs > pmax) {
p = i;
pmax = a[i][k].abs;
}
if (p != k) swap(a[k], a[p]);
auto akk = a[k][k];
foreach (i; k+1..n) {
auto aik = a[i][k];
foreach (j; k..n+1)
a[i][j] -= aik * (a[k][j] / akk);
}
}
a[n-1][n] /= a[n-1][n-1];
foreach_reverse (i; 0..n-1) {
auto ax = T(0);
foreach (j; i+1..n)
ax += a[i][j] * a[j][n];
a[i][n] = (a[i][n] - ax) / a[i][i];
}
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0