結果
問題 | No.463 魔法使いのすごろく🎲 |
ユーザー |
|
提出日時 | 2017-06-11 21:18:41 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 3 ms / 2,000 ms |
コード長 | 2,390 bytes |
コンパイル時間 | 1,044 ms |
コンパイル使用メモリ | 111,864 KB |
実行使用メモリ | 6,944 KB |
最終ジャッジ日時 | 2024-09-24 16:35:01 |
合計ジャッジ時間 | 2,189 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 36 |
ソースコード
#define _USE_MATH_DEFINES #include <cstdio> #include <iostream> #include <sstream> #include <fstream> #include <iomanip> #include <algorithm> #include <cmath> #include <complex> #include <string> #include <vector> #include <list> #include <queue> #include <stack> #include <set> #include <map> #include <bitset> #include <numeric> #include <limits> #include <climits> #include <cfloat> #include <functional> #include <iterator> using namespace std; const double EPS = 1.0e-10; int elementaryMatrix(vector<vector<double> >& mat) { const int n = mat.size(); const int m = mat[0].size(); int y = 0; int x = 0; while(y < n && x < m){ int tmp = y; for(int i=y+1; i<n; ++i){ if(abs(mat[i][x]) > abs(mat[tmp][x])) tmp = i; } if(abs(mat[tmp][x]) > EPS){ mat[y].swap(mat[tmp]); for(int i=y+1; i<n; ++i){ for(int j=m-1; j>=x; --j){ mat[i][j] -= mat[y][j] * (mat[i][x] / mat[y][x]); } } ++ y; } ++ x; } mat.resize(y); return y; } bool linearSystem(vector<vector<double> >& mat, vector<double>& x) { int n = mat[0].size() - 1; if(elementaryMatrix(mat) != n || mat[n-1][n-1] == 0) return false; x.resize(n); for(int i=n-1; i>=0; --i){ x[i] = mat[i][n]; for(int j=i+1; j<n; ++j) x[i] -= mat[i][j] * x[j]; x[i] /= mat[i][i]; } return true; } const double INF = DBL_MAX / 10; int main() { int n, m; cin >> n >> m; vector<int> c(n, 0); for(int i=1; i<n-1; ++i) cin >> c[i]; vector<vector<double> > mat(n, vector<double>(n+1, 0.0)); for(int i=0; i<n-1; ++i){ mat[i][i] = 1.0; mat[i][n] = c[i]; for(int j=1; j<=m; ++j){ int a = i + j; if(a > n - 1) a -= (a - (n - 1)) * 2; mat[i][a] -= 1.0 / m; } } mat[n-1][n-1] = 1.0; vector<double> x; linearSystem(mat, x); vector<double> dp(n, 0.0); for(int i=n-2; i>n-2-m; --i) dp[i] = c[i]; for(int i=n-2-m; i>=0; --i){ for(int j=1; j<=m; ++j) dp[i] += dp[i+j] / m; for(int j=1; j<=m; ++j) dp[i] = min(dp[i], x[i+j]); dp[i] += c[i]; } printf("%.10f\n", dp[0]); return 0; }