結果

問題 No.2222 Respawn
コンテスト
ユーザー stoq
提出日時 2022-05-28 04:01:03
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
MLE  
(最新)
AC  
(最初)
実行時間 -
コード長 1,136 bytes
コンパイル時間 1,990 ms
コンパイル使用メモリ 201,144 KB
最終ジャッジ日時 2025-01-29 16:33:42
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 10 MLE * 20
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;

// f(x) = a + bx
using linear = pair<double, double>;
linear operator+(linear f, linear g) {
  return {f.first + g.first, f.second + g.second};
}
linear operator+(linear f, double a) { return {f.first + a, f.second}; }
linear operator/(linear f, double a) { return {f.first / a, f.second / a}; }
// f(x) = 0
double solve(linear f) { return -f.first / f.second; }
// f(a)
double eval(linear f, double a) { return f.first + f.second * a; }

int main() {
  int n;
  string s;
  cin >> n >> s;
  s.push_back('#');
  vector E(n + 1, vector<double>(n));
  for (int j = n - 2; j >= 0; j--) {
    if (s[j] == '#') continue;
    vector<linear> f(n + 1);
    for (int i = n; i > j; i--) {
      if (i == n - 1) {
        f[i] = {0, 0};
      } else if (s[i] == '.') {
        f[i] = (f[i + 1] + f[i + 2] + E[i][i]) / 3 + 1;
      } else {
        f[i] = {0, 1};
      }
    }
    E[j][j] = solve((f[j + 1] + f[j + 2] + linear{0, -2}) / 3 + 1);
    for (int i = j + 1; i <= n; i++) E[i][j] = eval(f[i], E[j][j]);
  }
  cout << setprecision(20) << setiosflags(ios::fixed);
  cout << E[0][0] << "\n";
}
0