快速幂的大数运算_快速幂模

快速幂的大数运算_快速幂模快速幂运算1.什么是快速幂2.快速幂的“小数”运算3.高精度(大数)的快速幂1.什么是快速幂快速幂,是指在进行幂运算的时候,用一种快速方法得出答案。比如,要求2^100的值,那按照最简单的方式,就是一个一个2去相乘,然后最终得到答案,那么这样就要计算100次,非常浪费时间,那么快速幂就是使用一种技巧使得将其计算次数减少,快速得到答案。2.快速幂的“小数”运算对于系统内置类型的整型,暂且叫他“小数”,这个时候进行快速幂运算,代码如下:#include<cstdio>#include&l

大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。

Jetbrains全家桶1年46,售后保障稳定

1.什么是快速幂

快速幂,是指在进行幂运算的时候,用一种快速方法得出答案。比如,要求2^100的值,那按照最简单的方式,就是一个一个2去相乘,然后最终得到答案,那么这样就要计算100次,非常浪费时间,那么快速幂就是使用一种技巧使得将其计算次数减少,快速得到答案。

2.快速幂的“小数”运算

对于系统内置类型的整型,暂且叫他“小数”,这个时候进行快速幂运算,代码如下:

#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;
const long long int mod = 1000000000007;   //对答案取模
int main()
{ 
   
	long long int n;
	long long int ans = 1;
	long long int temp = 2;
	cin >> n;   //求2的n次方
	printf("2的%lld次幂对对1000000000007取模的最终值是:", n);
	while (n > 0)  //快速幂模板
	{ 
   
		if (n%2 == 1)
			ans = (ans%mod * temp%mod) % mod;
		n /= 2;
		temp = (temp % mod * temp % mod)%mod;
	}
	printf("%lld",ans%mod);
	return 0;
}

Jetbrains全家桶1年46,售后保障稳定

那么快速幂的原理是什么呢?
用一张图来表示
在这里插入图片描述

3.高精度(大数)的快速幂

上面的代码发现当n的值稍微大一点就不行了,但是用高精度运算就不要有这种限制。下面是代码

#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;

long long int ans[100000];  //存放答案
long long int temp[100000];  //存放2的2^n的值
long long int save[100000];    //临时数组

int ans_len = 1;
int temp_len = 1;
void count_1(long long int* ans, long long int* temp)  //计算数组ans*数组temp,实际上就是简单的高精度大数相乘了
{ 
   
	memset(save, 0, sizeof(save));
	for (int i = 0; i < ans_len; i++)  //先不考虑进制
	{ 
   
		for (int j = 0; j < temp_len; j++)
		{ 
   
			save[i + j] += ans[i] * temp[j];
		}
	}
	int weishu_max = ans_len + temp_len;   //最终答案最终这么多位数。可以那9999*9999这种极端证明
	for (int i = 0; i < weishu_max; i++)
	{ 
   
		save[i + 1] += save[i] / 10;
		save[i] %= 10;
	}
	if (save[weishu_max - 1] != 0)  //n位数*m位数要么就是m+n位数,要么就是m+n-1位数
		ans_len = weishu_max;
	else
	{ 
   
		ans_len = weishu_max - 1;
	}
	for (int i = 0; i < ans_len; i++)  //将save的复制给ans
		ans[i] = save[i];
}


void count_2(long long int* temp)   //算2的2^n次方,也就是temp自乘
{ 
   
	memset(save, 0, sizeof(save));
	for (int i = 0; i < temp_len; i++)
		for (int j = 0; j < temp_len; j++)
			save[i + j] += temp[i] * temp[j];
	int weishu_max = temp_len + temp_len;
	for (int i = 0; i < weishu_max; i++)
	{ 
   
		save[i + 1] += save[i] / 10;
		save[i] %= 10;
	}
	if (save[weishu_max - 1] != 0)
		temp_len = weishu_max;
	else
		temp_len = weishu_max - 1;
	for (int i = 0; i < temp_len; i++)
		temp[i] = save[i];

}

int main()
{ 
   
	int n;
	cin >> n;
	ans[0] = 1;
	temp[0] = 2;
	while (n > 0)
	{ 
   
		if (n % 2 == 1)
		{ 
   
			count_1(ans, temp);
		}
		n /= 2;
		count_2(temp);
	}
	for (int i = ans_len-1; i >=0; i--)
	{ 
   
		cout << ans[i];
	}
	return 0;
}

如果用不考虑进制的做法话,实际上求大数还是有一定的限制,因为一个数组元素最多临时保存了n个数相加,如果当这n个数的和大于long long int的范围,那就会出错。那么可以用考虑进制的做法,虽然麻烦一点,但是会完美无缺,具体做法之前的博客都有提到,可以去看一看。实际上也非常简单,写个乘法竖式总结一下规律就可以写了。

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请联系我们举报,一经查实,本站将立刻删除。

发布者:全栈程序员-站长,转载请注明出处:https://javaforall.net/215486.html原文链接:https://javaforall.net

(0)
全栈程序员-站长的头像全栈程序员-站长


相关推荐

  • 命令行中使用adb安装apk「建议收藏」

    命令行中使用adb安装apk「建议收藏」在你的android—IDE中找到D:\Softwave_Ghost\技术软件\IDE\adt_android_eclipse\adt-bundle-windows-x86-20130219\sdk\platform-tools目录下有adb.exe  cd到这个目录下然后可以进行adb操作 或者设置环境变量 我的电脑->属性->高级系统设置->环境变量->path编辑将上面那个目录复

    2022年5月10日
    29
  • java后端开发需要什么_从事Java后端开发,要学习哪些知识和技能?[通俗易懂]

    java后端开发需要什么_从事Java后端开发,要学习哪些知识和技能?[通俗易懂]很多小伙伴想转行做Java的后端,但是又不知道到底该学习些什么。今天就跟你们聊聊做Java的后端,需要学习和了解什么?1、首先要明确后端包括哪些职业DBA(数据库维护优化专家)Developer(程序猿)Architect(构架师)Scrummaster及类似(敏捷开发专家)ProjectManager(产品狗)Maintenance&ITsupport(通讯和服务器相关)当然这只是一个大…

    2022年7月7日
    21
  • vps自建_七牛云搭建图床

    vps自建_七牛云搭建图床目的:将一些小服务应用部署到免费的serverless/VPS上去环境:0成本实现方式:github+vercel/freewha效果:项目一:个人导航项目二:个人博客项目三:个人音乐服务器:背景:上面的项目以前我都是部署在家里群晖上,或者VPS上,但是FRPC和VPS的流量,延时、运维更新等问题,实际用起来很繁琐,最近两年serverless发展很火,于是就萌生了把他们部署到免费的VPS或者serverless产品上网络上有很多hexo博客部署到vercel、github.

    2025年6月11日
    3
  • 前端模块化-总结_前端模块化规范

    前端模块化-总结_前端模块化规范直接上demo。index.jsx:importReactfrom”react”;importReactDOMfrom”react-dom”;import{Button}from”antd”;importResizeableTablefrom”./ResizeableTable”;import”./styles.css”;constdata1=[{key:0,date:”2018-02-11″,amount:

    2025年6月12日
    2
  • 用LM方法的matlab求解,matlab中LM算法「建议收藏」

    用LM方法的matlab求解,matlab中LM算法「建议收藏」symsabyxreal;f=a*exp(-b*x);Jsym=jacobian(f,[ab]);%拟合用数据。参见《数学试验》,p190,例2data_1=[0.250.511.523468];obs_1=[19.2118.1515.3614.1012.899.327.455.243.01];%2.LM算法%初始猜测sa0=10;b0=0.5…

    2022年10月1日
    4
  • 数组长度计算_c语言计算数组长度的函数

    数组长度计算_c语言计算数组长度的函数(1)sizeof方法:sizeof(数组名)/sizeof(数组类型名)说明:数组占用字节除以数组类型所占字节,结果为数组元素个数(2)strlen说明:strlen,求字符串有效长度方法

    2022年8月5日
    2

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注

关注全栈程序员社区公众号