标题:这个程序运行了8个小时,求一亿以内的素数和超级素数.
取消只看楼主
budong998
Rank: 2
等 级:论坛游民
帖 子:25
专家分:15
注 册:2010-7-18
结帖率:75%
已结贴  问题点数:20 回复次数:6 
这个程序运行了8个小时,求一亿以内的素数和超级素数.
大家帮忙看看,我求一亿以内的素数和超级素数.从下午5:30开始,一直到零晨0点30才结束.我的电脑配置也可以.INTER E7500的CPU,G41的板子,前端总线1333,2G 1333的内存.大家看看这个程序还可以怎提高一下效率!
程序代码:
#include <stdio.h>
#include <math.h>
long int prime(long int i)     /*求素数的涵数*/
{
    int x;
    for(x=2;x<=sqrt(i);x++)
    if(i%x==0)return 1;    /*不是素数返回主调涵数1*/
    if(x>sqrt(i))return 0;    /*是素数返回主调涵数0*/
}
long int answer()/*此涵数调用prime()得到结果,并输出*/
{
    long int x;
    printf("\ninput a numbers:\n");
    scanf("\t%ld",&x);
    if(prime(x))             /*从被调涵数得到结果判断真假并输出相应结果*/
        printf("No\t");
    else
        printf("yes\t");
}
long int prime_super(long int x)     /*此涵数判断是否为超级涵数*/
{
    int a,c,b;
    a=prime(x);            /*调用判断素数的子涵数,得到被调涵数返回的0或者1*/
    b=0;
    if(x>19)
    {
        while(a==0&&b==0&&x>9)/*当第一次判断返回为0时,再循环直到被调涵数返回值为1或X小于是2位数。*/
        {
            x=x/10;                    /*将X去掉最后一位。*/
            if(x>1)b+=prime(x);       /*如果每次返回都是0,最高位不是1,那么B就一直是0*/
            else b++;
        }   
        if(b==0&&a==0)                /*所有返回值都为0时,这个数是超级素数,返回主调涵数0*/
        return 0;
    }
    else return 1;
}
long int in_n(long int a)           /*此涵数求N以内的所有超级素数或求N是否为素数*/
{
    long int i=2,x,c,d=0,e=0,y;
    printf("input a number:\n");
    input: scanf("%ld",&x);   
    if(x<=2)         /*如果输入的数小于11时返回重输入*/
    {
        printf("the numbers to small,please input a number:\n");
        goto input;
    }
    if(a==51||a==52||a==53)      /*当主调涵数传来的值为‘3’或‘4’或'5'时,要求求N以内的所有超级素数或素数或两个一起求*/
    {
        while(i<=x)              /*i从3开始,一直到输入的X大为止。*/
        {  
            if(a==51)c=prime_super(i);/*如果主调涵数传来的值为‘3’,要求求N以内的所有超级素数*/
            else c=prime(i);                /*如果主调涵数传来的值不为‘3’,要求求N以内的所有素数*/
            if(!c)                   /*返回为1时做以下程序*/
            {
                if(a==53)     /*如果主调涵数传来的值为'5',要求求N以内所有素数与超级素数.*/
                    { 
                        if(!(c=prime_super(i)))/*就在求素数的同时调用prime_super(i)涵数求超级涵数*/
                        {
                            printf("\b\b>>");d++;      /*用>标级所有素数中的超级素数,并记录个数*/
                        }
                    }
                printf("%ld  \t",i);          /*打印出是素的数*/
                e++;      /*记录多少个数*/
            }
            i++;
        }
        if(a==53)
            printf("\n\tin %ld within super prime numbers sum is %ld ",x,d);        /*打印出超级素数的个数*/
        printf("\n\tin %ld within prime numbers sum is %ld ",x,e);                     /*打印出素数的个数*/
        printf("\nIf necessary validate,press any key.exit press\"q\"\n");/*提示是否需要手动输入验证得到的任何一个素数的正确性。*/
        getchar();                       /*当输入为“Q”时,结束程序;否则继续*/
        while(getchar()!='q')
        answer();                         /*调用子涵数手动输入,去掉最后位验证*/
    }
    if(a==50)                         /*当主调涵数传来的值为‘2’时,要求求1个数是否为超级素数*/
    {
        if(!(c=prime_super(x)))
            printf("\n%ld is super prime Number\n",x);
        else printf("\n%ld no super primeNumber\n",x);
    }
}
main()
{
    int i;
    printf("\nanswer any a numbers whether prime Number,press \"1\";\nanswer any a numbers whether super prime Number,press \"2\";\nanswer in \"n\" within all super prime Number,press \"3\";");
    printf("\nanswer in \"n\" within all prime Number,press \"4\";\nanswer in \"n\" within all prime Number and super prime Number,press \"5\".\n");
    i=getchar();
    switch(i)
    {
        case '1': answer(); break;    /*按1,求任何一个数是否素数*/
        case '2': in_n(i); break;        /*按2,求任何一个数是否超级素数*/
        case '3': in_n(i); break;        /*按3,求任何一个数以内的所有超级素数*/
        case '4': in_n(i); break;        /*按4,求任何一个数以内的所有素数*/
        case '5': in_n(i); break;        /*按5,求任何一个数以内的所有素数及超级素数*/
        default: break;
    }
}
        
搜索更多相关主题的帖子: 素数 超级 运行 
2010-12-11 16:31
budong998
Rank: 2
等 级:论坛游民
帖 子:25
专家分:15
注 册:2010-7-18
得分:0 
回复 5楼 御坂美琴
我本来也是这样的做,先确定素数,然后再拆分确定是不是超级素数.
但是就光求一亿以内的所有素数,也要用两个小时才能求完.
2010-12-11 19:07
budong998
Rank: 2
等 级:论坛游民
帖 子:25
专家分:15
注 册:2010-7-18
得分:0 
回复 7楼 yu_hua
谢谢.你这种算法效率是很高.但在TC下数组付不了那么多个的元素.
2010-12-12 21:57
budong998
Rank: 2
等 级:论坛游民
帖 子:25
专家分:15
注 册:2010-7-18
得分:0 
回复 14楼 yu_hua
非常谢谢你的指导!
今天上班没事,我用你12楼给出的算法,重新改了那个程序.结果效率提高了四倍.同时算出一亿以内的素数和超级素数只用了一个半小时.以下是改写的程序,才学不久,写得很烂,不要见笑.
程序代码:
long int big(int a)
{
    long int sn,b=0,c,i,j,m=1,n,x,k[3000]={0,2};
    printf("input a number:\n");
    input: scanf("%ld",&x);   
        sn=sqrt(x);
    if(sn%2==0)sn+=1;
    for(i=3;i<=sn;i+=2)
    {
        j=prime(i);     /*prime 是常规方法求素数的子涵数,把输入的x的平方根以内的素数求出来放在K数组*/
        if(!j)          /*如果返回值为0是素数*/
        {
            m++;k[m]=i;
        }
    }
    for(i=1;i<=m;i++)
    {
        if(a=='4'&&(!(prime_super(k[i]))))          /*以下几个if语句,是看主调涵数传来的什么值,选择求素数还是超级素数,或两种一起求*/
        {
            printf("%ld\t",k[i]);  b++;
        }
        if(a=='5'&&(!(prime_super(k[i]))))
        {
            printf(">");                           /*如果是两个一起求,就在超级素数前标记*/
            b++;
        }
        if(a=='3'||a=='5')
            printf("%ld\t",k[i]); 
    }
    n=m;
    for(i=sn;i<=x;i+=2)               /*用x平方根到X的所有奇数对刚才数组求的的素数取余*/
    {
        for(j=1;j<=m;j++)
        if(i%k[j]==0)break;
        if(j>m)
        {
            if(a=='4'&&(!(prime_super(i))))       /*和上面几个if语句作用一样*/
            {
                printf("%ld\t",i);  b++;
            }
            if(a=='5'&&(!(prime_super(i))))
            {
                printf(">");
                b++;
            }
            n++;
            if(a=='3'||a=='5')
                printf("%ld\t",i);
        }
    }
    if(a=='4'||a=='5')
        printf("\n\tin %ld within super prime numbers sum is %ld ",x,b);   
    if(a=='3'||a=='5')
        printf("\n\tin %ld within prime numbers sum is %ld ",x,n);    
    printf("\nIf necessary validate,press any key.exit press\"q\"\n");
    while(getch()!='q')
        answer('q');       
}
2010-12-13 23:04
budong998
Rank: 2
等 级:论坛游民
帖 子:25
专家分:15
注 册:2010-7-18
得分:0 
回复 26楼 yu_hua
我用TC,把你写的那个代码拿来运行.运行了可能两分钟,算到78个的时候,程序死掉了.
2010-12-14 18:10
budong998
Rank: 2
等 级:论坛游民
帖 子:25
专家分:15
注 册:2010-7-18
得分:0 
回复 31楼 yu_hua
编译时出错.
2010-12-16 21:30
budong998
Rank: 2
等 级:论坛游民
帖 子:25
专家分:15
注 册:2010-7-18
得分:0 
回复 33楼 yu_hua
才学C语言,刚刚入一点门,连指针都还不会用.你说: TC 下编译连接
tcc  thisfile.c <回车>
我不懂要怎么样做.能不能麻烦你详细的讲讲,真诚的谢谢你这么多天来一直关注我的这个问题!
还有个请求:能把你的邮箱地址或QQ号或MSN给我吗?在你方便时,想给你学习.不懂的问题如果有你讲一讲,我想比我自己傻看几天资料来得快!
我的邮箱:budong998@
2010-12-17 17:37



参与讨论请移步原网站贴子:https://bbs.bccn.net/thread-328706-1-1.html




关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 0.142489 second(s), 9 queries.
Copyright©2004-2024, BCCN.NET, All Rights Reserved