P1045 麦森数

最新推荐文章于 2022-09-29 20:46:14 发布

原创最新推荐文章于 2022-09-29 20:46:14 发布 · 197 阅读

0 ·

CC 4.0 BY-SA版权

数学同时被 2 个专栏收录

8 篇文章

订阅专栏

高精度

5 篇文章

订阅专栏

本文介绍了一种使用高精度计算和快速幂算法解决特定数学问题的方法。通过定义两个数组f和res来存储中间结果，并利用快速幂的性质进行高效计算。算法首先初始化数组，然后通过双重循环实现乘法运算，避免了直接使用大数运算的复杂性。在每次迭代中，通过调整数组元素来完成进位操作，最后输出结果。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

在这里插入图片描述

#include<cstdio>
#include<cmath>
#include<cstring>
using namespace std;
int f[1001],p,res[1001],sav[1001];//乘法要开两倍长度
void result_1()
{
	memset(sav,0,sizeof(sav));
	for(register int i=1;i<=500;i+=1)
	 for(register int j=1;j<=500;j+=1)
	  sav[i+j-1]+=res[i]*f[j];//先计算每一位的值 不进位
	for(register int i=1;i<=500;i+=1)
	{
		sav[i+1]+=sav[i]/10;
		sav[i]%=10;
	 } 
	memcpy(res,sav,sizeof(res));
	//cstring 库里的赋值函数 把sav的值赋给res 
}
void result_2()
{
	memset(sav,0,sizeof(sav));
	for(register int i=1;i<=500;i+=1)
	 for(int j=1;j<=500;j++)
	  sav[i+j-1]+=f[i]*f[j];
	for(register int i=1;i<=500;i++)
	{
		sav[i+1]+=sav[i]/10;
		sav[i]%=10;
	}
	memcpy(f,sav,sizeof(f));
}
int main()
{
	scanf("%d",&p);
	printf("%d\n",(int)(log10(2)*p+1));
	res[1]=1;
	f[1]=2;//高精度赋初值
	while(p!=0)//快速幂模板
	{
		if(p%2==1) result_1();
		p/=2;
		result_2();
	 } 
	res[1]-=1;
	for(register int i=500;i>=1;i-=1)
	 if(i!=500&&i%50==0) printf("\n%d",res[i]);
	 else printf("%d",res[i]);
	return 0;	
 }