結果

問題 No.1710 Minimum OR is X
ユーザー keisuke6keisuke6
提出日時 2023-05-21 14:20:16
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 63 ms / 2,000 ms
コード長 1,630 bytes
コンパイル時間 2,205 ms
コンパイル使用メモリ 208,040 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-12-21 19:18:12
合計ジャッジ時間 3,747 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 4 ms
5,248 KB
testcase_01 AC 4 ms
5,248 KB
testcase_02 AC 9 ms
5,248 KB
testcase_03 AC 3 ms
5,248 KB
testcase_04 AC 4 ms
5,248 KB
testcase_05 AC 4 ms
5,248 KB
testcase_06 AC 3 ms
5,248 KB
testcase_07 AC 4 ms
5,248 KB
testcase_08 AC 14 ms
5,248 KB
testcase_09 AC 12 ms
5,248 KB
testcase_10 AC 20 ms
5,248 KB
testcase_11 AC 16 ms
5,248 KB
testcase_12 AC 20 ms
5,248 KB
testcase_13 AC 14 ms
5,248 KB
testcase_14 AC 18 ms
5,248 KB
testcase_15 AC 8 ms
5,248 KB
testcase_16 AC 17 ms
5,248 KB
testcase_17 AC 10 ms
5,248 KB
testcase_18 AC 46 ms
5,248 KB
testcase_19 AC 29 ms
5,248 KB
testcase_20 AC 20 ms
5,248 KB
testcase_21 AC 51 ms
5,248 KB
testcase_22 AC 63 ms
5,248 KB
testcase_23 AC 4 ms
5,248 KB
testcase_24 AC 4 ms
5,248 KB
testcase_25 AC 10 ms
5,248 KB
testcase_26 AC 4 ms
5,248 KB
testcase_27 AC 5 ms
5,248 KB
testcase_28 AC 6 ms
5,248 KB
testcase_29 AC 4 ms
5,248 KB
testcase_30 AC 3 ms
5,248 KB
testcase_31 AC 4 ms
5,248 KB
testcase_32 AC 7 ms
5,248 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