我也不知道这算什么类型题目 不知道是不是dp 简单点说就是先用欧拉筛枚举10的5次方以内的所有质数 然后对于每个数a[i] 在sqrt(a【i】)的时间内枚举他的所有质因子 同时如果前面出现过的这个质因子的话就累加 同时更新a【i】所有质因子出现过的次数的最大值 之后用一个并列的循环把到a【i】为止的前面所有为止的质因子出现数更新为 上一个并列循环算出的最大值 最后对所有出现过的质因子取最大就好了 注意一下最大值初始化为1 不要设为0
AC代码:
#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <queue>
#include <stack>
using namespace std;
#define maxn 100005
int a[maxn],prime[maxn],primesize,phi[maxn],b[maxn];
bool isprime[maxn];
void getlist(int listsize)
{
primesize=0;
memset(isprime,1,sizeof(isprime));
isprime[1]=false;
for(int i=2;i<=listsize;i++)
{
if(isprime[i]) prime[++primesize]=i;
for(int j=1;j<=primesize&&i*prime[j]<=listsize;j++)
{
isprime[i*prime[j]]=false;
if(i%prime[j]==0) break;
}
}
}
int main()
{
int n;scanf("%d",&n);
getlist(maxn-5);
//if(n==1){cout<<1<<endl;return 0;}
for(int i=1;i<=n;i++)
{
scanf("%d",&a[i]);
}
for(int i=1;i<=n;i++)
{
int flag=0;
int mx_i=0;
if(isprime[a[i]]||a[i]==1)
{
//cout<<a[i]<<'k'<<endl;
b[a[i]]++;
continue;
}
for(int j=1;j<=int(sqrt(a[i]));j++)
{
if(a[i]%j==0)
{
if(isprime[j]) {b[j]++;mx_i=max(mx_i,b[j]);}
if(isprime[a[i]/j]&&a[i]/j!=j) {b[a[i]/j]++;mx_i=max(mx_i,b[a[i]/j]);}
}
}
for(int j=1;j<=sqrt(a[i]);j++)
{
if(a[i]%j==0)
{
if(isprime[j]) {b[j]=max(mx_i,b[j]);}
if(isprime[a[i]/j]&&a[i]/j!=j) {b[a[i]/j]=max(mx_i,b[a[i]/j]);}
}
}
}
int ans=1;
for(int i=1;i<=n;i++)
{
if(isprime[i])
{
ans=max(ans,b[i]);
}
}
printf("%d\n",ans);
}