結果

問題 No.2633 Subsequence Combination Score
ユーザー 👑 potato167potato167
提出日時 2024-02-13 14:35:09
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 234 ms / 2,000 ms
コード長 1,072 bytes
コンパイル時間 4,813 ms
コンパイル使用メモリ 260,512 KB
最終ジャッジ日時 2025-02-19 05:53:14
ジャッジサーバーID
（参考情報） 		judge3 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 38
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#include<atcoder/all>
using namespace std;
using namespace atcoder;
using mint=modint998244353;
#define rep(i,a,b) for (int i=(int)(a);i<(int)(b);i++)

int main() {
    int N;
    cin >> N;
    vector<int> A(N);
    rep(i, 0, N){
    	cin >> A[i];
    }
    
    int L = A[0];
    vector<mint> fact(L + 1, 1), fact_inv(L + 1, 1);
    rep(i, 1, L + 1) fact[i] = fact[i - 1] * i;
    fact_inv[L] = fact[L].inv();
    for(int i = L; i > 0; i--){
    	fact_inv[i - 1] = fact_inv[i] * i;
    }
    
    mint ans = 0;
    vector<mint> dp(N);
    auto f = [&](auto self, int l, int r) -> void {
    	if(l + 1 == r){
    		dp[l] += fact[A[l]];
    		ans += dp[l] * fact_inv[A[l]];
    		return;
    	}
    	int m = (l + r) / 2;
    	self(self, l, m);
    	vector<mint> X(A[l] - A[m - 1] + 1), Y(A[l] - A[r - 1] + 1);
    	rep(i, l, m) X[A[l] - A[i]] += dp[i];
    	rep(i, 0, Y.size()) Y[i] += fact_inv[i];
    	X = convolution(X, Y);
    	rep(i, m, r) dp[i] += X[A[l] - A[i]];
    	self(self, m, r);
    };
    f(f, 0, N);
    cout << ans.val() << "\n";
}
0