64位乘法

来源：AcWing 90. 64位整数乘法

题目

求 a 乘 b 对 p 取模的值。

输入格式

第一行输入整数a，第二行输入整数b，第三行输入整数p。

输出格式

输出一个整数，表示a*b mod p的值。

数据范围

1≤a,b,p≤10¹⁸

输入样例：

3
4
5

输出样例：

2

题解

思路

指数可用二进制表示，即若干个2的次幂之和，a可单独乘每个2的次幂，然后再相加得到答案。

a * 2 ⁱ = ( a * 2 ^i-1 ) * 2

其中 a * 2ⁱ 为 x

代码

#include<bits/stdc++.h>
using namespace std;
long long a,b,p;
long long m_mul(long long x,long long y,long long z)
{
	long long ans=0;
	for(;y;y>>=1)
	{
		if(y&1) ans=(ans+x)%z;
		x=x*2%z;
	}
	return ans;
}
inline long long mr()
{
	long long x=0;char ch=getchar();
	while(ch<'0'||ch>'9') ch=getchar();
	while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+(ch^48);ch=getchar();}
	return x;
}
int main()
{
	a=mr();
	b=mr();
	p=mr();
	printf("%lld",m_mul(a,b,p));
	return 0;
}