結果
問題 | No.2669 Generalized Hitting Set |
ユーザー |
![]() |
提出日時 | 2024-03-08 23:57:49 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 1,565 ms / 4,000 ms |
コード長 | 1,315 bytes |
コンパイル時間 | 4,136 ms |
コンパイル使用メモリ | 252,268 KB |
最終ジャッジ日時 | 2025-02-20 03:16:23 |
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 41 |
ソースコード
#include <stdio.h> #include <bits/stdc++.h> #include <atcoder/all> using namespace atcoder; using mint = modint998244353; using namespace std; #define rep(i,n) for (int i = 0; i < (n); ++i) #define Inf32 1000000001 #define Inf64 1000000000000000001 long long comb[25][25]; int dp[1<<24]; int main(){ ios::sync_with_stdio(false); cin.tie(nullptr); comb[0][0] = 1; for(int i=1;i<=24;i++){ for(int j=0;j<=i;j++){ if(j!=0)comb[i][j] += comb[i-1][j-1]; comb[i][j] += comb[i-1][j]; } } int n,m,k; cin>>n>>m>>k; vector<long long> kei(n+1); kei[k] = 1; for(int i=k+1;i<=n;i++){ for(int j=k;j<i;j++){ kei[i] += comb[i][j] * kei[j]; } kei[i] = 1-kei[i]; } rep(i,m){ string s; cin>>s; int ii = 0; rep(j,n){ if(s[j]=='1')ii |= 1<<j; } dp[ii]++; } rep(i,n){ rep(j,1<<n){ if((j>>i)&1){ } else{ dp[j] += dp[j^(1<<i)]; } } } rep(i,1<<n)dp[i] *= kei[__builtin_popcount(i)]; rep(i,n){ rep(j,1<<n){ if((j>>i)&1){ } else{ dp[j^(1<<i)] += dp[j]; } } } vector<int> ans(m+1,Inf32); rep(i,1<<n){ //if(dp[i]<0)cout<<i<<' '<<dp[i]<<endl; int t = __builtin_popcount(i); ans[dp[i]] = min(ans[dp[i]],t); } for(int i=m;i>=1;i--)ans[i-1] = min(ans[i-1],ans[i]); rep(i,m){ cout<<ans[i+1]<<endl; } return 0; }