【dp】取数游戏III

题目描述

小 C 刚学了辗转相除法，正不亦乐乎，这小 P 又出来捣乱，给小 C 留了个 难题。 给 N 个数，用 a1,a2…an来表示。现在小 P 让小 C 依次取数，第一个数可以 随意取。假使目前取得 aj，下一个数取ak(k>j)，则ak必须满足gcd(aj,ak)≥L。 到底要取多少个数呢？自然是越多越好！ 不用多说，这不仅是给小 C 的难题，也是给你的难题。

---

输入

第一行包含两个数N 和 L。 接下来一行，有 N 个数用空格隔开，依次是 a1,a2…an。

---

输出

仅包含一行一个数，表示按上述取法，最多可以取的数的个数。

---

样例输入

 

5 6
7 16 9 24 6

样例输出

3

---

提示

选取 3个数16、24、6。gcd(16,24)=8，gcd(24,6)=6。 

2≤L≤ai≤1 000 000； 
30% 的数据N≤1000； 
100% 的数据 N≤50 000

---

【题解】：

题意：如提示所说的，这个题目一开始没有写出来，后来是看了网上的博客来写出来的。

我觉得dp不简单，因为我根本没有想到可以枚举因数来做，通过因数的桥梁来搭建。

 

原始做法：我们可以枚举每一个位上，然后枚举所有在它之前的位置，

然后全部取Gcd看看能否满足题目所给的L  O(N^2)

 

DP做法：首先我们一定要想一下，gcd(a,b) | a  ,gcd(a,b) | b,我们开设两个数组，

一个叫dp[ i ] , 就是第 i 位 为末尾，能够取到最长的个数。

一个叫Last [ i ] ，就是   就是 因数为 i  在哪 取得的方案数最多。

然后我们枚举每一个a[ i ] 的所有因数，通过因数的桥梁来搭建关系，

他们的关系是：

dp[ i ] = max ( dp[ i ] , Last [ j ]  +1 )

就是说更新该位置上的值，这个位置上的值，决定它的是前面所有数的因子和他有关的。

 还要更新Last [ j ] 。具体可以看到代码理解一下

 

#include<bits/stdc++.h>
using namespace std;
const int N=2000005;
int a[N],Last[N],dp[N],n,L,ans=-1;
int main()
{
    scanf("%d%d",&n,&L);
    for(int i=1;i<=n;i++){
        scanf("%d",&a[i]);
    }
    for(int i=1;i<=n;i++){
        dp[i]=1;
        for(int j=1;j*j<=a[i];j++){
            if(a[i]%j) continue;
            if(j>=L)      { dp[i]=max(Last[j]+1,dp[i]);}
            if(a[i]/j>=L) { dp[i]=max(Last[a[i]/j]+1,dp[i]);}
        }
        for(int j=1;j*j<=a[i];j++){
            if(a[i]%j) continue;
            if(j>=L) { Last[j]=max(Last[j],dp[i]);}
            if(a[i]/j>=L) { Last[a[i]/j]=max(Last[a[i]/j],dp[i]);}
        }
        ans=max(ans,dp[i]);
    }
 
    printf("%d\n",ans);
    return 0;
}