fjut3467(gcd+dp)

题目链接

求一个序列最长不互质子序列

第一层循环遍历序列，第二层遍历该数的因子，转移方程为f[i] = max(f[i], pre[prime[j]])

f[i]表示到i位时的最长子序列长度，prime[j]是a[i]的质因子，pre[prime[j]]表示上一个公约数含有prime[i]的数的最长子序列长度

这里需要提前打表打出每个数的因子，我一开始只能想到On根号n的打表，其实这里可以利用素数打表的过程先找出质数，再去找出含有这个质数的约数。而且int范围内的数的质因子不会超过10个因为2*3*5*7*11*...29就超过int了。

 

#include <cstdio>
#include <cstring>
#include <algorithm>

using namespace std;
const int mx = 2e5 + 5;
int prime[mx][10], vis[mx], a[mx], pre[mx], f[mx];

void init() {
	for (int i = 2; i <= 100000; i++) {
		if (!vis[i]) {
			prime[i][++prime[i][0]] = i;
			for (int j = 2; i*j <= 100000; j++) {
				prime[i*j][++prime[i*j][0]] = i;
				vis[i*j] = 1;
			}
		}
	}
}

int main() {
	init();
	memset(f, 0, sizeof(f));
	memset(pre, 0, sizeof(pre));
	
	int n;
	scanf("%d", &n);
	//for (int i = 1; i <= n; i++) f[i] = 1;
	for (int i = 1; i <= n; i++) {
		scanf("%d", &a[i]);
		for (int j = 1; j <= prime[a[i]][0]; j++) {
			f[i] = max(f[i], pre[prime[a[i]][j]]+1);
		}
		for (int j = 1; j <= prime[a[i]][0]; j++) {
			pre[prime[a[i]][j]] = f[i];
		}
	}
	int ans = 0;
	for (int i = 1; i <= n; i++)
		ans = max(ans, f[i]);
	printf("%d\n", ans>0?ans:1);
 	return 0;
}