题意:给出一个数,用k个不同的素数组成他,素数不能重复,求有多少种方法
分析:首先打一个表,然后用类似01背包的方法来写
dp[i][j]表示的是数字i用j个数字来组成的方法数,状态转移方程如下
dp[i][j] = sum{dp[i-prime][k-1]};
for(int i=0;i<tot;i++){
for(int j=maxn-1;j>=prime[i];j--){
for(int k=1;k<=14;k++){
dp[j][k] += dp[j-prime[i]][k-1];
}
}
}
关于这里循环的顺序说一下,如果仅看转移方程,dp[6][2]会转移到dp[3][1],但是其实是转移不到的,为了避免这种情况,我们将素数正向枚举(第一层循环),数逆向枚举(第二层循环)。这样的化在计算dp[6][2]的时候dp[3][1]还没有被更新,保证了每个素数只被使用一次。这个思路和01背包的滚动数组类似。
代码:
#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <map>
#include <vector>
#include <queue>
#include <set>
#define FRER() freopen("in.txt","r",stdin)
#define FREW() freopen("out.txt","w",stdout)
#define go int T;cin>>T;for(int kase=1;kase<=T;kase++)
#define debug cout<<"****************"<<endl
#define mod 100000007
using namespace std;
const int maxn = 1121 ,inf = 0x3f3f3f3f;
typedef long long ll;
int n,k;
int tot;
bool vis[maxn];
int prime[maxn];
ll dp[maxn][20];
void cal_prime(){
tot = 0;
for(int i=2;i<maxn;i++){
for(int j=2*i;j<maxn;j+=i)
vis[j] = 1;
}
for(int i=2;i<maxn;i++)
if(!vis[i]) prime[tot++] = i;
}
int main(){
FRER();
FREW();
cal_prime();
memset(dp, 0, sizeof(dp));
dp[0][0] = 1;
for(int i=0;i<tot;i++){
for(int j=maxn-1;j>=prime[i];j--){
for(int k=1;k<=14;k++){
dp[j][k] += dp[j-prime[i]][k-1];
}
}
}
while(~scanf("%d%d",&n,&k)&&n&&k){
printf("%lld\n",dp[n][k]);
}
}