模板—慢速乘

Posted al-ca

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了模板—慢速乘相关的知识,希望对你有一定的参考价值。

用于模数很大直接乘会爆longlong的情况。

#include<iostream>
#include<cstring>
#include<cstdio>
#define LL long long
using namespace std;
LL a,b,p;
LL mull(LL a,LL b,LL p)

	LL res=0;
	while(b)
	
		if(b&1)res=(res+a)%p;
		a=(a+a)%p;
		b=b>>1;
	
	return res;

signed main()

	cin>>a>>b>>p;
	LL ans=mull(a,b,p);
	cout<<ans<<endl;

好处是只要a,b在long long内,无论它们乘起来多大,都可以做。

劣势是时间复杂度是log(b)的,比起正常乘法来太慢了。

其实有O(1)的方法, 只适用于a*b没有超过long long太多的情况(即a,b并不算大,大概均在10^12左右)

设一个常数t。

令x1=a/t,x2=a%t

    y1=b/t,y2=b%t

显然a*b=(x1*t+x2)*(y1*t+y2)=x1*y1*t*t+x1*y2*t+x2*y1*t+x2*y2。但是感觉并不如慢速乘好用……

以上是关于模板—慢速乘的主要内容,如果未能解决你的问题,请参考以下文章

慢速乘模板(支持负数与零

模板—十进制快速幂

BSGS模板(慢速)

模板快速乘

快速乘模板

快速乘模板