結果

問題 No.1710 Minimum OR is X
ユーザー pockyny
提出日時 2021-10-16 15:44:50
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 56 ms / 2,000 ms
コード長 1,308 bytes
コンパイル時間 861 ms
コンパイル使用メモリ 70,612 KB
最終ジャッジ日時 2025-01-25 01:35:17
ジャッジサーバーID
(参考情報)
judge3 / judge3
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 30
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#include <iostream>
#include <string>
using namespace std;
typedef long long ll;
const int MX = 100010;
ll f[MX],inv[MX],fi[MX];
constexpr ll mod = 998244353;
void solve(){
inv[1] = 1;
for(int i=2;i<MX;i++){
inv[i] = mod - (mod/i)*inv[mod%i]%mod;
}
f[0] = fi[0] = 1;
for(int i=1;i<MX;i++){
f[i] = f[i-1]*i%mod;
fi[i] = fi[i-1]*inv[i]%mod;
}
}
ll nck(ll n, ll k){
if(n<0 || k<0 || n<k) return 0;
return f[n]*fi[k]%mod*fi[n-k]%mod;
}
ll pw(ll a,ll x){
ll ret = 1;
while(x){
if(x&1) (ret *= a) %= mod;
(a *= a) %= mod; x /= 2;
}
return ret;
}
ll dp[210][210];
int main(){
solve();
ll i,j,l,n,m,k; cin >> n >> m >> k;
string s; cin >> s;
dp[0][n] = 1;
for(i=1;i<=m;i++){
if(s[i - 1]=='0'){
for(j=k;j<=n;j++){
for(l=j;l<=n;l++){
(dp[i][j] += dp[i - 1][l]*nck(l,j)%mod*pw(2,n - l)%mod) %= mod;
}
}
}
if(s[i - 1]=='1'){
for(j=k;j<=n;j++){
for(l=0;l<k;l++){
(dp[i][j] += dp[i - 1][j]*pw(2,n - j)%mod*nck(j,l)) %= mod;
}
}
}
}
ll ans = 0;
for(j=k;j<=n;j++) (ans += dp[m][j]) %= mod;
cout << ans << endl;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0