历届试题:拿糖果

问题描述
妈妈给小B买了N块糖！但是她不允许小B直接吃掉。
假设当前有M块糖，小B每次可以拿P块糖，其中P是M的一个不大于根号下M的质因数。这时，妈妈就会在小B拿了P块糖以后再从糖堆里拿走P块糖。然后小B就可以接着拿糖。
现在小B希望知道最多可以拿多少糖。
输入格式
一个整数N
输出格式
最多可以拿多少糖
样例输入
15
样例输出
6
数据规模和约定
N <= 100000
 

#include<bits/stdc++.h>
using namespace std;
#define M 100005
int prime[M];

int cnt=0;
int d[M];

bool is(int n){//判断n是不是质数 
	for(int i=2;i<=n/2;i++){
		if(n%i==0)
			return false;
	}
	return true;
}
void init(int n){//找到所有不大于根号下n的质因数
	int k=sqrt(n);
	for(int i=2;i<=k;i++){
		if(is(i)){
			prime[i]=1;
		}
	}	
}
int dfs(int n) {  
    if(d[n]) return d[n];  
    if(n <= 3) return d[n] = 0;  
    int m = sqrt(n);  
    for(int i = 2; i <= m; ++i) {  
        if(prime[i] && n % i == 0) { //质因数   
            d[n] = max(d[n], i + dfs(n - 2*i));  
        }  
    }  
    return d[n];  
}  

int main(){
	int n;
	cin>>n;
	init(n);
	cout<<dfs(n)<<endl; 
	return 0;
}