結果

問題 No.309 シャイな人たち (1)
ユーザー tjake
提出日時 2015-12-02 16:15:32
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
WA  
実行時間 -
コード長 2,374 bytes
コンパイル時間 792 ms
コンパイル使用メモリ 75,224 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-09-14 07:56:57
合計ジャッジ時間 19,122 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 4 WA * 9
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<iostream>
#include<string>
#include<vector>
#include<queue>
#include<stack>
#include<map>
#include<set>
#include<algorithm>
#include<functional>
#include<cstdio>
#include<cstdlib>
#include<cmath>
using namespace std;

#define mind(a,b) (a>b?b:a)
#define maxd(a,b) (a>b?a:b)
#define absd(x) (x<0?-(x):x)
#define pow2(x) ((x)*(x))
#define rep(i,n) for(int i=0; i<n; ++i)
#define repr(i,n) for(int i=n-1; i>=0; --i)
#define repl(i,s,n) for(int i=s; i<=n; ++i)
#define replr(i,s,n) for(int i=n; i>=s; --i)
#define repf(i,s,n,j) for(int i=s; i<=n; i+=j)
#define repe(e,obj) for(auto e : obj)

#define SP << " " <<
#define COL << " : " <<
#define COM << ", " <<
#define ARR << " -> " <<
#define PNT(STR) cout << STR << endl
#define POS(X,Y) "(" << X << ", " << Y << ")"
#define DEB(A) " (" << #A << ") " << A
#define DEBREP(i,n,val) for(int i=0; i<n; ++i) cout << val << " "; cout << endl
#define ALL(V) (V).begin(), (V).end()
#define INF 1000000007
#define INFLL 10000000000000000007LL
#define EPS 1e-9

typedef unsigned int uint;
typedef unsigned long ulong;
typedef unsigned long long ull;
typedef long long ll;
typedef long double ld;
typedef pair<int, int> P;
//typedef pair<ll, ll> P;
typedef pair<P, int> PI;
typedef pair<int, P> IP;
typedef pair<P, P> PP;
typedef priority_queue<P, vector<P>, greater<P> > pvqueue;

#define N 13
#define N2 1<<13

int bit_count(int v) {
  int cnt = 0;
  for(;v!=0;v&=v-1) ++cnt;
  return cnt;
}

int main() {
  int r, c;
  cin >> r >> c;
  int p[N][N], s[N][N];
  double dp[N][N2];
  rep(i, r) rep(j, c) cin >> p[i][j];
  rep(i, r) rep(j, c) cin >> s[i][j];
  rep(i, c+1) rep(j, (1<<c)) dp[i][j] = 0.0;

  dp[0][0] = 1.0;
  double ans = 0.0;
  rep(i, r) {
    rep(S, (1<<c)) {
      if(dp[i][S]==0.0) continue;
      rep(T, (1<<c)) {
        double res = dp[i][S];
        rep(j, c) {
          int x = (S>>j)&1, y = (j>0 ? (T>>(j-1))&1 : 0), z = (T>>(j+1))&1;
          if(x+y+z >= s[i][j]) {
            if((T>>j)&1) {
              res *= (double)p[i][j]/100.0;
            } else {
              res *= (double)(100-p[i][j])/100.0;
            }
          } else {
            if((T>>j)&1) {
              res = 0.0;
              break;
            }
          }
        }
        dp[i+1][T] += res;
      }
    }
    rep(T, (1<<c)) {
      ans += dp[i+1][T] * bit_count(T);
    }
  }
  printf("%.10f\n", ans);
  return 0;
}
0