結果

問題 No.309 シャイな人たち (1)
ユーザー pekempeypekempey
提出日時 2015-12-02 12:03:33
言語 C++11
(gcc 11.4.0)
結果
TLE  
実行時間 -
コード長 1,717 bytes
コンパイル時間 1,195 ms
コンパイル使用メモリ 153,028 KB
実行使用メモリ 12,216 KB
最終ジャッジ日時 2023-10-12 08:53:37
合計ジャッジ時間 11,703 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
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>
#define GET_MACRO(a, b, c, NAME, ...) NAME
#define rep(...) GET_MACRO(__VA_ARGS__, rep3, rep2)(__VA_ARGS__)
#define rep2(i, a) rep3 (i, 0, a)
#define rep3(i, a, b) for (int i = (a); i < (b); i++)
#define repr(...) GET_MACRO(__VA_ARGS__, repr3, repr2)(__VA_ARGS__)
#define repr2(i, a) repr3 (i, 0, a)
#define repr3(i, a, b) for (int i = (b) - 1; i >= (a); i--)
using namespace std;
typedef long long ll;
template<class T1, class T2> inline bool chmin(T1 &a, T2 b) { return b < a && (a = b, true); }
template<class T1, class T2> inline bool chmax(T1 &a, T2 b) { return a < b && (a = b, true); }

double dp[1 << 11];
double P[11][11];
int S[11][11], pts[11];

int main() {
	int r, c;
	cin >> r >> c;
	cin.ignore();
	rep (i, r) rep (j, c) {
		cin >> P[i][j];
		P[i][j] /= 100;
	}
	cin.ignore();
	rep (i, r) rep (j, c) cin >> S[i][j];

	double ans = 0;
	vector<double> curr(1 << c);
	curr[0] = 1;
	queue<int> q;
	rep (i, r) {
		vector<double> next(1 << c);
		rep (s, 1 << c) {
			rep (t, 1 << c) {
				memset(pts, 0, sizeof(pts));
				double pr = 1;
				rep (j, c) {
					if (t & 1 << j) pr *= P[i][j];
					else pr *= 1 - P[i][j];
					if (t & 1 << j) pts[j] += 4 - S[i][j];
					if (s & 1 << j) pts[j]++;
				}
				int u = 0;
				rep (j, c) if (pts[j] >= 4) q.push(j);
				while (!q.empty()) {
					int j = q.front(); q.pop();
					u |= 1 << j;
					if (j - 1 >= 0) {
						pts[j - 1]++;
						if (pts[j - 1] == 4) q.push(j - 1);
					}
					if (j + 1 < c) {
						pts[j + 1]++;
						if (pts[j + 1] == 4) q.push(j + 1);
					}
				}
				next[u] += curr[s] * pr;
			}
		}
		curr = next;
		rep (s, 1 << c) {
			ans += __builtin_popcount(s) * curr[s];
		}
	}
	printf("%.20f\n", ans);
	return 0;
}
0