21天好习惯第一期－11

D、小L的数列

D

思路：
引理：对于质数 pp，若存在 x<y<zx<y<z 且 x,y,zx,y,z 都是 pp 的倍数，那么选出的序列一定不会让 x,zx,z相邻。

先筛出1e5以内的所有素数，然后对序列排序
根据定理，对每个数zz，分解质因数(直接枚举质数来加速)pp，然后dp[z]=max(dp[z],dp[y]+1)dp[z]=max(dp[z],dp[y]+1)，yy与zz有相同的质因数pp且最靠近它(y<zy<z)

#include<bits/stdc++.h>
using namespace std;
const int maxn=1e5+7;
typedef long long ll;
bool vis[maxn];
vector<int>prime;
inline void initial() {
    for (int i = 2; i < maxn ; ++i) {
        if (!vis[i])
            prime.emplace_back(i);
        for (int j = 0,size=prime.size(); j < size; ++j) {
            if (i * prime[j] >= maxn) break;
            vis[i * prime[j]] = true;
            if (i % prime[j] == 0) break;
        }
    }
}
int  dp[maxn] , last[maxn];
bool a[maxn];
inline int max(int a,int b) {
    return a>b?a:b;
}
inline int read() {
    char c=getchar();
    int x=0;
    bool f=0;
    for(; !isdigit(c); c=getchar())f^=!(c^45);
    for(; isdigit(c); c=getchar())x=(x<<1)+(x<<3)+(c^48);
    if(f)x=-x;
    return x;
}
int main() {
    initial();
    int T = read() , n , k;
    while(T--) {
        n = read();
        memset(last , 0 , sizeof(last));
        memset(a , 0 , sizeof(a));
        for(int i = 1 ; i <= n ; ++i) {
            k = read();
            a[k] = 1;
        }
        int m (0),res=1;
        dp[1] = 1;
        for(int i = 2 ,y; i <= 100000 ; ++i) {
            if(a[i] == 0)   continue;
            dp[y=i]=1;
            for(auto x : prime) {
                if(x > y)   break;
                if(!vis[y]) {
                    if(last[y])
                        dp[i] = max(dp[i] , dp[last[y]] + 1);
                    last[y] = i;
                    break;
                }
                if(y % x == 0) {
                    if(last[x])
                        dp[i] = max(dp[i] , dp[last[x]] + 1);
                    last[x] = i;
                    while(y % x == 0)   y /= x;
                }
            }
            res = max(res , dp[i]);
        }
        printf("%d\n" , res);
    }
    return 0;
}