【题解 && 海亮集训 && 动态规划】 先进序列（sequence）

///这个月似乎好久没写博客了。。赶紧补一篇

---

题目描述：

萌猪大统领温温正在为他看到的所有事物评选先进。
温温已经评选出了 n 个先进整数 ai，现在，他决定在此基础上评选一下先进序列。
一个整数序列 x1,x2,…xk 是先进的，当且仅当它满足以下三个条件
1、这个序列是严格递增的。
2、序列中，任意两个相邻的元素都不是互质的，即 gcd(xi, xi+1) > 1 对 1 ≤ i < k 成立。
3、序列中的所有元素都是先进整数。
现在温温想知道，最长的先进序列的长度是多少。

---

Solution:

暴力做法大家都会，这里不在多述。

这题有一个坑点，那就是题目中没有说给出的序列是有序的，但是给出的样例却是有序的（坑啊。。。）
题目中gcd这一个条件有些麻烦，直接暴力求解肯定会超时，我们需要想优化。

• 设$dp[i]$ 表示已i结尾的最长上升序列长度
• 设 $d[i]$ 表示所有拥有质因数$i$ 的数所对应的位置的$dp$ 的最大值。

两者间可以“优势互补”，由于题目中要求相邻的数不能互质，所以dp便可以有d更新而来。

枚举i，对于当前数a[i]，枚举它的质因数，可有$max(d[p{]}_{p\in 质因数a[i]})$+1直接转移到dp，既满足最大值，又满足相邻之间不互质
同理，d也可以由dp转移而来，具体细节参考代码。

---

Code

#include<bits/stdc++.h>
using namespace std;
int n,cnt=0;
int a[1001000];
int dp[1001000];
int d[1001000];
int pri[1000100];
bool isp[1001000];
vector < int > Pr[101000];

int main(){
	freopen("sequence.in","r",stdin);
	freopen("sequence.out","w",stdout);
	scanf("%d",&n);
	for (int i=1;i<=n;i++)
	  scanf("%d",&a[i]),dp[i]=1;
	sort(a+1,a+n+1);
	memset(isp,1,sizeof(isp));
	isp[1]=0;
	for (int i=2;i<=200000;i++)
	  if (isp[i])
	    for (int j=2;j<=200000/i;j++)
	      isp[i*j]=0;
	for (int i=1;i<=200000;i++)
	  if (isp[i]) pri[++cnt]=i;
	for (int i=1;i<=n;i++){
	    int x=a[i];
	    for (int j=1;j<=cnt && pri[j]<=x && x!=1;j++)
	      if (x%pri[j] == 0){
		      while (x%pri[j]==0) x/=pri[j];
		      Pr[i].push_back(pri[j]);
		  }
		if (x!=1) Pr[i].push_back(x);
	}
//	for (int i=1;i<=n;i++){
//	    for (int j=0;j<Pr[i].size();j++)
//	      cout<<Pr[i][j]<<' ';
//	    cout<<endl;
//	}
	for (int i=0;i<Pr[1].size();i++)
	  d[Pr[1][i]] = dp[1];
//	for (int i=1;i<=10;i++)cout<<d[i]<<' ';
//	cout<<endl;
	for (int i=2;i<=n;i++){
	    for (int j=0;j<Pr[i].size();j++)
	      dp[i] = max(dp[i],d[Pr[i][j]]+1);
	    for (int j=0;j<Pr[i].size();j++)
	      d[Pr[i][j]] = max(dp[i],d[Pr[i][j]]);
	}
//	for (int i=1;i<=n;i++) cout<<dp[i]<<' ';
//	cout<<endl;
	int maxx=0;
	for (int i=1;i<=n;i++) maxx=max(maxx,dp[i]);
	printf("%d",maxx);
    return 0;
}