結果

問題 No.3615 Ge Gusser
コンテスト
ユーザー D M
提出日時 2026-08-19 15:36:03
言語 C++23
(gcc 15.2.0 + boost 1.90.0)
コンパイル:
g++-15 -O2 -lm -std=c++23 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
WA  
実行時間 -
コード長 2,020 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 4,819 ms
コンパイル使用メモリ 379,216 KB
実行使用メモリ 1,003,904 KB
最終ジャッジ日時 2026-08-19 15:36:19
合計ジャッジ時間 16,033 ms
ジャッジサーバーID
（参考情報） 		judge3_0 / judge2_0
このコードへのチャレンジ
（要ログイン）
サブタスク 配点 結果
サンプル 0 % AC * 3
小課題1 40 % AC * 8
小課題2 40 % AC * 15
小課題3 20 % AC * 24 WA * 3
合計 2.5 * 80% = 200 点
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <bits/stdc++.h>
using namespace std;
#define rep(i, x, limit) for (long long i = (long long)x; i < (long long)limit; i++)
#define REP(i, x, limit) for (long long i = (long long)x; i <= (long long)limit; i++)
#define all(x) x.begin(), x.end()
#define rall(x) x.rbegin(), x.rend()
#define el '\n'
#define spa " "
#define Yes cout << "Yes" << el
#define No cout << "No" << el
#define YES cout << "YES" << el
#define NO cout << "NO" << el
#define inp(x) for(auto &i:x)cin>>i
#define eps (1e-10)
#define Equals(a,b) (fabs((a) - (b)) < eps )
#define debug(x) cerr << #x << " = " << x << el
using ll = long long;
using ull = unsigned long long;
using i128 = __int128_t;
using pii = pair<int,int>;
using pll = pair<ll,ll>;
using vi = vector<int>;
using vl = vector<ll>;
using vvl = vector<vector<ll>>;
using vs = vector<string>;
using vb = vector<bool>;
const double pi = 3.141592653589793238;
const int inf = 1073741823;
const ll infl = 1LL << 60;
const string ABC = "ABCDEFGHIJKLMNOPQRSTUVWXYZ";
const string abc = "abcdefghijklmnopqrstuvwxyz";
const ll MOD = 998244353;
#include<atcoder/all>
using namespace atcoder;
using mint = modint998244353;
using vm = vector<mint>;

int main() {
  ios::sync_with_stdio(false);
  cin.tie(nullptr);
  int n,m;cin>>n>>m;
  vs S(n);inp(S);
  vector<vi> J(1LL<<n,vi(m,0));
  rep(i,1,1<<n){
    rep(j,0,n){
      if((i>>j)&1){
        ll pre=i^(1<<j);
        rep(k,0,m){
          J[i][k]=J[pre][k]|(S[j][k]=='x');
        }
        break;
      }
    }
  }
  vb B(1<<n,false);
  rep(i,0,1<<n){
    int cnt=0;
    rep(j,0,m){
      if(J[i][j]==0)cnt++;
    }
    if(cnt==1)B[i]=true;
  }
  vector<double> P(1<<n);
  P[0]=1;
  vi CC(1<<n);
  rep(i,1,1<<n){
    int cnt=0;
    rep(j,0,m){
      if((i>>j)&1)cnt++;
    }
    CC[i]=cnt;
    rep(j,0,m){
      if((i>>j)&1){
        ll pre=i^(1<<j);
        if(B[pre]==false){
        P[i]+=P[pre]/(n-cnt+1);}
      }
    }
  }
  double ans=0;
  rep(i,0,1<<n){
    if(B[i])ans+=P[i]*CC[i];
  }
  cout<<fixed<<setprecision(10)<<ans<<el;
}
0