四个素数之和问题.docxVIP

问题描述:欧拉证明素数有无穷多个。但每个整数能表示成4个素数之和吗？输入：每行输入一个整数N(N=10 000 000)，它是你需要把它表示成4个素数之和的数。输出：对于每个输入行都有一个输出行，每行包含符合要求的4个素数。如果该素数不能表示为4个素数之和，那么输出一行不能。也可能有很多组解，任何合理的解都将被接受。输入样例：输出样例:24 3 11 3 7363 7 13 134611 11 17 7代码：#includectime#includemath.h#includeiostreamusingnamespace std;constint MAXLEN=10001;long modular_power1(long a,long r,long m)//求模幂a^r mod m的值{long d,t;d=1;t=a;while(r0){if((r%2)==1)d=(d*t)%m;r/=2;t=t*t%m;}return d;}//素数测试算法代码(Miller-Rabin素数测试算法)bool Miller_Rabin(longlong a,longlong n){longlong d,t,bitnum,bit;bit=130;bitnum=n-1;while((bitbitnum)==0)bit=bit1;d=1;while(bit0){t=d;d=(d*d)%n;if((d==1)(t!=1)(t!=n-1))returntrue;if((bitbitnum)!=0)d=(d*a)%n;bit=bit1;}if(d!=1)returntrue;returnfalse;}void G_maini(int n,int *a,int *b)//用到上面的素数测试算法，此函数用来验证哥德巴赫猜想(每个大//于4的偶数都可以写成两个奇素数之和){int a1,a2,b1,b2,k;if((n3)||(n%2==1))coutErrorendl;elseif(n==4||n==6){*a=n/2;*b=n/2;}elseif(n==8){*a=3;*b=5;}elseif(n==10){*a=3;*b=7;}else{srand(time(NULL));for(k=3;k=n/2;k=k+2)//取奇数k作为素数测试的对象{if(k==3)//3是素数，测试n-3是否为素数{a1=1;while((a1==1)||(a1=n-k-1))a1=(int(rand()%(n-k))+1);//随机选取a1if(Miller_Rabin(a1,n-k))continue;//第次Miller_Rabin素数测试else{b1=1;while((b1==1)||(b1=n-k-1)||b1==a1)b1=(int(rand()%(n-k))+1);//随机选取b1if(Miller_Rabin(b1,n-k))continue;//第次Miller_Rabin素数测试else{*a=k;*b=n-k;break;}}}else{a1=1;a2=1;while((a1==1)||(a1=k-1))a1=(int(rand()%k)+1);//随机选取a1while((a2==1)||(a2=n-k-1))a2=(int(rand()%(n-k))+1);//随机选取a2if(Miller_Rabin(a1,k)||Miller_Rabin(a2,n-k))continue;//第次Miller_Rabin素数测试else{b1=1;b2=1;while((b1==1)||(b1=k-1)||b1==a1)b1=(int(rand()%k)+1);//随机选取b1while((b2==1)||(b2=k-1)||b2==a1)b2=(int(rand()%(n-k))+1);//随机选取b2if(Miller_Rabin(b1,k)||Miller_Rabin(b2,n-k))continue;//第次Miller_Rabin素数测试else{*a=k;*b=n-k;break;}}}}}//return 0;}bool prime(int n)//判断n是否是素数，若是则返回true{int j;if(n==2)returntrue;if(n=1||n%2==0)returnfalse; j=3;while(j=int(sqrt(longdouble(n)))) {if(n%j==0)returnfalse; j+=2; } returntrue;}void OddNumber(int num,int *a,int *b,int *a1,int *a