紫书p343
令d[i][j]为数i分成j个不同素数的个数 d[i][j]+=d[i-prime[k]][j-1] prime[k]<i
因为不允许有重复元素,所以计算的顺序就很关键,先从各个素数从小到大开始枚举,然后i和j都是逆序枚举,如果不逆序的话,就可能导致重复出现,类似01背包。
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <vector>
#include <map>
#include <cmath>
#include <set>
#include <queue>
using namespace std;
const int INF=1e9+10;
typedef long long ll;
bool isprime[2000];//用true表示i是素数
int prime[2000];//第i个素数
int get_prime(int n){
int p=0;
for(int i=0;i<=n;i++) isprime[i]=true;
isprime[0]=isprime[1]=false;
for(int i=2;i<=n;i++){
if(isprime[i]){
prime[p++]=i;
for(int j=2*i;j<=n;j+=i) isprime[j]=false;
}
}
return p;
}
ll d[1500][20];
int main(){
int p=get_prime(1500);
int n,m;
d[0][0]=1LL;
for(int k=0;prime[k]<1120;k++){
for(int j=14;j>=1;j--){
for(int i=1120;i>=1;i--){
if(prime[k]>i) continue;
d[i][j]+=1LL*d[i-prime[k]][j-1];
}
}
}
while(scanf("%d %d",&n,&m)&&n&&m){
printf("%lld\n",d[n][m]);
}
return 0;
}