結果

問題 No.992 最長増加部分列の数え上げ
ユーザー leaf_1415leaf_1415
提出日時 2020-02-14 22:34:09
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 176 ms / 2,000 ms
コード長 2,609 bytes
コンパイル時間 1,119 ms
コンパイル使用メモリ 76,400 KB
実行使用メモリ 28,760 KB
最終ジャッジ日時 2024-10-06 12:58:15
合計ジャッジ時間 7,479 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 9 ms
19,524 KB
testcase_01 AC 7 ms
20,224 KB
testcase_02 AC 9 ms
20,360 KB
testcase_03 AC 8 ms
19,636 KB
testcase_04 AC 64 ms
22,620 KB
testcase_05 AC 47 ms
22,144 KB
testcase_06 AC 78 ms
23,208 KB
testcase_07 AC 57 ms
23,268 KB
testcase_08 AC 33 ms
21,340 KB
testcase_09 AC 59 ms
22,292 KB
testcase_10 AC 78 ms
22,964 KB
testcase_11 AC 98 ms
24,292 KB
testcase_12 AC 25 ms
21,056 KB
testcase_13 AC 56 ms
22,164 KB
testcase_14 AC 55 ms
23,220 KB
testcase_15 AC 25 ms
20,432 KB
testcase_16 AC 151 ms
27,128 KB
testcase_17 AC 33 ms
21,028 KB
testcase_18 AC 54 ms
22,240 KB
testcase_19 AC 95 ms
23,904 KB
testcase_20 AC 172 ms
27,856 KB
testcase_21 AC 168 ms
28,012 KB
testcase_22 AC 174 ms
27,988 KB
testcase_23 AC 174 ms
27,864 KB
testcase_24 AC 171 ms
27,888 KB
testcase_25 AC 176 ms
27,860 KB
testcase_26 AC 176 ms
27,892 KB
testcase_27 AC 175 ms
27,732 KB
testcase_28 AC 171 ms
27,820 KB
testcase_29 AC 168 ms
27,876 KB
testcase_30 AC 110 ms
27,856 KB
testcase_31 AC 111 ms
27,864 KB
testcase_32 AC 112 ms
27,856 KB
testcase_33 AC 110 ms
27,860 KB
testcase_34 AC 111 ms
27,732 KB
testcase_35 AC 97 ms
27,476 KB
testcase_36 AC 98 ms
27,628 KB
testcase_37 AC 99 ms
27,504 KB
testcase_38 AC 97 ms
27,468 KB
testcase_39 AC 99 ms
27,456 KB
testcase_40 AC 118 ms
28,296 KB
testcase_41 AC 122 ms
28,268 KB
testcase_42 AC 121 ms
28,760 KB
testcase_43 AC 121 ms
28,240 KB
testcase_44 AC 121 ms
28,372 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <map>
#include <vector>
#include <algorithm>
#define llint long long
#define inf 1e18
#define mod 1000000007

using namespace std;

struct SegTree{
	int size;
	vector<llint> seg;

	SegTree(){}
	SegTree(int size){
		this->size = size;
		seg.resize(1<<(size+1));
	}

	void init()
	{
		for(int i = 0; i < (1<<(size+1)); i++) seg[i] = -inf;
	}

	void update(int i, llint val)
	{
		i += (1 << size);
		seg[i] = val;
		while(i > 1){
			i /= 2;
			seg[i] = max(seg[i*2], seg[i*2+1]);
		}
	}

	llint query(int a, int b, int k, int l, int r)
	{
		if(b < l || r < a) return -inf;
		if(a <= l && r <= b) return seg[k];
		llint lval = query(a, b, k*2, l, (l+r)/2);
		llint rval = query(a, b, k*2+1, (l+r)/2+1, r);
		return max(lval, rval);
	}
	llint query(int a, int b)
	{
		return query(a, b, 1, 0, (1<<size)-1);
	}
};

llint n;
llint a[200005];
vector<llint> comp;
SegTree seg(18);
llint dp[200005], dp2[200005];
vector<llint> vec[200005], svec[200005];

int main(void)
{
  ios::sync_with_stdio(0);
  cin.tie(0);

  cin >> n;
  for(int i = 1; i <= n; i++) cin >> a[i];
  for(int i = 1; i <= n; i++) comp.push_back(a[i]);
  comp.push_back(-inf);
  sort(comp.begin(), comp.end());
  comp.erase(unique(comp.begin(), comp.end()), comp.end());
  for(int i = 1; i <= n; i++) a[i] = lower_bound(comp.begin(), comp.end(), a[i]) - comp.begin();

  seg.init();
  seg.update(0, 0);
  dp[0] = 0;
  for(int i = 1; i <= n; i++){
    dp[i] = seg.query(0, a[i]-1) + 1;
    seg.update(a[i], dp[i]);
    //for(int j = 0; j < 10; j++) cout << seg.query(j, j) << " "; cout << endl;
  }
  llint lis = seg.query(0, (int)comp.size());

  //for(int i = 0; i <= n; i++) cout << dp[i] << " "; cout << endl;

  for(int i = 0; i <= n; i++) vec[dp[i]].push_back(i);
  svec[0].push_back(1);
  dp2[0] = 1;

  for(int i = 1; i <= n; i++){
    vector<llint> &cvec = vec[dp[i]-1];
    llint p = lower_bound(cvec.begin(), cvec.end(), i) - cvec.begin();

    //cout << p << endl;

    llint ub = p, lb = -1, mid;
    while(ub-lb>1){
      mid = (ub+lb)/2;
      if(a[cvec[mid]] < a[i]) ub = mid;
      else lb = mid;
    }
    //cout << ub << endl;
    vector<llint> &csvec = svec[dp[i]-1];
    dp2[i] = csvec[p-1];
    if(ub > 0) dp2[i] += mod - csvec[ub-1] + mod, dp2[i] %= mod;
    if(svec[dp[i]].size() == 0) svec[dp[i]].push_back(dp2[i]);
    else svec[dp[i]].push_back((svec[dp[i]].back() + dp2[i]) % mod);
  }

  //for(int i = 0; i <= n; i++) cout << dp2[i] << " "; cout << endl;

  llint ans = 0;
  for(int i = 0; i <= n; i++){
    if(dp[i] == lis) ans += dp2[i], ans %= mod;
  }
  cout << ans << endl;

  return 0;
}
0