saber_jk的博客

一个学生的学习记录

【百炼oj】2706:麦森数

描述形如2p-1的素数称为麦森数,这时P一定也是个素数。但反过来不一定,即如果P是个素数。2p-1不一定也是素数。到1998年底,人们已找到了37个麦森数。最大的一个是P=3021377,它有909526位。麦森数有许多重要应用,它与完全数密切相关。任务:从文件中输入P (1000<P<3100000) ,计算2p-1的位数和最后500位数字(用十进制高精度数表示)输入文件中只包含一个整数P(1000<P<3100000)输出第1行:十进制高精度数2p-1的位数。第2-11行:十进制高精度数2p-1的最后500位数字。(每行输出50位,共输出10行,不足500位时高位补0)不必验证2p-1与P是否为素数。样例输入
1279
样例输出
386
00000000000000000000000000000000000000000000000000
00000000000000000000000000000000000000000000000000
00000000000000104079321946643990819252403273640855
38615262247266704805319112350403608059673360298012
23944173232418484242161395428100779138356624832346
49081399066056773207629241295093892203457731833496
61583550472959420547689811211693677147548478866962
50138443826029173234888531116082853841658502825560
46662248318909188018470682222031405210266984354887
32958028878050869736186900714720710555703168729087
来源联赛复赛试题2003
#include<stdio.h>
#include<string.h>
#define len 500
int result[len],sz[len];
int main(){
    int mult(int*a,int*b);
    int p,i,j;
    scanf("%d",&p);
    printf("%d",(int)(p*log10(2))+1);


    memset(sz,0,sizeof(sz));
    memset(result,0,sizeof(result));
    result[0]=1;
    sz[0]=2;
    while(p>0){
        if(p&1)
            mult(result,sz);
        p>>=1;
        mult(sz,sz);
    }
    result[0]--;
    for(i=len-1;i>=0;i--){
        if((i+1)%50==0)printf("\n");
         printf("%d",result[i]);
        }

}


int mult(int*a,int*b){
    int i,j,k,l;
    int c[len];
    memset(c,0,sizeof(c));
    for(i=0;i<len;i++){
        k=0;
        for(j=0;j<len-i;j++){
            l=c[i+j]+a[j]*b[i]+k;
            c[i+j]=l%10;
            k=l/10;
        }
    }
    for(i=0;i<len;i++)a[i]=c[i];
}
阅读更多
个人分类: 代码
上一篇【百炼oj】4033:铺地毯
下一篇【百炼oj】2818:密码
想对作者说点什么? 我来说一句

没有更多推荐了,返回首页

关闭
关闭