UPC：位运算64位整数乘法

问题 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)
    }
}