codeforces 1475 G Strange Beauty

题目的意思是给出n个数，把他们排序后，为了使得相邻的两个数之前的数是之后的数的因子，问最少需要删除多少个其中的数。

这个题可以先设一个函数，f(x),然后再慢慢的去找f(x)要表示什么东西，显然，删除最少的数等于总数n减去最大的满足条件的那个序列长度，那么f(x)就应该表示长度，然后，这个x显然不能表示成排序后数组的下标，因为发生联系的是数与数之间的倍数关系，那么x表示为这个数列中的任意一个数更容易建立联系，然后就是找有关于这个函数的表达，那么就从最初的头开始找。
先假设只有两个数x1,x2，如果序列为1,1,2,则f(1)=cnt(1)=2, f(2)=f(1)+cnt(2),
序列若为 1 ，1 ，3 ,5，那么f(1)=2, f(3)=f(1)+cnt(3) =3,f(5)=f(1)+cnt(5)=3;
那么 f(x)=f(t)+cnt(x) ,x%t==0, 那么关键是这个t怎么去选，由于f(x/2) 能包含f(x/2/3), 那么只要考虑f(x/2), 那么就是求出 所有的这个数的素因子p，然后分别考虑f(x/p), 因为设A为少一个p因子的x的其他因子的集合，那么所有A的并集就能表示出所有的x的真因子。也就是x/p(p为所有的素因子）,涵盖了所有x的真因子的倍数。
但是x/p在原序列中可能没有，为了递推下去，那就这样，设f(x)表示若以x结尾则由原序列有的数所能组成的最大的合法序列长度，这样，问题就解决了。复杂度n*根号下n.

#include<bits/stdc++.h>
#define int long long 
#define ioss ios::sync_with_stdio(0),cin.tie(0)
using namespace std;
template<class cl>void read(cl& x)
{
    x=0;int f=0;char ch;ch=getchar();
    while(!isdigit(ch)){f=f|(ch=='-'),ch=getchar();}
    while(isdigit(ch)){x=(x<<1)+(x<<3)+(ch^48);ch=getchar();}
    x=f?-x:x;
    return ;
}
template<class cl>void put(cl x)
{
    if(x<0) putchar('-'),x=-x;
    if(x>9) put(x/10);
    putchar(x%10+'0');
    return;
}
const int N=2e5+10;
int dp[N];
int a[N];
int cnt[N];
int g[N];
int prime[N];
int cntt=0;//注意cnt与cntt 
void fip()
{
   g[1]=g[0]=1;
   for(int i=2;i<=N;i++)//i 从二开始  ，从一错误  易错
   {
      if(g[i]==0)
      {
          prime[++cntt]=i;
         for(int j=i+i;j<=N;j=j+i)
         {
             g[j]=1;
         }
      }
   }
}
int mx(int x)
{
    int z=0;
    int tt;
    int xx=x;//这里要备份.
    for(int i=1;prime[i]<=x/prime[i];i++) //这里是i=1,不是i=2,易错。
    {
        if(x%prime[i]==0)
        {
            tt=xx/prime[i];
            z=max(z,dp[tt]);
            while(x%prime[i]==0)
            {
                x=x/prime[i];
            }
        }
    }
    if(x^1)
    {
        z=max(z,dp[xx/x]);
    }
    return z+cnt[xx];
}
signed main()
{     
   fip();
   //cout<<prime[2]<<endl;
    
    int t;
    read(t);
    int n;
    while(t--)
    {
        for(int i=0;i<N;i++)
        {
            cnt[i]=0;
            dp[i]=0;
        }
        read(n);
        for(int i=1;i<=n;i++)
        {
            read(a[i]);
            cnt[a[i]]++;
        }
        sort(a+1,a+1+n);
         dp[1]=cnt[1];
     
       int res=0;
       res=max(dp[1],res);
       for(int i=2;i<=a[n];i++)
       {
           dp[i]=mx(i);//是dp[i] ,不是dp[2].
           res=max(res,dp[i]);
       }
     /*  for(int i=1;i<=n;i++)
       {
           cout<<dp[a[i]]<<" ";
       }*/
     //  cout<<endl;
       res=n-res;
       put(res);
       puts("");
    }
    return 0;
}