64位整数乘法[CH0102]

题目:求$a\ast b\mathrm{m}\mathrm{o}\mathrm{d}p$的值，其中 $1$≤$a$,$b$,$p$≤$10^{18}$
方法一
$快速幂变形，将b用二进制表示，即$$b=c_{k-1}{2}^{k-1}$+$c_{k-1}{2}^{k-2}$ +$...$+$c_0{2}^0$
则$a\ast b=c_{k-1}{2}^{k-1}\ast a$+$c_{k-1}{2}^{k-2}\ast a$ +$...$+$c_0{2}^0\ast a$
$\because a\ast 2^k=(a\ast 2^{k-1})\ast 2$;
$\therefore 只要观察到c_k=1,则ans=ans\ast a\mathrm{m}\mathrm{o}\mathrm{d}p$
$每一次位移，都要a=a\ast a\mathrm{m}\mathrm{o}\mathrm{d}p$
代码

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#define ll long long
using namespace std;
ll mul(ll a,ll b,ll p)
{
	a%=p;ll ans=0;
	while(b)
	{
		if(b&1)ans=(ans+a)%p;
		b>>=1;a=(a<<1)%p;
	}
	return ans;
}
int main()
{
	ll a,b,c;scanf("%lld%lld%lld",&a,&b,&c);
	printf("%lld\n",mul(a,b,c));
	return 0;
}

方法二
用longlong自然溢出
代码

#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
ll mul(ll a,ll b,ll p)
{
	a%=p,b%=p;
	ll c=(long double)a*b/p;//向下取整
	ll ans=(a*b-c*p+p)%p;
	return ans;
}
int main()
{
	ll a,b,p;scanf("%lld%lld%lld",&a,&b,&p);
	printf("%lld\n",mul(a,b,p));
	return 0;
}