找回密码
 注册
搜索
热搜: 回贴
  • 前程无忧官网首页 有什么好的平台可以
  • 最新的销售平台 互联网营销的平台有哪
  • 制作网页的基本流程 网页制作和网页设
  • 【帝国CMS】输出带序号的列表(数字排
  • 网站建设公司 三一,中联,极东泵车的
  • 织梦 建站 织梦网站模版后台怎么更改
  • 云服务官网 哪些网站有免费的简历模板
  • 如何建网站要什么条件 建网站要用什么
  • 吉林市移动公司电话 吉林省退休人员网
  • 设计类毕业论文 网站设计与实现毕业论
查看: 1262|回复: 2

[求助]怎么编这个歌德巴赫猜想

[复制链接]
发表于 2009-10-31 02:06:51 | 显示全部楼层 |阅读模式 IP:江苏扬州
四、编程验证哥德巴赫猜想之一:任何大于5的奇数都可表示为3个素数之和(且不唯一)。请编程验证从9到999999间的所有奇数是否符合这个猜想。
发表于 2009-10-31 02:06:52 | 显示全部楼层 IP:江苏扬州
我只能证到99999再大一位运算量太惊人了.
#include <iostream.h> #include <stdlib.h> main() { long i,j[20000]={0},x; int m=0,y=1; cout<<"找出全部素数..."<<endl; for(x=1;x<=99999;x+=2) { y=1; for(i=2;i<=x/2;i++) if(x%i==0) {y=0;break;} if(y==1) { j[m]=x; m++; } } m=0; cout<<"开始验证哥德巴赫猜想..."<<endl; for(x=9;x<=99999;x+=2) { for(int z=0;j[z]<x;z++) for(int q=z+1;j[q]<x;q++) for(int w=q+1;j[w]<x;w++) if(j[z]+j[q]+j[w]==x) { //cout<<j[z]<<"\t"<<j[q]<<"\t"<<j[w]<<"\t"<<x<<endl; m=1; goto _NEXT; } _NEXT: if(m!=1) { cout<<x<<endl; cout<<"不成立!"<<endl; exit(0); } } cout<<"成立!"<<endl; }
回复

使用道具 举报

发表于 2009-10-31 02:06:53 | 显示全部楼层 IP:江苏扬州
这样行吗? #include <iostream.h>
long NextPrime(long int * s, long max, long v) //取得比v大的下一个素数 { for(long n=v+1;n<max;n++) if(s[n]!=0) return n; return 0; }
int main() { long int max; cout<<"输入MAX(如999999)="; cin>>max;
long int * s =new long int[max]; //标记2-max之间的所有素数 for(long int i=0;i<max;i++) {s[i]=i;}
for(long int n=2;n<max;n++) { for(long pos=2*n;pos<max;pos+=n) { s[pos]=0; } }
//验证猜想 for (long value=7;value<max;value+=2) { long v1=3; long v2=2,v3=value-v1-v2; while(v2>0 && v2<value && v3>=v2 && v3>1) { if(v3>1 && s[v3]!=0) { cout<<value<<"="<<v1<<"+"<<v2<<"+"<<v3<<endl; break; } v2=NextPrime(s,max,v2); v3=value-v1-v2; } }
delete[] s; }
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

QQ|小黑屋|最新主题|手机版|微赢网络技术论坛 ( 苏ICP备08020429号 )

GMT+8, 2024-9-29 17:24 , Processed in 0.232860 second(s), 13 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5

© 2001-2023 Discuz! Team.

快速回复 返回顶部 返回列表