Codeforces 265D Good Sequences 贪心+DP

最新推荐文章于 2021-02-24 22:29:56 发布

原创最新推荐文章于 2021-02-24 22:29:56 发布 · 497 阅读

1 ·

本内容遵循CC 4.0 BY-SA版权协议

收录于

CODEFORCES_ 同时被 2 个专栏收录

232 篇文章

订阅专栏

Dynamic Programming

143 篇文章

订阅专栏

本文介绍了一种O(nlogn)复杂度的算法，用于解决从n个合法数中找出最长序列的问题，要求序列中任意两个连续数满足递增且最大公约数大于1。通过预处理素因子和动态规划实现。

点击打开链接

题意: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;
}