蓝桥杯 拿糖果

  算法提高 拿糖果  

时间限制：1.0s   内存限制：256.0MB

   

问题描述

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

输入格式

　　一个整数N

输出格式

　　最多可以拿多少糖

样例输入

15

样例输出

6

数据规模和约定

　　N <= 100000

http://lx.lanqiao.cn/problem.page?gpid=T376

动态规划  每次都有一些选择 一步一步下去 数据规模是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;
}