結果

問題 No.1963 Subset Mex
ユーザー 👑 tute7627tute7627
提出日時 2022-04-18 22:33:44
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 1,367 bytes
コンパイル時間 2,066 ms
コンパイル使用メモリ 200,700 KB
最終ジャッジ日時 2025-01-28 19:13:47
ジャッジサーバーID
(参考情報)
judge5 / judge4
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 2 WA * 2
other AC * 20 WA * 6
権限があれば一括ダウンロードができます

ソースコード

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

#include<bits/stdc++.h>
#include<atcoder/modint>
using namespace std;
using modint = atcoder::modint998244353;
const int SZ = 105;
modint dp[SZ][SZ][SZ];
int main(){
int N;
cin >> N;
vector<int> S(N);
for(int i = 0; i < N; i++){
cin >> S[i];
}
sort(S.begin(), S.end());
vector<int> C(SZ);
for(int i = 0; i < N; i++){
C[S[i]]++;
}
memset(dp, 0, sizeof(dp));
dp[SZ - 1][0][0] = 1;
modint ans = 0;
for(int i = SZ - 2; i >= 0; i--){
int require = 0;
vector<int>f(i + 1);
f[max(0, i - 1)] = 1;
for(int j = i - 1; j >= 1; j--){
int use = min(f[j], C[j]);
f[j - 1] = 2 * f[j] - use;
f[j - 1] = min(f[j - 1], 10000);
require -= C[j] - use;
}
require += f[0] - C[0];
require++;
for(int j = 0; j < SZ; j++){
for(int k = 0; k < SZ; k++){
for(int num = 0; num < SZ; num++){
int use = min(j + num, C[i]);
int nj = 2 * j + num - use;
int nk = k + C[i] - use;
if(i == 0 && nj > 0 && j + num <= C[i] + k){
ans += dp[i+1][j][k];
}
if(nj < SZ && nk < SZ){
dp[i][nj][nk] += dp[i+1][j][k];
}
if((nk != 0 || S[0] != i) && nj == 0 && num > 0 && nk >= require - 1){
ans += dp[i+1][j][k];
}
}
}
}
}
cout << ans.val() + 1 << endl;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0