G - G CodeForces - 264B

我也不知道这算什么类型题目 不知道是不是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);
}