算法提高 拿糖果
时间限制:1.0s 内存限制:256.0MB
问题描述
妈妈给小B买了N块糖!但是她不允许小B直接吃掉。
假设当前有M块糖,小B每次可以拿P块糖,其中P是M的一个不大于根号下M的质因数。这时,妈妈就会在小B拿了P块糖以后再从糖堆里拿走P块糖。然后小B就可以接着拿糖。
现在小B希望知道最多可以拿多少糖。
假设当前有M块糖,小B每次可以拿P块糖,其中P是M的一个不大于根号下M的质因数。这时,妈妈就会在小B拿了P块糖以后再从糖堆里拿走P块糖。然后小B就可以接着拿糖。
现在小B希望知道最多可以拿多少糖。
输入格式
一个整数N
输出格式
最多可以拿多少糖
样例输入
15
样例输出
6
数据规模和约定
N <= 100000
动态规划 每次都有一些选择 一步一步下去 数据规模是10万每次最多也就几十个质因子
剩下糖果数一样的时候选择小B拿到最多的就行
#include<bits/stdc++.h>
using namespace std;
int dp[100010],bj[100010],vis[100010];
int ans;
struct node
{
int x,y;
friend bool operator <(node a,node b)
{
return a.x<b.x;
}
} ft,et;
int pd(int h)
{
for(int i=2; i*i<=h; i++)
if(h%i==0) return 0;
return 1;
}
void solve(int n)
{
ans=0;
priority_queue<node>qe;
ft.x=n,ft.y=0;
qe.push(ft);
while(!qe.empty())
{
ft=qe.top();
qe.pop();
if(ft.y<dp[ft.x])
continue;
if(vis[ft.x]||ft.x<=1)//如果是质数或者小于一就只能停止了
{
ans=max(ans,ft.y);
continue;
}
for(int i=1; bj[i]*bj[i]<=ft.x; i++)//
{
if(ft.x%bj[i]==0)//bj[i]不能为零否则运行出错
{
et.x=ft.x-2*bj[i];
et.y=ft.y+bj[i];
if(dp[et.x]<et.y)
{
dp[et.x]=et.y;
qe.push(et);
}
}
}
}
}
int main()
{
memset(bj,0,sizeof(bj));
memset(vis,0,sizeof(vis));
int k=1;
for(int i=2; i*i<=101000; i++)//多一个质数出来
{
if(pd(i))
{
bj[k++]=i;
vis[i]=1;//记录质数
}
}
int n;
while(scanf("%d",&n)!=EOF)
{
memset(dp,0,sizeof(dp));
solve(n);
printf("%d\n",ans);
}
return 0;
}