2019分解成若干个素数之和的方案数

2019分解成若干个素数之和的方案数

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
bool vis[3005];
ll prime[3005],cnt,dp[3555][3555];
ll n=2019,ans;
//2019分解成若干个素数之和的方案数 
int main() {
	//欧拉线性筛2019内的素数 
	for(int i=2; i<=n; i++) {
		if(!vis[i]) prime[cnt++]=i;
		for(int j=0;j<cnt&&i*prime[j]<=n;j++){
			vis[i*prime[j]]=1;
			if(i%prime[j]==0) break;
		}
	}
//dp[i][j]表示把i分解成若干个素数之和，其中最大素数为j的方案数 
//dp[i][j]=Σdp[i-j][k],(k<=j) 
//	for(int i=0; i<cnt;i++) printf("%d ",prime[i]);
//    printf("\n");
    memset(dp,0,sizeof(dp));
    for(int i=0;i<cnt;i++) dp[prime[i]][prime[i]]=1;
    for(int i=2;i<=n;i++){
    	for(int j=0;prime[j]<=i&&j<cnt;j++){
    		for(int k=0;k<j&&prime[k]<=i-prime[j];k++){
    			dp[i][prime[j]]+=dp[i-prime[j]][prime[k]]; 
			}
		}
	}
	ll ans=0;
	for(int i=0;i<cnt;i++) ans+=dp[n][prime[i]];
	printf("%lld",ans);
}