2020ICPC·小米网络选拔赛第一场 Intelligent Warehouse（数学）

最新推荐文章于 2022-09-29 16:23:55 发布

帅气的大马猴

最新推荐文章于 2022-09-29 16:23:55 发布

阅读量185

点赞数

文章标签：数学

本文链接：https://blog.csdn.net/weixin_45187630/article/details/109289908

版权

2020ICPC·小米网络选拔赛第一场 Intelligent Warehouse（数学）
题意：给你一个数组a包含N个数，然后找一个最长的子数组，要求每个子数组里面的任意两个数能够整除。求个子数组的长度；

题解：
1 4 2 8 5 7 样例很清楚最长就是1 2 4 8 答案为4，首先我们预处理1e7的素数，然后再去标记数组a的出现次数，然后每次我们都去枚举每个素数的倍数积然后更新最大，
dp[i*prime[j]=max(dp[i],dp[i * prime[j])
然后最后加上1的数量

#pragma GCC optimize(2)
#pragma GCC optimize(3)
#include<bits/stdc++.h>
#define ll long long
#define ull unsigned long long
#define fio ios_base::sync_with_stdio(0); cin.tie(0); cout.tie(0);
#define mse(a,b) memset(a,b,sizeof a)
#define pb push_back
using namespace std;
const int mod=1e9+7;
const int maxx=1e7+10;
using namespace std;
const long double PI = 3.14159265358979323846;
#define fio ios_base::sync_with_stdio(0); cin.tie(0); cout.tie(0);
inline int read(){ int x=0,f=1; char ch=getchar(); while (!isdigit(ch)) { if (ch=='-') f=-1; ch=getchar(); } while (isdigit(ch)) { x=x*10+ch-48;  ch=getchar();  } return x*f;}
//     ll cc = ((1ll << 62) - 1 + (1ll << 62));
ll prime[maxx];
ll vis[maxx];
ll cnt=0;
void init()
{
    mse(vis,0);
    for(int i=2;i<1e7;i++)
    {
        if(!vis[i]) prime[++cnt]=i;
        for(int j=1;j<=cnt&&prime[j]*i<1e7;j++)
        {
            vis[prime[j]*i]=1;
            if(i%prime[j]==0) break;
        }
    }
}
ll dp[maxx];
ll ans[maxx];
ll dou[maxx];
signed main()
{
   int n;
   init();
   //for(int i=1;i<=cnt;i++) cout<<prime[i]<<endl;
   n=read();;
   for(int i=1;i<=n;i++) ans[i]=read(),dou[ans[i]]++;
   
   ll sum=0;
   for(int i=2;i<1e7;i++)
   {
       dp[i]+=dou[i];
       for(int j=1;j<=cnt&&prime[j]*i<1e7;j++)
       {
           dp[i*prime[j]]=max(dp[i*prime[j]],dp[i]);
       }
       sum=max(sum,dp[i]);
   }
   cout<<sum+dou[1]<<endl;
   return 0;
}