結果
問題 | No.1710 Minimum OR is X |
ユーザー |
![]() |
提出日時 | 2021-10-15 22:03:22 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 117 ms / 2,000 ms |
コード長 | 1,882 bytes |
コンパイル時間 | 3,968 ms |
コンパイル使用メモリ | 183,232 KB |
最終ジャッジ日時 | 2025-01-25 00:49:39 |
ジャッジサーバーID (参考情報) |
judge2 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 30 |
ソースコード
#include <cstdio> #include <cstring> #include <iostream> #include <string> #include <cmath> #include <bitset> #include <vector> #include <map> #include <set> #include <queue> #include <deque> #include <algorithm> #include <complex> #include <unordered_map> #include <unordered_set> #include <random> #include <cassert> #include <fstream> #include <utility> #include <functional> #include <time.h> #include <stack> #include <array> #include <list> #include <atcoder/all> #define popcount __builtin_popcount using namespace std; using namespace atcoder; typedef long long ll; typedef pair<int, int> P; using mint=modint998244353; mint f[2000010], invf[2000010]; void fac(int n){ f[0]=1; for(ll i=1; i<=n; i++) f[i]=f[i-1]*i; invf[n]=f[n].inv(); for(ll i=n-1; i>=0; i--) invf[i]=invf[i+1]*(i+1); } mint comb(int x, int y){ if(!(0<=y && y<=x)) return 0; return f[x]*invf[y]*invf[x-y]; } int main() { int n, m, k; cin>>n>>m>>k; fac(n+m); string x; cin>>x; bool ok[202][202]={}; mint dp[202][202]; auto solve=[&](auto solve, int i, int c)->mint{ if(ok[i][c])return dp[i][c]; ok[i][c]=1; if(i==m-1){ mint s=0; for(int j=0; j<k && j<=c; j++){ s+=comb(c, j); } if(x[i]=='0'){ return dp[i][c]=mint(2).pow(c)-s; }else{ return dp[i][c]=s; } } mint res=0; if(x[i]=='0'){ for(int j=k; j<=c; j++){ res+=solve(solve, i+1, j)*comb(c, j)*(mint(2).pow((m-1-i)*(c-j))); } }else{ mint s=0; for(int j=0; j<k && j<=c; j++){ s+=comb(c, j); } res+=s*solve(solve, i+1, c); } return dp[i][c]=res; }; cout<<solve(solve, 0, n).val()<<endl; return 0; }