結果

問題 No.2794 I Love EDPC-T
ユーザー 箱星箱星
提出日時 2023-02-06 03:02:17
言語 C++17(gcc12)
(gcc 12.3.0 + boost 1.87.0)
結果
AC  
実行時間 1,995 ms / 3,000 ms
コード長 2,515 bytes
コンパイル時間 4,839 ms
コンパイル使用メモリ 255,156 KB
実行使用メモリ 14,344 KB
最終ジャッジ日時 2024-12-20 21:34:21
合計ジャッジ時間 8,922 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 31
権限があれば一括ダウンロードができます

ソースコード

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;
    
    if (s == "<") {
        cout << 1 << endl;
        return 0;
    }
    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++;
    }
    auto c = [](vector<mint> l, vector<mint> r) { return l.size() > r.size(); };
    priority_queue<vector<mint>, vector<vector<mint>>, decltype(c)> que(c);
    vector<mint> f(cnt + 1);
    f[cnt] = 1;
    que.push(f);
    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++;
    	}
    	que.push(g);
    }
    mint con = 1;
    while (que.size() >= 1) {
    	vector<mint> f = que.top();
    	if (f.size() == 1) {
    		que.pop();
    		con *= f[0];
    	} else {
    		break;
    	}
    }
    vector<mint> fcon(1, con);
    que.push(fcon);
    while (que.size() >= 2) {
    	vector<mint> f1 = que.top();
    	que.pop();
    	vector<mint> f2 = que.top();
    	que.pop();
    	que.push(convolution(f1, f2));
    }
    f = que.top();
    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