UVA Sum of Different Primes (数论+dp)

题意：给出一个数，用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]);
    }
}