結果

問題 No.979 Longest Divisor Sequence
ユーザー i`d like to eat papico
提出日時 2024-02-09 15:38:08
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 1,398 ms / 2,000 ms
コード長 1,599 bytes
コンパイル時間 2,597 ms
コンパイル使用メモリ 79,648 KB
最終ジャッジ日時 2025-02-19 02:58:21
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 16
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <vector>
#include <algorithm>
using namespace std;

#define MAXN 300009
#define MAXA 300009

vector<int> calc_divisors(int N) {
    // 答えを表す集合
    vector<int> res;

    // 各整数 i が N の約数かどうかを調べる
    for (long long i = 1; i * i <= N; ++i) {
        // i が N の約数でない場合はスキップ
        if (N % i != 0) continue;

        // i は約数である
        res.push_back(i);

        // N ÷ i も約数である (重複に注意)
        if (N / i != i) res.push_back(N / i);
    }

    // 約数を小さい順に並び替えて出力
    sort(res.begin(), res.end());
    return res;
}

int main(void) {
	int N;
	cin >> N;
	//step l_memoに最大長を格納 初期値は0
	vector<int> DP(N, 0), A(N), l_memo(MAXA,0);
	for (int i = 0; i < N; i++)  cin >> A[i];
	//step2 DPを埋める
	for (int i = 0; i < N; i++) {
        vector<int> divisors = calc_divisors(A[i]);
        //divisorsの中で最大長を求める
        int max_len = 0;
        for (int j = 0; j < divisors.size() - 1; j++) {
            if (max_len < l_memo[divisors[j]]) max_len = l_memo[divisors[j]];
        }
        DP[i] = max_len + 1;
        //既にもっと長いのが入っている可能性がある(実際には無いと思われる)
        if (l_memo[A[i]] < DP[i]) {
            l_memo[A[i]] = DP[i];
        }
	}
    //step3 LDSを見つける(O(N))
    int ans = 0;
    for (int i = 0; i < N; i++) {
        if (ans < DP[i]) ans = DP[i];
        //cout << DP[i] << " ";
    }
    cout <<ans << endl;
    return 0;
}
0