結果
問題 |
No.2218 Multiple LIS
|
ユーザー |
![]() |
提出日時 | 2025-09-16 19:27:54 |
言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 114 ms / 3,000 ms |
コード長 | 2,301 bytes |
コンパイル時間 | 3,596 ms |
コンパイル使用メモリ | 191,992 KB |
実行使用メモリ | 15,236 KB |
最終ジャッジ日時 | 2025-09-16 19:28:02 |
合計ジャッジ時間 | 7,672 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 39 |
ソースコード
#include <iostream> #include <vector> #include <algorithm> #include <cmath> #include <atcoder/all> using namespace std; using namespace atcoder; using ll = long long; using mint = modint998244353; using vi = vector<int>; using vvi = vector<vi>; using vvvi = vector<vvi>; using vll = vector<ll>; using vvll = vector<vll>; using vvvll = vector<vvll>; using vmi = vector<mint>; using vvmi = vector<vmi>; using vvvmi = vector<vvmi>; #define all(a) (a).begin(), (a).end() #define rep2(i, m, n) for (int i = (m); i < (n); ++i) #define rep(i, n) rep2(i, 0, n) #define drep2(i, m, n) for (int i = (m)-1; i >= (n); --i) #define drep(i, n) drep2(i, n, 0) class Primes{ private: int size_p; vector<bool> is_p; vector<int> primes; vector<int> mobius; void eratos(int n){ for(int i=2; i<=n; i++){ if(minfactor[i] < i)continue; primes.emplace_back(i); mobius[i] = -1; is_p[i] =true; for(int j = 2*i; j<=n; j+=i){ minfactor[j] = min(i, minfactor[j]); mobius[j] = ((j/i) % i == 0) ? 0 : -mobius[j]; } } } public: vector<int> minfactor; Primes(int n){ size_p = n; rep(i, size_p+1)is_p.emplace_back(false); rep(i, size_p+1)minfactor.emplace_back(i); rep(i, size_p+1)mobius.push_back(1); mobius[1] = 1; eratos(size_p); } bool isprime(int n){ return is_p[n]; } vector<int> listPrimes(){ return primes; } int getMobius(int n){return mobius[n];} }; int main(){ int n; cin >> n; int mx = 0; vi a(n); rep(i, n){ cin >> a[i]; mx = max(mx, a[i]); } Primes pr(mx+1); vvi cnt(mx+1); rep2(i, 1, mx+1){ for(int j = i; j <= mx; j += i){ cnt[j].push_back(i); } } vi dp(mx+1, 0); rep(i, n){ int tmp = 0; for(auto j : cnt[a[i]]){ tmp = max(tmp, dp[j]); } dp[a[i]] = max(dp[a[i]], tmp+1); } int ans = 0; rep(i, mx+1)ans = max(dp[i], ans); cout << ans << endl; return 0; }