质数拆分

质数拆分

题目

本题为填空题，只需要算出结果后，在代码中使用输出语句将所填结果输出即可。

将 2019拆分为若干个两两不同的质数之和，一共有多少种不同的方法？

注意交换顺序视为同一种方法，例如 2 + 2017 = 2019 与 2017 + 2 = 2019 视为同一种方法。

运行限制
最大运行时间：1s
最大运行内存: 128M

埃氏筛(O(nloglogn))
如果$n<10^6$，可以用
一个合数总是可以分解成若干个质数的乘积，换个角度去理解，也就是说合数是某个质数的倍数。此时如果把质数的倍数都去掉，那么剩下的就是质数了。

bool isPrime[10000005];//标记数组 用来表示数字是否是质数 true-是质数 false-不是质数

void aiPrime(int n){// 埃氏筛处理n内的质数
	memset(isPrime,true,sizeof(isPrime));//所有数字，默认标记为质数 
	isPrime[1]=false;//修改1的状态，1不是质数
	for(int i=2;i<=n;i++){//从2开始，枚举范围n内的每个数字 
		if(isPrime[i]){//如果i是质数,将n范围内所有i的倍数，标记为非质数 
			for(int j=2;j*i<=n;j++){ 
				isPrime[i*j]=false;// 标记倍数为非质数 
			}
			
		}
	}
}

线性欧拉筛（O(n)）

int prime[200010];//存放素数
int num=0;
int visit[1000010];//0表示是素数，1表示不是素数
void getprime(int n){
	memset(visit,0,sizeof(visit));
	for(int i=2;i<=n;i++){
		if(visit[i]==0){
			prime[num]=i;
			num++;
		}
		for(int j=0;j<num&&prime[j]*i<=n;j++){
			visit[prime[j]*i]=1;
			if(i%prime[j]==0){
				break;
			}
		}
	}
} 

别问为什么要放质数筛的代码，问就是想复习一遍
……
又是被自己蠢哭的一天
以为是一道很水的题，结果发现自己看错题了/(ㄒoㄒ)/~~
好吧，这道题考的是dp

#include<iostream>
#include<cmath>
#include<cstring>
#include<algorithm> 
using namespace std;
typedef long long ll;
int prime[2500];
int num=1;
int visit[2500];
void getprime(int n){
	memset(visit,0,sizeof(visit));
	for(int i=2;i<=n;i++){
		if(visit[i]==0){
			prime[num]=i;
			num++;
		}
		for(int j=1;j<num&&prime[j]*i<=n;j++){
			visit[prime[j]*i]=1;
			if(i%prime[j]==0){
				break;
			}
		}
	}
} 
ll dp[2500][2500];
int main(){
	getprime(2019);
	memset(dp,0,sizeof(dp));
	dp[0][0]=1;
	for(int i=1;i<num;i++){
		for(int j=0;j<=2019;j++){
			if(j-prime[i]>=0){
				dp[i][j]=dp[i-1][j]+dp[i-1][j-prime[i]];
			}
			else{
				dp[i][j]=dp[i-1][j];
			}
		}
	}
	cout<<dp[num-1][2019]<<endl;
	return 0;
}

这样写也行，我还是比较喜欢第二种写法，不过第一种比较好理解

#include<iostream>
#include<cmath>
#include<cstring>
#include<algorithm> 
using namespace std;
typedef long long ll;
int prime[2500];
int num=1;
int visit[2500];
void getprime(int n){
	memset(visit,0,sizeof(visit));
	for(int i=2;i<=n;i++){
		if(visit[i]==0){
			prime[num]=i;
			num++;
		}
		for(int j=1;j<num&&prime[j]*i<=n;j++){
			visit[prime[j]*i]=1;
			if(i%prime[j]==0){
				break;
			}
		}
	}
} 
ll dp[2500];
int main(){
	getprime(2019);
	memset(dp,0,sizeof(dp));
	dp[0]=1;
	for(int i=1;i<num;i++){
		for(int j=2019;j>=prime[i];j--){
			dp[j]+=dp[j-prime[i]];
		}
	}
	cout<<dp[2019]<<endl;
	return 0;
}

总结一下，就是，这道题感觉出题人是不是在引导着我们看错题啊，给的例子也很误导人啊，算了，是我自己太蠢了