Codeforces 265D Good Sequences 贪心+DP

点击打开链接

题意:n个合法的数a[i],从其中选出k个数满足,a[i]<a[i+1], gcd(a[i],a[i+1])>1 问k最大为?
n,a[i]<=1e5


O(nlogn)晒出a[i]的素因子,设d[x]:以x结尾的最长合法序列,根据di=gcd(a[i],a[i-1])转移状态
a[i-1]:为a[i]某个素因子的倍数 a[i-1]=k*v[i]贪心取最大的一个显然最优.因为d[kx]>=d[x](kx存在时)
用f[x]记录当前x的最大倍数. d[x]=max(d[f[v[j]]+1,d[x]). 

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll inf=2e15;
const int M=1e5;
const int N=2e5+20;
int n,a[N],d[N],vis[N],f[N];
vector<int> v[N];
void init()
{
	memset(vis,0,sizeof(vis));
	for(int i=2;i<=M;i++)
	{
		if(vis[i]==0)
		{
			for(int j=i+i;j<=M;j+=i)
				v[j].push_back(i),vis[j]=1;
		}
	}
}
int main()
{
	init();
	while(scanf("%d",&n)!=EOF)
	{
		for(int i=1;i<=n;i++)
			scanf("%d",&a[i]);	
		memset(d,0,sizeof(d));
		memset(f,0,sizeof(f));
		//O(nd)
		int ans=1;
		for(int i=1;i<=n;i++)
		{
			d[a[i]]=1;
			for(int j=0;j<v[a[i]].size();j++)
			{
				int x=v[a[i]][j];
				d[a[i]]=max(d[a[i]],d[f[x]]+1);	
				f[x]=a[i];
				//for(int k=x;k<a[i];k+=x)
				//	d[a[i]]=max(d[a[i]],d[k]+1);
				//cout<<a[i]<<' '<<x<<' '<<f[x]<<endl;
			}
			f[a[i]]=a[i];
			ans=max(ans,d[a[i]]);
		}
		cout<<ans<<endl;
	}	
	return 0;
}