結果

問題 No.1710 Minimum OR is X
ユーザー keisuke6keisuke6
提出日時 2023-05-21 14:20:16
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 62 ms / 2,000 ms
コード長 1,630 bytes
コンパイル時間 2,217 ms
コンパイル使用メモリ 205,648 KB
実行使用メモリ 4,820 KB
最終ジャッジ日時 2023-08-23 15:07:08
合計ジャッジ時間 4,970 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge15
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 4 ms
4,380 KB
testcase_01 AC 4 ms
4,380 KB
testcase_02 AC 8 ms
4,552 KB
testcase_03 AC 3 ms
4,376 KB
testcase_04 AC 3 ms
4,376 KB
testcase_05 AC 4 ms
4,376 KB
testcase_06 AC 3 ms
4,376 KB
testcase_07 AC 3 ms
4,376 KB
testcase_08 AC 14 ms
4,552 KB
testcase_09 AC 12 ms
4,380 KB
testcase_10 AC 20 ms
4,376 KB
testcase_11 AC 15 ms
4,376 KB
testcase_12 AC 19 ms
4,376 KB
testcase_13 AC 15 ms
4,388 KB
testcase_14 AC 17 ms
4,384 KB
testcase_15 AC 8 ms
4,464 KB
testcase_16 AC 18 ms
4,392 KB
testcase_17 AC 10 ms
4,432 KB
testcase_18 AC 44 ms
4,616 KB
testcase_19 AC 29 ms
4,616 KB
testcase_20 AC 20 ms
4,664 KB
testcase_21 AC 52 ms
4,820 KB
testcase_22 AC 62 ms
4,608 KB
testcase_23 AC 3 ms
4,380 KB
testcase_24 AC 3 ms
4,380 KB
testcase_25 AC 10 ms
4,396 KB
testcase_26 AC 3 ms
4,400 KB
testcase_27 AC 4 ms
4,376 KB
testcase_28 AC 6 ms
4,376 KB
testcase_29 AC 4 ms
4,404 KB
testcase_30 AC 3 ms
4,380 KB
testcase_31 AC 5 ms
4,584 KB
testcase_32 AC 6 ms
4,404 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
#define int long long
const int mod = 998244353;
vector<long long> fact, fact_inv, inv;
/*  init_nCk :二項係数のための前処理
    計算量:O(n)
*/
void init_nCk(int SIZE) {
    fact.resize(SIZE + 5);
    fact_inv.resize(SIZE + 5);
    inv.resize(SIZE + 5);
    fact[0] = fact[1] = 1;
    fact_inv[0] = fact_inv[1] = 1;
    inv[1] = 1;
    for (int i = 2; i < SIZE + 5; i++) {
        fact[i] = fact[i - 1] * i % mod;
        inv[i] = mod - inv[mod % i] * (mod / i) % mod;
        fact_inv[i] = fact_inv[i - 1] * inv[i] % mod;
    }
}
/*  nCk :MODでの二項係数を求める(前処理 int_nCk が必要)
    計算量:O(1)
*/
long long nCk(int n, int k) {
    if(n < k || n < 0 || k < 0) return 0;
    return fact[n] * (fact_inv[k] * fact_inv[n - k] % mod) %mod;
}
int int_mod(int a, int x){
	if(x == 0) return 1;
	else if(x%2 == 0){
		int t = int_mod(a,x/2);
		return (t*t)%mod;
	}
	else{
		return (int_mod(a,x-1)*a)%mod;
	} 
}
signed main(){
	init_nCk(50000);
	int N,M,K;
	cin>>N>>M>>K;
	string S;
	cin>>S;
	vector<vector<int>> dp(M+1,vector<int>(N+1,0));
	dp[0][N] = 1;
	for(int i=0;i<M;i++)for(int j=0;j<=N;j++){
		if(S[i] == '1'){
			int x = int_mod(2,N-j), y = 0;
			for(int k=0;k<K;k++){
				y += nCk(j,k);
				y %= mod;
			}
			dp[i+1][j] = ((x*y)%mod)*dp[i][j]%mod;
		}
		else{
			for(int k=j;k>=K;k--){
				dp[i+1][k] += dp[i][j]*int_mod(2,N-j)%mod*nCk(j,k)%mod;
				dp[i+1][k] %= mod;
			}
		}
	}
	/*for(vector<int> s:dp){
		for(int x:s){
			cout<<x<<' ';
		}
		cout<<endl;
	}*/
	int ans = 0;
	for(int x:dp[M]){
		ans += x;
		ans %= mod;
	}
	cout<<ans<<endl;
}
0