来啦来啦,寒假复健第一题cf1475g
一道div3的g题,dp,作为一个忙期末1个月没做题的废人,先整点最爱的Dp康复一下。
题目不是很难,给定一个序列,问最少剔除几个数后能使序列中的每个数之间都呈倍数关系。
很明显,倍数关系是传递的(其实太久没做题想了好一会儿),所以从小到大排序后,就变成找一个最长的倍数关系的链。
那么我们开一个dp数组,dp[i]表示以i为结尾的数最长的长度是多少(有点像LIS的思想)。那么对于每个数只要分解出因子就可以了,所以时间复杂度是O(nsqrt(n))的,n=20W完全可行。
然后在状态转移的时候小心一点要先把大的因子转了(主要是防止1与本身的时候重复更新)。
下附代码:
1 #include<bits/stdc++.h> 2 using namespace std; 3 const int INF=0X3f3f3f3f; 4 int a[200005],n; 5 int dp[200005]; 6 int main(){ 7 int T; 8 scanf("%d",&T); 9 while (T--){ 10 scanf("%d",&n); 11 memset(dp,-INF,sizeof(dp)); 12 for (int i=1; i<=n; i++){ 13 scanf("%d",&a[i]); 14 } 15 sort(a+1,a+1+n); 16 dp[1]=0; 17 for (int i=1; i<=n; i++){ 18 for (int j=1; j*j<=a[i]; j++){ 19 if (a[i]%j==0){ 20 dp[a[i]]=max(dp[a[i]],dp[a[i]/j]+1); 21 if (a[i]/j!=j) 22 dp[a[i]]=max(dp[a[i]],dp[j]+1); 23 } 24 } 25 } 26 int res=0; 27 for (int i=1; i<=200000; i++){ 28 res=max(res,dp[i]); 29 } 30 res=n-res; 31 printf("%d\n",res); 32 } 33 }