結果

問題 No.3615 Ge Gusser
コンテスト
ユーザー Unbakedbread
提出日時 2026-07-24 18:27:23
言語 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
結果
AC  
実行時間 1,263 ms / 3,000 ms
+ 440µs
コード長 880 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 2,115 ms
コンパイル使用メモリ 340,272 KB
実行使用メモリ 36,480 KB
最終ジャッジ日時 2026-08-06 13:34:50
合計ジャッジ時間 8,067 ms
ジャッジサーバーID
（参考情報） 		judge3_0 / judge2_1
このコードへのチャレンジ
（要ログイン）
サブタスク 配点 結果
サンプル 0 % AC * 3
小課題1 40 % AC * 8
小課題2 40 % AC * 15
小課題3 20 % AC * 27
合計 2.5 * 100% = 250 点
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

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

# pragma GCC target("avx2")
# pragma GCC optimize("O3")
# pragma GCC optimize("unroll-loops")

#define rep(i, n) for(int i = 0; i < (int)(n); ++i)

int main() {
	ios::sync_with_stdio(false);
	cin.tie(nullptr);
	cout.tie(nullptr);
	
	int N, M;
	cin >> N >> M;
	vector<string> S(N);
	rep(i, N) cin >> S[i];
	
	double ans = 0;
	vector<double> dp(1 << N, 0);
	dp[0] = 1;
	
	rep(bit, 1 << N) {
		if(dp[bit] < 1e-13) continue;
		
		vector<bool> vec(M, 1);
		rep(i, N) if(bit >> i & 1) rep(j, M) if(S[i][j] == 'x') vec[j] = 0;
		int cnt = 0;
		rep(i, M) if(vec[i]) cnt += 1;
		if(cnt == 1) {
			ans += __builtin_popcount(bit) * dp[bit];
			continue;
		}
		
		dp[bit] /= (N - __builtin_popcount(bit));
		rep(i, N) {
			if(bit >> i & 1) continue;
			dp[bit | (1 << i)] += dp[bit];
		}
	}
	
	cout << fixed << setprecision(15) << ans << "\n";
}
0