SGU 116 Index of super-prime（数论+dp）

Description
定义超级素数：n是素数且n在素数序列中的位置也是素数，先给出一整数n，问n是否能够被表示成几个超级素数的和，如果可以则输出最少的超级素数数量然后从小到大输出这几个素数，否则输出0
Input
一个不超过10000的整数n
Output
如果n能够被表示成几个超级素数的和则第一行输出组成n的最少的超级素数的数量，第二行升序输出这几个超级素数，否则输出0
Sample Input
6
Sample Output
2
3 3
Solution
筛素数得到素数数组和标记数组，之后将素数数组中下标为素数的素数拿出来放在num数组里表示超级素数，之后就是dp了，dp[i]=min(dp[i-num[j]+1,dp[i])，答案就是dp[n]
Code

#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
#define maxn 10001
#define INF 0x3f3f3f3f
int res,prime[maxn],is_prime[maxn],num[maxn],cnt;
int n,dp[maxn],pre[maxn];
void get_prime()
{
    res=cnt=0;
    memset(is_prime,0,sizeof(is_prime));
    is_prime[0]=is_prime[1]=1;
    for(int i=2;i<maxn;i++)
        if(!is_prime[i])
        {
            prime[++res]=i;
            for(int j=2*i;j<maxn;j+=i)
                is_prime[j]=1;
        }
    for(int i=1;i<=res;i++)
        if(!is_prime[i])
            num[cnt++]=prime[i];
}
int main()
{
    get_prime();
    while(~scanf("%d",&n))
    {
        for(int i=1;i<=n;i++)dp[i]=INF;
        for(int i=1;i<=n;i++)
            for(int j=cnt-1;j>=0;j--)
                if(i-num[j]>=0&&dp[i-num[j]]+1<dp[i])
                {
                    dp[i]=dp[i-num[j]]+1;
                    pre[i]=num[j];
                }
        if(dp[n]!=INF)
        {
            printf("%d\n%d",dp[n],pre[n]);
            for(int i=n-pre[n];i>0;i-=pre[i])
                printf(" %d",pre[i]);
            printf("\n");
        }
        else printf("0\n");
    }
    return 0;
}