蓝桥杯 算法提高 拿糖果

题意：

问题描述

　　妈妈给小B买了N块糖！但是她不允许小B直接吃掉。
　　假设当前有M块糖，小B每次可以拿P块糖，其中P是M的一个不大于根号下M的质因数。这时，妈妈就会在小B拿了P块糖以后再从糖堆里拿走P块糖。然后小B就可以接着拿糖。
　　现在小B希望知道最多可以拿多少糖。　N <= 100000

首先可以想到的是拿走P块糖之后，问题转化为当前有M-2P块糖，小B最多还能拿多少糖？

这是一个最优子结构问题，递推式为

#include<iostream>
#include<cstring>
using namespace std;
const int N=1e5+10;
int dp[N];
int tag[N];
int prime[N];
int cnt=0;
void getp()
{
    memset(tag,1,sizeof(tag));
    tag[0]=tag[1]=0;
    for(int i=2;i<N;++i)
    {
        if(tag[i])
            prime[cnt++]=i;
        for(int j=0;j<cnt&&i*prime[j]<N;++j)
        {
            tag[i*prime[j]]=0;
            if(i%prime[j])
                break;
        }
    }

}
int main()
{
    getp();
    int n;
    cin>>n;
    memset(dp,0,sizeof(dp));        //最小的素数是2，只有当N>=4的时候才能拿
    for(int i=4;i<=n;++i)
        for(int j=0;prime[j]*prime[j]<=i;++j)
        {
            if(i%prime[j]==0)       //审题不当，说的是M的素因子。。。要是因子qwq
                dp[i]=max(dp[i],dp[i-2*prime[j]]+prime[j]);
        }
//    for(int i=4;i<=n;++i)
//        cout<<i<<' '<<dp[i]<<endl;
    cout<<dp[n]<<endl;


}