結果

問題 No.2794 I Love EDPC-T
ユーザー 箱星
提出日時 2023-02-03 08:44:16
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 1,855 bytes
コンパイル時間 3,457 ms
コンパイル使用メモリ 234,936 KB
最終ジャッジ日時 2025-02-10 08:30:06
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 22 WA * 1 TLE * 8
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <atcoder/convolution>
#include <atcoder/modint>
#include <bits/stdc++.h>

using namespace std;
using namespace atcoder;

using mint = modint998244353;

int main() {
    int n;
    string s;
    cin >> n;
    cin >> s;
    
    for (int i = 0; i < n - 2; i++) {
    	if (s[i] == '<' && s[i + 1] == '<') {
    		cout << 0 << endl;
    		return 0;
    	}
    }
    
    const int fac_size = 600000;
    vector<mint> fac(fac_size);
    vector<mint> inv_fac(fac_size);
    fac[0] = 1;
    inv_fac[0] = 1;
    for (int i = 1; i < fac_size; i++) {
    	fac[i] = fac[i - 1] * i;
    }
    inv_fac[fac_size - 1] = 1 / fac[fac_size - 1];
    for (int i = fac_size - 2; i >= 1; i--) {
    	inv_fac[i] = inv_fac[i + 1] * (i + 1);
    }
    
    vector<int> v;
    int streak = 0;
    for (int i = 0; i < n - 1; i++) {
    	if (s[i] == '>') {
    		streak++;
    	} else if (streak > 0) {
    		if (v.empty() && s[0] == '>') {
    			v.push_back(streak);
    		} else if (streak >= 2) {
    			v.push_back(streak - 1);
    		}
    		streak = 0;
    	}
    }
    if (streak > 0) {
    	v.push_back(streak);
    }
    if (streak == n - 1) {
    	cout << 0 << endl;
    	return 0;
    }
    int cnt = 0;
    for (int i = 0; i < n - 1; i++) {
    	if (s[i] == '<') cnt++;
    }
    vector<mint> f(cnt + 1);
    f[cnt] = 1;
    for (auto e : v) {
    	vector<mint> g;
    	int i = 0;
    	while (true) {
    		if (e - i < i) break;
    		g.push_back(fac[e - i] * inv_fac[i] * inv_fac[e - 2 * i]);
    		i++;
    	}
    	f = convolution(f, g);
    }
    mint ans = 0;
    for (int i = cnt; i < n - 1; i++) {
    	if (i >= f.size()) break;
    	mint y = f[i] * (fac[2 * n - 2 * i] * inv_fac[n - i] * inv_fac[n - i]) / (n - i + 1);
    	if (i % 2 == cnt % 2) {
    		ans += y;
    	} else {
    		ans -= y;
    	}
    }
    cout << ans.val() << endl;
    return 0;
}
0