原題鏈接
思路：
看到這題第一眼覺得要用些數論的知識優化，決定先寫個暴力，
假設dp[i]表示列舉到第i項的時候最長的好的序列長度，轉移仿照LIS的n^2轉移，過了80，
然后就不會優化了，
宇巨說可以列舉因子進行優化：
///dp[i]表示以數i為結尾的最長的好的序列長度
///dp1[i]表示因數包括i的數結尾的最長好的序列長度
這樣轉移的時候將第二層for改成sqrt(n)的復雜度，并且在程序中用dp1更新dp，列舉完因子后，再用此數的dp更新因子的dp1，
ps:
1.要對原陣列進行排序，因為第一個要求
2.注意因子為1時不進行轉移
事實證明,dp的狀態選擇多么重要（逃
代碼:
#pragma GCC optimize(2)
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<ll,ll>PLL;
typedef pair<int,int>PII;
typedef pair<double,double>PDD;
#define I_int ll
inline ll read()
{
ll x=0,f=1;
char ch=getchar();
while(ch<'0'||ch>'9')
{
if(ch=='-')f=-1;
ch=getchar();
}
while(ch>='0'&&ch<='9')
{
x=x*10+ch-'0';
ch=getchar();
}
return x*f;
}
#define read read()
#define closeSync ios::sync_with_stdio(0);cin.tie(0);cout.tie(0)
#define multiCase int T;cin>>T;for(int t=1;t<=T;t++)
#define rep(i,a,b) for(int i=(a);i<=(b);i++)
#define repp(i,a,b) for(int i=(a);i<(b);i++)
#define per(i,a,b) for(int i=(a);i>=(b);i--)
#define perr(i,a,b) for(int i=(a);i>(b);i--)
ll ksm(ll a,ll b,ll p)
{
ll res=1;
while(b)
{
if(b&1)res=res*a%p;
a=a*a%p;
b>>=1;
}
return res;
}
#define PI acos(-1)
#define x first
#define y second
const int maxn=1e6+7,inf=0x3f3f3f3f;
int n,a[maxn],dp[maxn],dp1[maxn],w[maxn];
int gcd(int a,int b){
return b==0?a:gcd(b,a%b);
}
int main()
{
int t=read;
while(t--){
n=read;
rep(i,1,n) a[i]=read;
sort(a+1,a+1+n);
///dp[i]表示以數i為結尾的最長的好的序列長度
///dp1[i]表示因數包括i的數結尾的最長好的序列長度
memset(dp,0,sizeof dp);
memset(dp1,0,sizeof dp1);
int res=1;
for(int i=1;i<=n;i++){
dp[a[i]]=1;
int idx=0;
if(a[i]==1) continue;
for(int j=1;j*j<=a[i];j++){
if(a[i]%j==0){
int x1=a[i]/j,x2=j;
if(x1==x2){
if(x1!=1) dp[a[i]]=max(dp[a[i]],dp1[x1]+1),w[++idx]=x1;
}
else{
if(x1!=1) dp[a[i]]=max(dp[a[i]],dp1[x1]+1),w[++idx]=x1;
if(x2!=1) dp[a[i]]=max(dp[a[i]],dp1[x2]+1),w[++idx]=x2;
}
}
}
for(int j=1;j<=idx;j++) dp1[w[j]]=max(dp1[w[j]],dp[a[i]]);///列舉a[i]的因子
}
for(int i=1;i<=n;i++) res=max(res,dp[a[i]]);
printf("%d\n",res);
}
return 0;
}
轉載請註明出處，本文鏈接：https://www.uj5u.com/yidong/267157.html
標籤：其他