結果

問題 No.309 シャイな人たち (1)
コンテスト
ユーザー cympfh
提出日時 2015-12-02 14:52:52
言語 C++11（廃止可能性あり）
(gcc 13.3.0 + boost 1.89.0)
結果
WA  
実行時間 -
コード長 2,675 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 1,605 ms
コンパイル使用メモリ 172,272 KB
実行使用メモリ 6,948 KB
最終ジャッジ日時 2024-09-14 07:52:06
合計ジャッジ時間 2,172 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 1 WA * 12
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#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;
}

//           ^   >  v   <
int dx[] = { -1, 0, 1,  0 };
int dy[] = {  0, 1, 0, -1 };

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;
  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];

  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))); // for copy

  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;
  }

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

  rep (i, h) {
    rep (j, w) {
      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 = pt[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];
        }
      }
    }
  }

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

  return 0;
}
0