結果

問題 No.309 シャイな人たち (1)
ユーザー cympfhcympfh
提出日時 2015-12-02 15:36:19
言語 C++11
(gcc 11.4.0)
結果
TLE  
実行時間 -
コード長 3,037 bytes
コンパイル時間 1,710 ms
コンパイル使用メモリ 177,024 KB
実行使用メモリ 59,596 KB
最終ジャッジ日時 2024-09-14 07:53:29
合計ジャッジ時間 12,196 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 TLE -
testcase_01 -- -
testcase_02 -- -
testcase_03 -- -
testcase_04 -- -
testcase_05 -- -
testcase_06 -- -
testcase_07 -- -
testcase_08 -- -
testcase_09 -- -
testcase_10 -- -
testcase_11 -- -
testcase_12 -- -
testcase_13 -- -
testcase_14 -- -
testcase_15 -- -
権限があれば一括ダウンロードができます

ソースコード

diff # 

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

#define rep(i,n) for(int i=0;i<(n);++i)
#define loop for(;;)
#define trace(var) cerr<<">>> "<<#var<<" = "<<var<<endl;
#define all(v) begin(v),end(v)
#define pb push_back
#define inf (1e9)
#define eps (1e-9)
using Integer = long long;
using Real = long double;
const Real PI = acosl(-1);
using P = pair<int, int>;

template<class S, class T> inline
ostream& operator<<(ostream&os, pair<S,T> p) {
  return os << '(' << p.first << ", " << p.second << ')';
}

template<class S, class T, class U> inline
ostream& operator<<(ostream&os, tuple<S,T,U> t) {
  return os << '('
    << get<0>(t) << ", "
    << get<1>(t) << ", "
    << get<2>(t) << ')';
}

template<class T> inline
ostream& operator<<(ostream&os, set<T> v) {
  os << "(set";
  for (T item: v) os << ' ' << item;
  os << ")";
  return os;
}

template<class T> inline
ostream& operator<<(ostream&os, vector<T> v) {
  if (v.size() == 0) { return os << "(empty)"; }
  os << v[0];
  for (int i=1, len=v.size(); i<len; ++i) os << ' ' << v[i];
  return os;
}

template<class T> inline
istream& operator>>(istream&is, vector<T>&v) {
  rep (i, v.size()) is >> v[i];
  return is;
}

using vi = vector<int>;
using vvi = vector<vi>;
using vd = vector<double>;
using vvd = vector<vd>;
using vb = vector<bool>;

int main() {
  cerr << fixed << setprecision(5);

  int h, w;
  cin >> h >> w;

  vector<vvd> pt(5, vvd(h, vd(w,0))); // pt[k][i][j] = Pr [point(i,j) == k]
  vector<vvd> qt(5, vvd(h, vd(w,0)));

  {
    vvi p(h, vi(w, 0));
    vvi s(h, vi(w, 0));
    rep (i, h) rep (j, w) cin >> p[i][j];
    rep (i, h) rep (j, w) cin >> s[i][j];
    rep (i, h) rep (j, w) {
      pt[4-s[i][j]][i][j] = p[i][j] / 100.0;
      pt[0][i][j] = (100 - p[i][j]) / 100.0;
      qt[4-s[i][j]][i][j] = p[i][j] / 100.0;
      qt[0][i][j] = (100 - p[i][j]) / 100.0;
    }
  }

  stack<pair<int, int>> s;

  rep (i, h) rep (j, w) {
    s.push({ i, j });
  }

  vi dx = { 0, -1, 0 };
  vi dy = {-1,  0, 1 };

  while (not s.empty()) {
    int i = s.top().first;
    int j = s.top().second; s.pop();
    if (i < 0 or j < 0 or i>=h or j>=w) continue;
    double memo = qt[4][i][j];

    vd q(4, 0); // q[n] == n人いる確率
    q[0] = 1.0;
    rep (t, 3) {
      int i2 = i + dx[t];
      int j2 = j + dy[t];
      if (i2 < 0 or i2 >= h or j2 < 0 or j2 >= w) continue;
      {
        double p =
          max<double>(
              pt[4][i2][j2],
              qt[4][i2][j2]
              );
        q[4] += q[3] * p;
        q[3] = q[3] * (1-p) + q[2] * p;
        q[2] = q[2] * (1-p) + q[1] * p;
        q[1] = q[1] * (1-p) + q[0] * p;
        q[0] *= (1-p);
      }
    }

    rep (k, 5) qt[k][i][j] = 0;
    rep (k, 5) {
      rep (t, 4) {
        int k2 = k + t; k2 = min(k2, 4);
        qt[k2][i][j] += pt[k][i][j] * q[t];
      }
    }

    if (qt[4][i][j] > memo) {
      s.push({ i+1, j });
      s.push({ i, j-1 });
      s.push({ i, j+1 });
    }
  }

  double ans = 0;
  rep (i, h) rep (j, w) ans += qt[4][i][j];
  cout << ans << endl;

  return 0;
}
0