結果
問題 | No.992 最長増加部分列の数え上げ |
ユーザー | 🍮かんプリン |
提出日時 | 2020-10-02 03:33:04 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 223 ms / 2,000 ms |
コード長 | 4,780 bytes |
コンパイル時間 | 1,650 ms |
コンパイル使用メモリ | 187,212 KB |
実行使用メモリ | 38,468 KB |
最終ジャッジ日時 | 2024-07-07 16:55:47 |
合計ジャッジ時間 | 10,947 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 1 ms
5,248 KB |
testcase_01 | AC | 2 ms
5,376 KB |
testcase_02 | AC | 2 ms
5,376 KB |
testcase_03 | AC | 2 ms
5,376 KB |
testcase_04 | AC | 78 ms
16,968 KB |
testcase_05 | AC | 55 ms
13,364 KB |
testcase_06 | AC | 93 ms
20,004 KB |
testcase_07 | AC | 71 ms
15,520 KB |
testcase_08 | AC | 37 ms
9,856 KB |
testcase_09 | AC | 76 ms
15,812 KB |
testcase_10 | AC | 101 ms
20,028 KB |
testcase_11 | AC | 123 ms
24,448 KB |
testcase_12 | AC | 24 ms
7,680 KB |
testcase_13 | AC | 66 ms
15,060 KB |
testcase_14 | AC | 63 ms
15,176 KB |
testcase_15 | AC | 24 ms
7,680 KB |
testcase_16 | AC | 195 ms
35,136 KB |
testcase_17 | AC | 35 ms
9,856 KB |
testcase_18 | AC | 65 ms
14,768 KB |
testcase_19 | AC | 125 ms
24,112 KB |
testcase_20 | AC | 218 ms
38,452 KB |
testcase_21 | AC | 219 ms
38,384 KB |
testcase_22 | AC | 223 ms
38,436 KB |
testcase_23 | AC | 217 ms
38,332 KB |
testcase_24 | AC | 215 ms
38,384 KB |
testcase_25 | AC | 218 ms
38,456 KB |
testcase_26 | AC | 213 ms
38,464 KB |
testcase_27 | AC | 212 ms
38,384 KB |
testcase_28 | AC | 213 ms
38,360 KB |
testcase_29 | AC | 215 ms
38,432 KB |
testcase_30 | AC | 181 ms
38,376 KB |
testcase_31 | AC | 175 ms
38,468 KB |
testcase_32 | AC | 173 ms
38,440 KB |
testcase_33 | AC | 179 ms
38,372 KB |
testcase_34 | AC | 178 ms
38,404 KB |
testcase_35 | AC | 181 ms
38,376 KB |
testcase_36 | AC | 176 ms
38,416 KB |
testcase_37 | AC | 175 ms
38,400 KB |
testcase_38 | AC | 176 ms
38,352 KB |
testcase_39 | AC | 175 ms
38,380 KB |
testcase_40 | AC | 183 ms
38,416 KB |
testcase_41 | AC | 172 ms
38,376 KB |
testcase_42 | AC | 171 ms
38,356 KB |
testcase_43 | AC | 174 ms
38,468 KB |
testcase_44 | AC | 177 ms
38,376 KB |
コンパイルメッセージ
main.cpp: In function 'std::ostream& operator<<(std::ostream&, const BinaryIndexedTree<T>&)': main.cpp:56:5: warning: no return statement in function returning non-void [-Wreturn-type] 56 | } | ^
ソースコード
/** * @FileName a.cpp * @Author kanpurin * @Created 2020.10.02 03:32:58 **/ #include "bits/stdc++.h" using namespace std; typedef long long ll; template< typename T > struct BinaryIndexedTree { std::vector< T > data; BinaryIndexedTree() {} BinaryIndexedTree(int sz) { data.assign(++sz, 0); } inline T sum(int k) const { T ret = 0; for (k; k > 0; k -= k & -k) ret += data[k]; return (ret); } inline T sum(int left, int right) const { return sum(right) - sum(left); } inline void add(int k, T x) { for (++k; k < data.size(); k += k & -k) data[k] += x; } int lower_bound(ll k) const { if (k <= 0) return 0; int res = 0; int N = 1; while (N < (int)data.size()) N *= 2; for (int i = N / 2; i > 0; i /= 2) { if (res + i < (int)data.size() && data[res + i] < k) { k -= data[res + i]; res += i; } } return res; } friend std::ostream& operator<<(std::ostream &os, const BinaryIndexedTree &bit) { os << "[ "; for (int i = 0; i < bit.data.size() - 1; i++) { os << bit.sum(i, i + 1); if (i < bit.data.size() - 2) os << ", "; } os << " ]"; } }; template< int MOD > struct mint { public: long long x; mint(long long x = 0) :x((x%MOD+MOD)%MOD) {} mint(std::string &s) { long long z = 0; for (int i = 0; i < s.size(); i++) { z *= 10; z += s[i] - '0'; z %= MOD; } this->x = z; } mint& operator+=(const mint &a) { if ((x += a.x) >= MOD) x -= MOD; return *this; } mint& operator-=(const mint &a) { if ((x += MOD - a.x) >= MOD) x -= MOD; return *this; } mint& operator*=(const mint &a) { (x *= a.x) %= MOD; return *this; } mint& operator/=(const mint &a) { long long n = MOD - 2; mint u = 1, b = a; while (n > 0) { if (n & 1) { u *= b; } b *= b; n >>= 1; } return *this *= u; } mint operator+(const mint &a) const { mint res(*this); return res += a; } mint operator-() const {return mint() -= *this; } mint operator-(const mint &a) const { mint res(*this); return res -= a; } mint operator*(const mint &a) const { mint res(*this); return res *= a; } mint operator/(const mint &a) const { mint res(*this); return res /= a; } friend std::ostream& operator<<(std::ostream &os, const mint &n) { return os << n.x; } friend std::istream &operator>>(std::istream &is, mint &n) { long long x; is >> x; n = mint(x); return is; } bool operator==(const mint &a) const { return this->x == a.x; } mint pow(long long k) const { mint ret = 1; mint p = this->x; while (k > 0) { if (k & 1) { ret *= p; } p *= p; k >>= 1; } return ret; } }; constexpr int MOD = 1e9 + 7; int main() { constexpr int INF = 1e9 + 6; int n;cin >> n; vector<pair<int,int>> a(n); vector<pair<int,int>> lis(n); int max_lis = 0; vector<int> count_lis(n+1,0); { vector<int> v(n,INF); for (int i = 0; i < n; i++) { cin >> a[i].first; a[i].second = i; auto lb = lower_bound(v.begin(), v.end(), a[i].first); *lb = a[i].first; lis[i].first = lb - v.begin() + 1; lis[i].second = count_lis[lis[i].first]++; max_lis = max(max_lis,lis[i].first); } } vector<BinaryIndexedTree<mint<MOD>>> bit(n+1); bit[0] = BinaryIndexedTree<mint<MOD>>(1); bit[0].add(0,1); for (int i = 0; i < n; i++) { bit[i+1] = BinaryIndexedTree<mint<MOD>>(count_lis[i+1]); } sort(a.begin(), a.end(), [](pair<int,int> a,pair<int,int> b){ return a.first != b.first ? a.first < b.first : a.second > b.second; }); mint<MOD> ans = 0; vector<set<pair<int,int>>> st(n+1); st[0].insert({-1,0}); for (int i = 0; i < n; i++) { int idx = a[i].second; auto k = st[lis[idx].first-1].lower_bound({idx,0}); k--; mint<MOD> t = bit[lis[idx].first-1].sum(k->second+1); bit[lis[idx].first].add(lis[idx].second,t); st[lis[idx].first].insert({idx,lis[idx].second}); } cout << bit[max_lis].sum(count_lis[max_lis]) << endl; return 0; }