poj2706麦森数

最新推荐文章于 2021-06-06 20:13:05 发布

yanerhao

最新推荐文章于 2021-06-06 20:13:05 发布

阅读量886

点赞数

分类专栏：算法区

本文链接：https://blog.csdn.net/yanerhao/article/details/47007091

版权

算法区专栏收录该内容

96 篇文章 6 订阅

订阅专栏

Description
形如2p-1的素数称为麦森数，这时P一定也是个素数。但反过来不一定，即如果P是个素数。2p-1不一定也是素数。到1998年底，人们已找到了37个麦森数。最大的一个是P=3021377，它有909526位。麦森数有许多重要应用，它与完全数密切相关。
任务：从文件中输入P (1000<P<3100000) ,计算2p-1的位数和最后500位数字（用十进制高精度数表示）
Input
文件中只包含一个整数P(1000<P<3100000)
Output
第1行：十进制高精度数2p-1的位数。
第2-11行：十进制高精度数2p-1的最后500位数字。（每行输出50位，共输出10行，不足500位时高位补0）
不必验证2p-1与P是否为素数。
Sample Input
1279
Sample Output
386
00000000000000000000000000000000000000000000000000
00000000000000000000000000000000000000000000000000
00000000000000104079321946643990819252403273640855
38615262247266704805319112350403608059673360298012
23944173232418484242161395428100779138356624832346
49081399066056773207629241295093892203457731833496
61583550472959420547689811211693677147548478866962
50138443826029173234888531116082853841658502825560
46662248318909188018470682222031405210266984354887

32958028878050869736186900714720710555703168729087

代码：

#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#include<string.h>
void mul(int*x,int*y){
int a[126];
memset(a,0,sizeof(a));
int i,j;
for(i=0;i<126;i++)
    {  for(j=0;j<126;j++)
           if(i+j<126)a[i+j]+=x[j]*y[i];
for(j=0;j<125;j++)
	if(a[j]>=10000){a[j+1]+=a[j]/10000;a[j]%=10000;}
}
memcpy(x,a,sizeof(a));}
int main(){
int p;
	int aa[126];
    int aapow[126];
	memset(aa,0,sizeof(aa));
	memset(aapow,0,sizeof(aapow));
	aa[0]=1;
	aapow[0]=2;
	int i,j;
	int k;
	scanf("%d",&p);
	k=(int)p*log10(2.0)+1;
	printf("%d\n",k);
	while(p){
	if(p&1)
	  mul(aa,aapow);
	p>>=1;
	if(p)
	  mul(aapow,aapow);}
	aa[0]--;
	for(i=124;i>=0;i--){
	if(i%25==12)printf("%02d\n%02d",aa[i]/100,aa[i]%100);
	else{printf("%04d",aa[i]);
	  if(i%25==0)printf("\n");}}
return 0;}