快速幂（2进制）

最新推荐文章于 2023-07-29 22:07:29 发布

xiaojuruo_haokee

最新推荐文章于 2023-07-29 22:07:29 发布

阅读量176

点赞数

分类专栏：好渴鹅学数学文章标签： c++ 开发语言

本文链接：https://blog.csdn.net/xiaojuruo_haokee/article/details/126422291

版权

好渴鹅学数学专栏收录该内容

1 篇文章 0 订阅

订阅专栏

先看代码：

inline unsigned long long power(unsigned long long a,unsigned long long b)
{
	unsigned long long sum = 1,base = a;
	while(b)
	{
		if(b&1==1)
		{
			(sum *= base) %= p;
		}
		base = base*base % p;
		b = b>>1;
	}
	return sum;
}

这就是快速幂。它的原理很简单，就是求出指数的二进制的最后一位，来判断是不是偶数，然后右移一位把二进制位的最后一位抛掉。这样就能找到指数的二进制的每一位，最后进行运算即可。。。

P3197 [HNOI2008]越狱 这一题就是用的快速幂来运算的，但这题 $m^{n}$ 求出来要减去不能越狱的情况，即 $m×m-1^{n-1}$ ,最后输出即可。。。（~~好像漏了什么~~）

Code:

#include<bits/stdc++.h>
using namespace std;
const unsigned long long p=100003;
inline unsigned long long power(unsigned long long a,unsigned long long b);

signed main(signed argc,char *argv[])
{
	unsigned long long m,n;
	cin >> m >> n;
	cout << (power(m,n) % p - (m * power(m-1,n-1)) % p + p )% p << endl; 
	return 0;
}

inline unsigned long long power(unsigned long long a,unsigned long long b)
{
	unsigned long long sum = 1,base = a;
	while(b)
	{
		if(b&1==1)
		{
			(sum *= base) %= p;
		}
		base = base*base % p;
		b = b>>1;
	}
	return sum;
}