結果

問題 No.992 最長増加部分列の数え上げ
ユーザー square1001square1001
提出日時 2020-02-14 21:55:25
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 2,319 bytes
コンパイル時間 1,160 ms
コンパイル使用メモリ 87,296 KB
実行使用メモリ 11,488 KB
最終ジャッジ日時 2024-10-06 11:34:51
合計ジャッジ時間 6,822 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,816 KB
testcase_01 AC 2 ms
6,816 KB
testcase_02 AC 1 ms
6,816 KB
testcase_03 AC 2 ms
6,820 KB
testcase_04 WA -
testcase_05 WA -
testcase_06 WA -
testcase_07 WA -
testcase_08 WA -
testcase_09 WA -
testcase_10 WA -
testcase_11 WA -
testcase_12 WA -
testcase_13 WA -
testcase_14 WA -
testcase_15 WA -
testcase_16 WA -
testcase_17 WA -
testcase_18 WA -
testcase_19 WA -
testcase_20 WA -
testcase_21 WA -
testcase_22 WA -
testcase_23 WA -
testcase_24 WA -
testcase_25 WA -
testcase_26 WA -
testcase_27 WA -
testcase_28 WA -
testcase_29 WA -
testcase_30 AC 112 ms
11,136 KB
testcase_31 AC 115 ms
11,136 KB
testcase_32 AC 115 ms
11,136 KB
testcase_33 AC 112 ms
11,360 KB
testcase_34 AC 104 ms
11,136 KB
testcase_35 AC 109 ms
11,008 KB
testcase_36 AC 109 ms
11,136 KB
testcase_37 AC 107 ms
11,136 KB
testcase_38 AC 109 ms
11,008 KB
testcase_39 AC 108 ms
11,136 KB
testcase_40 WA -
testcase_41 WA -
testcase_42 WA -
testcase_43 WA -
testcase_44 WA -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <vector>
#include <iostream>
#include <algorithm>
#include <unordered_map>
using namespace std;
const int mod = 1000000007;
const int inf = 1012345678;
int sz, seg[524295];
void update(int pos, int val) {
	pos += sz;
	seg[pos] = val;
	while (pos > 1) {
		pos >>= 1;
		seg[pos] = max(seg[pos * 2], seg[pos * 2 + 1]);
	}
}
int query(int a, int b, int k, int l, int r) {
	if (r <= a || b <= l) return -inf;
	if (a <= l && r <= b) return seg[k];
	int lc = query(a, b, k * 2, l, (l + r) >> 1);
	int rc = query(a, b, k * 2 + 1, (l + r) >> 1, r);
	return max(lc, rc);
}
int N, bit[200009];
void add(int pos, int val) {
	for (int i = pos + 1; i <= N; i += i & (-i)) {
		bit[i] += val;
	}
}
int sum(int r) {
	int ans = 0;
	for (int i = r; i >= 1; i -= i & (-i)) {
		ans += bit[i];
	}
	return ans;
}
int main() {
	cin.tie(0);
	ios_base::sync_with_stdio(false);
	cin >> N;
	vector<int> A(N);
	for (int i = 0; i < N; ++i) {
		cin >> A[i];
	}
	vector<int> perm(N);
	for (int i = 0; i < N; ++i) {
		perm[i] = i;
	}
	sort(perm.begin(), perm.end(), [&](int i, int j) { return A[i] != A[j] ? A[i] < A[j] : i > j; });
	sz = 1;
	while (sz < N) sz *= 2;
	vector<int> dp(N);
	for (int i = 0; i < N; ++i) {
		dp[perm[i]] = query(0, perm[i], 1, 0, sz) + 1;
		if (dp[perm[i]] <= 0) dp[perm[i]] = 1;
		update(perm[i], dp[perm[i]]);
	}
	vector<int> p2(N);
	for (int i = 0; i < N; ++i) {
		p2[i] = i;
	}
	sort(p2.begin(), p2.end(), [&](int i, int j) { return dp[i] != dp[j] ? dp[i] < dp[j] : i < j; });
	vector<int> p2inv(N);
	for (int i = 0; i < N; ++i) {
		p2inv[p2[i]] = i;
	}
	vector<int> seq(N + 1, inf);
	for (int i = 0; i < N; ++i) {
		seq[dp[i]] = min(seq[dp[i]], p2inv[i]);
	}
	vector<int> ways(N);
	for (int i = 0; i < N; ++i) {
		if (dp[perm[i]] == 1) {
			ways[perm[i]] = 1;
		}
		else {
			int cl = seq[dp[perm[i]] - 1] - 1, cr = seq[dp[perm[i]]];
			while (cr - cl > 1) {
				int cm = (cl + cr) >> 1;
				if (p2[cm] >= perm[i]) cr = cm;
				else cl = cm;
			}
			ways[perm[i]] = sum(cr) - sum(seq[dp[perm[i]] - 1]) + mod;
			if (ways[perm[i]] >= mod) ways[perm[i]] -= mod;
		}
		add(p2inv[perm[i]], ways[perm[i]]);
	}
	int len = *max_element(dp.begin(), dp.end());
	int ans = 0;
	for (int i = 0; i < N; ++i) {
		if (dp[i] == len) {
			ans += ways[i];
			if (ans >= mod) ans -= mod;
		}
	}
	cout << ans << endl;
	return 0;
}
0