问题 B: 64位整数乘法
时间限制: 1 Sec 内存限制: 128 MB
提交: 111 解决: 56
[提交] [状态] [讨论版] [命题人:admin]
题目描述
求 a 乘 b 对 p 取模的值,其中 1≤a,b,p≤10^18。
输入
第一行a,第二行b,第三行p。
输出
一个整数,表示a*b mod p的值。
样例输入
2 3 9
样例输出
6
#include <iostream>
#include <cstdio>
#include <cassert>
#include <cstring>
#define LL long long
using namespace std;
LL Bit_64(LL a, LL b, LL c)//位运算64位a*b%c
{
LL ans = 0;
LL sum = 1;//记录结果正负性
if(a<0){
a = -a;
sum = -sum;
}
if(b<0){
b = -b;
sum = -sum;
}
if(a<b) swap(a, b);
while(b){
if(b&1)
ans = (ans+a) % c;
a = (2*a) % c;
b >>= 1;
}
return sum * ans;
}
LL mul(LL x, LL y, LL z)//计算x*y mod z的数值
{
LL a=(long double) x*y/z;//long double +-1.2e^4932
LL ans=x*y-a*z;
ans=(ans+z)%z;
return ans;
}
int main()
{
LL a, b, c;
while(~scanf("%lld%lld%lld", &a, &b, &c)){
assert(Bit_64(a, b, c) > 0);
printf("%lld\n", Bit_64(a, b, c)); // mul(a, b, c)
}
}