快速乘法

最新推荐文章于 2022-03-30 15:55:32 发布

斌同学斌同學

最新推荐文章于 2022-03-30 15:55:32 发布

阅读量218

点赞数

分类专栏：数论文章标签： acm 数论 math

本文链接：https://blog.csdn.net/GreatJames/article/details/76847162

版权

数论专栏收录该内容

34 篇文章 0 订阅

订阅专栏

/* ***********************************************
┆  ┏┓　　　┏┓ ┆
┆┏┛┻━━━┛┻┓ ┆
┆┃　　　　　　　┃ ┆
┆┃　　　━　　　┃ ┆
┆┃　┳┛　┗┳　┃ ┆
┆┃　　　　　　　┃ ┆
┆┃　　　┻　　　┃ ┆
┆┗━┓　马　┏━┛ ┆
┆　　┃　勒　┃　　┆　　　　　　
┆　　┃　戈　┗━━━┓ ┆
┆　　┃　壁　　　　　┣┓┆
┆　　┃　的草泥马　　┏┛┆
┆　　┗┓┓┏━┳┓┏┛ ┆
┆　　　┃┫┫　┃┫┫ ┆
┆　　　┗┻┛　┗┻┛ ┆
************************************************ */
#include <iostream>
#include <set>
#include <map>
#include <stack>
#include <cmath>
#include <queue>
#include <cstdio>
#include <bitset>
#include <string>
#include <vector>
#include <iomanip>
#include <cstring>
#include <algorithm>
#include <functional>
#define PI acos(-1)
#define eps 1e-8
#define inf 0x3f3f3f3f
#define debug(x) cout<<"---"<<x<<"---"<<endl
typedef long long ll;
using namespace std;


const long long M = 1000000007;
long long quickmul_1(long long a, long long b)
{
    long long ret = 0;
    for (; b; b >>= 1, a = (a << 1) % M)
        if (b & 1)
        {
            ret = (ret + a) % M;
        }
    return ret;
}

///What the fuck is this??!! It works!!
long long quickmul_2(long long a, long long b)
{
    return (a * b - (long long)(a / (long double)M * b + 1e-3) * M + M) % M;
}

int main()
{
    int a, b;
    while (cin >> a >> b)
    {
        cout << quickmul_1(a, b) << endl;
        cout << quickmul_2(a, b) << endl;
    }
    return 0;
}


/************************************************
┆  ┏┓　　　┏┓ ┆
┆┏┛┻━━━┛┻┓ ┆
┆┃　　　　　　　┃ ┆
┆┃　　　━　　　┃ ┆
┆┃　┳┛　┗┳　┃ ┆
┆┃　　　　　　　┃ ┆
┆┃　　　┻　　　┃ ┆
┆┗━┓　  　┏━┛ ┆
┆　　┃　  　┃　　┆　　　　　　
┆　　┃　  　┗━━━┓ ┆
┆　　┃　 AC代马 　　┣┓┆
┆　　┃　        　　┏┛┆
┆　　┗┓┓┏━┳┓┏┛ ┆
┆　　　┃┫┫　┃┫┫ ┆
┆　　　┗┻┛　┗┻┛ ┆
************************************************ */