結果

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

ソースコード

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], SZ);
      require -= C[j] - use;
    }
    require += f[0];

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