結果

問題 No.979 Longest Divisor Sequence
ユーザー ok
提出日時 2020-01-31 22:38:37
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
TLE  
(最新)
AC  
(最初)
実行時間 -
コード長 2,649 bytes
コンパイル時間 1,029 ms
コンパイル使用メモリ 97,296 KB
実行使用メモリ 144,812 KB
最終ジャッジ日時 2024-09-17 09:27:18
合計ジャッジ時間 5,100 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 15 TLE * 1
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<iostream>
#include<string>
#include<iomanip>
#include<cmath>
#include<vector>
#include<algorithm>
#include <functional>
#include <limits>

using namespace std;

#define int long long
#define endl "\n"

constexpr long long INF = (long long)1e18;
constexpr long long MOD = 1'000'000'007; 

string yn(bool f){return f?"Yes":"No";}
string YN(bool f){return f?"YES":"NO";}


template<typename T>
class segment_tree{
	int n;
	T fval;
	function<T(T,T)> operation;
	vector<T> dat;
	
	T _query(int a, int b, int k, int l, int r){
		if(r <= a || b <= l){
			return fval;
		}
		
		if(a <= l && r <= b){
			return dat[k];
		}else {
			T vl = _query(a, b, k * 2 + 1, l, (l + r) / 2 );
			T vr = _query(a, b, k * 2 + 2, (l + r) / 2, r);
			return operation(vl, vr);
		}
	}
	
public:
	
	segment_tree(){
		
	}
	
	segment_tree(int n_, T val, function<T(T,T)> fun){
		init(n_, val, fun);
	}
	
	~segment_tree(){
		
	}
	
	void init(int n_, T val, function<T(T,T)> fun){
		fval = val;
		
		operation = fun;
		
		n = 1;
		
		while(n < n_) n *= 2;
		
		dat.resize(2 * n - 1);
	
		for(int i = 0; i < 2 * n - 1; i++) dat[i] = fval;
	}

	void update(int k, T a){
		k += n - 1;
		dat[k] = a;
		while(k > 0){
			k = (k - 1) / 2;
			dat[k] = operation(dat[k*2 + 1],dat[k*2 + 2]);
		}
	}
	
	T query(int a, int b){
		return _query(a, b, 0, 0, n);
	}
	
};


signed main(){
	cin.tie(nullptr);
	ios::sync_with_stdio(false);
	cout<<fixed<<setprecision(10);
	
	int N, maxA = 0;
	int ans = 0;
	vector<int> A;
	vector<segment_tree<int>> seg;
	
	cin>>N;
	
	A.resize(N);
	
	for(int i = 0; i < N; i++){
		cin>>A[i];
		maxA = max(maxA, A[i]);
	}	
	
	seg.resize(maxA+1);
	
	for(int i = 1; i <= maxA; i++){
		// cout<<i<<endl;
		seg[i].init((maxA+i)/i+4,0, [](int first, int second){ return max(first,second);});
	}
	
	for(int i = N-1; i >= 0; i--){
		int maximum = seg[A[i]].query(1, (maxA+A[i])/A[i] + 1);
		// int maximum = 0;
		// for(int j = 0; j < (maxA+A[i])/A[i]; j++){
			// cout<<j<<" "<<seg[A[i]].query(j, j+ 1)<<endl;
		// }
		// cout<<i<<" max = "<<maximum<<endl;
		
		for(int j = 1; j * j <= A[i]; j++){
			if(A[i]%j == 0){
				{
					// int hoge = seg[j].query(A[j]/j-1, A[j]/j);
					
					// hoge = max(hoge, maximum+1);
					
					// seg[j].update(A[i]/j-1, hoge);
					seg[j].update(A[i]/j-1, maximum+1);
				}
				
				{
					int z = A[i]/j;
					// int hoge = seg[z].query(A[j]/z-1, A[j]/z);
					
					// hoge = max(hoge, maximum+1);
					
					// seg[z].update(A[i]/z-1, hoge);
					seg[z].update(A[i]/z-1, maximum+1);
				}
			}
		}
	}
	 
	for(int i = 1; i <= maxA; i++){
		ans = max(ans, seg[i].query(0, (maxA+i)/i + 2));
	}
	
	cout<<ans<<endl;
	
	return 0;
}
0