結果

問題 No.992 最長増加部分列の数え上げ
ユーザー chocorusk
提出日時 2020-02-15 04:56:21
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 200 ms / 2,000 ms
コード長 2,144 bytes
コンパイル時間 1,358 ms
コンパイル使用メモリ 123,984 KB
実行使用メモリ 16,256 KB
最終ジャッジ日時 2024-10-06 13:41:36
合計ジャッジ時間 8,657 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 42
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <cstdio>
#include <cstring>
#include <iostream>
#include <string>
#include <cmath>
#include <bitset>
#include <vector>
#include <map>
#include <set>
#include <queue>
#include <deque>
#include <algorithm>
#include <complex>
#include <unordered_map>
#include <unordered_set>
#include <random>
#include <cassert>
#include <fstream>
#include <utility>
#include <functional>
#include <time.h>
#include <stack>
#include <array>
#define popcount __builtin_popcount
using namespace std;
typedef long long int ll;
typedef pair<int, ll> P;
const ll MOD=1e9+7;
template<typename Monoid>
struct SegmentTree{
	using F=function<Monoid(Monoid, Monoid)>;
	int sz;
	vector<Monoid> seg;
	const F f;
	const Monoid e;

	SegmentTree(int n, const F f, const Monoid &e): f(f), e(e){
		sz=1;
		while(sz<n) sz<<=1;
		seg.resize(2*sz, e);
	}

	SegmentTree(int n, const F f, const Monoid &e, vector<Monoid> v): f(f), e(e){
		sz=1;
		while(sz<n) sz<<=1;
		seg.resize(2*sz, e);
		for(int i=0; i<n; i++) seg[i+sz]=v[i];
		for(int i=sz-1; i>=1; i--){
			seg[i]=f(seg[2*i], seg[2*i+1]);
		}
	}

	void update(int k, const Monoid &x){
		k+=sz;
		seg[k]=x;
		while(k>1){
			k>>=1;
			seg[k]=f(seg[2*k], seg[2*k+1]);
		}
	}

	Monoid query(int a, int b){
		a+=sz, b+=sz;
		Monoid ret=e;
		for(;a<b; a>>=1, b>>=1){
			if(b&1) ret=f(ret, seg[--b]);
			if(a&1) ret=f(ret, seg[a++]);
		}
		return ret;
	}

	Monoid operator[](const int &k) const{
		return seg[k+sz];
	}
};
int main()
{
	int n; cin>>n;
	int a[200020];
	vector<int> va(n);
	for(int i=0; i<n; i++){
		cin>>a[i];
		va[i]=a[i];
	}
	sort(va.begin(), va.end());
	va.erase(unique(va.begin(), va.end()), va.end());
	int m=va.size()+1;
	for(int i=0; i<n; i++) a[i]=lower_bound(va.begin(), va.end(), a[i])-va.begin()+1;
	P dp[200020];
	dp[0]=P(0, 1);
	auto f=[](P a, P b){
		if(a.first!=b.first) return max(a, b);
		else return P(a.first, (a.second+b.second)%MOD);
	};
	SegmentTree<P> seg(m, f, P(-1, 0));
	seg.update(0, dp[0]);
	for(int i=0; i<n; i++){
		P p=seg.query(0, a[i]);
		dp[i+1]=P(p.first+1, p.second);
		seg.update(a[i], f(seg[a[i]], dp[i+1]));
	}
	P ans=seg.query(0, m);
	cout<<ans.second<<endl;
	return 0;
}
0