結果

問題 No.1621 Sequence Inversions
ユーザー 👑 Nachia
提出日時 2021-07-22 21:45:50
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 26 ms / 3,000 ms
コード長 1,041 bytes
コンパイル時間 1,226 ms
コンパイル使用メモリ 84,756 KB
最終ジャッジ日時 2025-01-23 06:27:15
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 26
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <vector>
#include <algorithm>
using namespace std;
using ll = long long;
using ull = unsigned long long;
#define rep(i,n) for(int i=0; i<(n); i++)

const ull MOD = 998244353;

int N,K;
vector<int> A;
vector<int> C;
vector<vector<ull>> dp;

int main(){
  cin >> N >> K;
  A.resize(N);
  rep(i,N) cin >> A[i];
  sort(A.begin(),A.end());
  C.push_back(1);
  rep(i,N-1) if(A[i] == A[i+1]) C.back()++; else C.push_back(1);

  dp.push_back(vector<ull>(K+1,0));
  dp[0][0] = 1;
  for(int c=0; c<N; c++){
    vector<ull> buf(K+1,0);
    for(int i=0; i<=K; i++){
      buf[i] += dp.back()[i];
      if(i+c+1 <= K) buf[i+c+1] += MOD - dp.back()[i];
    }
    rep(i,K) buf[i+1] = (buf[i] + buf[i+1]) % MOD;
    dp.emplace_back(move(buf));
  }
  
  vector<ull> ans = dp.back();
  for(int c : C) if(c != 1){
    for(int k=0; k<=K; k++){
      for(int dk=1; dk<=c*(c-1)/2 && k+dk<=K; dk++){
        ans[k+dk] += (MOD - ans[k]) * dp[c][dk];
        ans[k+dk] %= MOD;
      }
    }
  }

  cout << ans[K] << "\n";
  return 0;
}
0