結果
問題 | No.813 ユキちゃんの冒険 |
ユーザー |
![]() |
提出日時 | 2019-04-16 02:27:04 |
言語 | C++11 (gcc 13.3.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,097 bytes |
コンパイル時間 | 825 ms |
コンパイル使用メモリ | 64,312 KB |
実行使用メモリ | 11,664 KB |
最終ジャッジ日時 | 2025-01-02 06:41:31 |
合計ジャッジ時間 | 1,531 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 6 WA * 20 |
ソースコード
#include <iostream> #include <math.h> #include <stdio.h> using namespace std; int n; double p, q; double mat[1005][1005]; void swap(int i, int j) { for(int k = 1; k <= n+1; k++){ double t = mat[i][k]; mat[i][k] = mat[j][k]; mat[j][k] = t; } } int main(void) { cin >> n >> p >> q; for(int i = 1; i <= n; i++){ for(int j = 1; j <= n; j++){ if(i==j) mat[i][j] = 1; if(i+1==j) mat[i][j] = -q; if(i==j+1) mat[i][j] = -p; } } mat[1][n+1] = p; for(int i = 1; i <= n; i++){ double max_val = 0; int max_j; for(int j = i; j <= n; j++){ if(fabs(mat[j][i]) > max_val){ max_val = fabs(mat[j][i]); max_j = j; } } if(max_val < 1e-9){ cout << "Non-Regular" << endl; return 0; } swap(i, max_j); double div = mat[i][i]; for(int j = i; j <= i+1; j++) mat[i][j] /= div; mat[i][n+1] /= div; for(int j = 1; j <= i+1; j++){ double coe = mat[j][i]; if(i == j) continue; for(int k = i; k <= i+1; k++){ mat[j][k] -= coe * mat[i][k]; } mat[j][n+1] -= coe * mat[i][n+1]; } } printf("%.11f\n", mat[1][n+1]); return 0; }