『动态规划·数论优化』先进序列

题目描述

题解

先排序，很容易得到$f[i]=f[j]+1,gcd(a[i],a[j])=1.$

我们可以对每一个数分解质因数，寻出到数组内，更新最大的$f_i$值。

然后每次只需要在质因数内找最大值更新即可。

由于分解质因数，时间复杂度:$O(nlogn)$

代码如下：

#include <bits/stdc++.h>

using namespace std;
const int N = 300000;

int n, ans = 0;
int a[N], b[N], f[N], Max[N];

vector <int> c[N];

int main(void)
{
	freopen("sequence.in","r",stdin);
	freopen("sequence.out","w",stdout);
	scanf("%d", &n);
	for (int i=1;i<=n;++i) 
		scanf("%d", a+i), b[i] = a[i];
	sort(a+1,a+n+1);
	sort(b+1,b+n+1);
	for (int i=1;i<=n;++i)
	{
		for (int j=2;j*j<=a[i];++j)
			if (a[i] % j == 0) {
				c[i].push_back(j);
				while (b[i] % j == 0) b[i] /= j;
			}
		if (b[i] > 1) c[i].push_back(b[i]);
	}
	for (int i=1;i<=n;++i)
	{
		for (int j=0;j<c[i].size();++j)
		{
			int t = c[i][j];
			f[i] = max(f[i],Max[t]+1);	
		}	
		for (int j=0;j<c[i].size();++j)
		{
			int t = c[i][j];
			Max[t] = max(Max[t],f[i]);
		}
		ans = max(ans,f[i]);
	} 
	printf("%d\n", ans);
	return 0;
}