結果
問題 | No.2669 Generalized Hitting Set |
ユーザー |
|
提出日時 | 2024-03-08 22:52:56 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 1,153 ms / 4,000 ms |
コード長 | 1,103 bytes |
コンパイル時間 | 652 ms |
コンパイル使用メモリ | 67,900 KB |
実行使用メモリ | 70,456 KB |
最終ジャッジ日時 | 2024-09-29 20:21:58 |
合計ジャッジ時間 | 14,247 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 41 |
ソースコード
#include<iostream>#include<cassert>using namespace std;int N,M,K;signed A[1<<24],comb[25][25],coef[25][25],c[25];int ans[3<<17];int main(){ios::sync_with_stdio(false);cin.tie(nullptr);cin>>N>>M>>K;for(int i=0;i<=N;i++){comb[i][0]=comb[i][i]=1;for(int j=1;j<i;j++)comb[i][j]=comb[i-1][j-1]+comb[i-1][j];}for(int k=K;k<=N;k++)coef[k][k]=1;for(int k=N-1;k>=K;k--){for(int i=k+1;i<=N;i++){for(int j=1;j<=N;j++)coef[k][j]-=coef[i][j]*comb[i][k];}}for(int k=K;k<=N;k++){for(int i=k;i>=K;i--)c[k]+=coef[i][k];}for(int i=0;i<M;i++){string s;cin>>s;int S=0;for(int j=0;j<N;j++)if(s[j]=='1')S|=1<<j;A[S]++;}for(int k=0;k<N;k++)for(int i=0;i<1<<N;i++)if(i>>k&1)A[i^1<<k]+=A[i];for(int i=0;i<1<<N;i++){int pc=__builtin_popcount(i);A[i]*=c[pc];}for(int k=0;k<N;k++)for(int i=0;i<1<<N;i++)if(i>>k&1)A[i]+=A[i^1<<k];for(int i=0;i<=M;i++)ans[i]=1e9;for(int i=0;i<1<<N;i++){ans[A[i]]=min(ans[A[i]],__builtin_popcount(i));}for(int i=M;i>=1;i--)ans[i-1]=min(ans[i-1],ans[i]);for(int i=1;i<=M;i++)cout<<ans[i]<<"\n";}