結果
問題 | No.425 ジャンケンの必勝法 |
ユーザー |
|
提出日時 | 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 |
ソースコード
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];}}