結果

問題 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;
}
0