0x01位运算——64位整数乘法

0x01 位运算——64位整数乘法

• 问题简析

C++内置的最高整数类型是64位，所以显然不能通过直接运算得到答案，我们需要一些特殊的处理方法，在这里提供两种方法。

1. 思路一

类似于快速幂的思想，把b拆分成k位二进制。
$b=C_{k-1}\ast 2^{k-1}+C_{k-2}\ast 2^{k-2}...+C_0\ast 2^0。$
$a\ast b=C_{k-1}\ast 2^{k-1}\ast a+C_{k-2}\ast 2^{k-2}\ast a...+C_0\ast 2^0\ast a。$
$若算出a\ast b^{i-1}modp，则计算(a\ast b^{i-1})\ast 2modp的过程中不会超过2\ast 10^{18}，longlong可以胜任。$

• 时间复杂度分析

$分解出的项数为k=\lceil lo{g}_2(b+1)\rceil$
$时间复杂度:Olo{g}_2(n)$

• 代码实现

#include<bits/stdc++.h>

typedef long long ll;
typedef unsigned long long ull;

const int N = 1e5+10,M = N,INF = 0x3f3f3f3f,P = 998244353;


int main()
{
    std::ios::sync_with_stdio(false);
    ll a,b,p;
    std::cin>>a>>b>>p;
    ll res = 0;
    for(; b ; b >>= 1)
    {
        if(b&1)res = (res + a)%p;
        a = a * 2 % p;
    }
    std::cout<<res;
    return 0;
}

2. 思路二

$利用a\ast bmodp=a\ast b-\lfloor a\ast b/p\rfloor \ast p，\lfloor \rfloor 为向下取整$
$首先，假设a<p,b<p,那么a\ast b/p下取整之后一定也小于p，我们可以使用浮点数完成中间计算。$
$longdouble的浮点数有效数字在18$~$19位之间，可以胜任。$
$a\ast b和\lfloor a\ast b/p\rfloor \ast p可能会很大,但是易知它们的差值一定在0$~$p-1之间。$
$因为p的范围在longlong以内，所以只需要关心a\ast b和\lfloor a\ast b/p\rfloor \ast p在longlong范围内的位，更高的位一定是相同的。$

• 时间复杂度分析

$O(1)$

• 代码实现

#include<bits/stdc++.h>

typedef long long ll;
typedef unsigned long long ull;

const int N = 1e5+10,M = N,INF = 0x3f3f3f3f,P = 998244353;


int main()
{
    std::ios::sync_with_stdio(false);
    ll a,b,p;
    std::cin>>a>>b>>p;
    a %= p,b %= p;
    ll c = (long double)a * b / p;
    ll res = a * b - c * p;
    if(res < 0)res += p;
    else if(res >= p)res -= p;
    std::cout<<res;
    return 0;
}

• 注意点

1. 为什么a和b要mod p？
   因为我们需要保证$a\ast b/p$范围在long long内
2. 为什么用别的浮点数类型(double)会错？
   因为double的有效位数过小，在乘p以后导致错误
   例子：0.002 * 10000和0.0018 * 10000显然不同
   精度足够以后：
   0.0018 * 10000 和0.001811 * 10000不会影响
3. 为什么最后要有这一段代码？

if(res < 0)res += p;
else if(res >= p)res -= p;

因为在第一行对a和b mod p后，由模的性质可以知道最后答案可能会与原来的答案产生p的差值。

• 题目：
  64位整数乘法