快速幂，大数处理，递归优化

最新推荐文章于 2022-04-21 10:36:22 发布
最新推荐文章于 2022-04-21 10:36:22 发布
阅读量846 收藏
点赞数
分类专栏： # 数学-组合 # 递归优化/大数处理
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_26437925/article/details/52103213
版权

快速幂

例如求3^75
O(n)求法，直接将3连续乘75次

75的2进制 1001011
3^75 = 3^64 * 3^8 * 3^2 * 3^1 ，乘法次数将减少

对比程序

#include <cstdio>
#include <cstdlib>
#include <iostream>
#include <vector>
#include <string>
#include <algorithm>
#include <iterator>
#include <set>
#include <map>
#include <ctime>

using namespace std;

const int INF = 0x7fffffff;
const int MIN = -INF - 1;

typedef long long LL;

// 时间复杂度 O(n)
LL f1(LL a, int n)
{
    LL rs = 1;
    for (int i = 0; i < n; i++)
    {
        rs *= a;
    }
    return rs;
}

// 时间复杂度 O(logn)
LL f2(LL x, int n)
{
    int Max = 1000000007;
    LL res = 1;
    while (n>0)
    {
        if (n & 1) // 判断最低位是否是1
            res = (res*x) % Max;
        x = (x*x) % Max;// x变大
        n >>= 1; // 右移一位（或者除以2）
    }
    return res;
}

int main()
{
    int n = 1e7;
    clock_t start, finish;
    double totaltime;
    start = clock();

    cout << f2(1,n) << endl;

    finish = clock();
    //totaltime = (double)(finish - start) / CLOCKS_PER_SEC;
    totaltime = (double)(finish - start) ;
    cout << totaltime << endl;

    //
    start = clock();

    cout << f1(1, n) << endl;

    finish = clock();
    totaltime = (double)(finish - start);
    cout << totaltime << endl;

    return 0;
}

数值的整数次方

题目地址 （牛客网 剑指offer）

https://www.nowcoder.com/practice/1a834e5e3e1a4b7ba251417554e07c00?tpId=13&tqId=11165&tPage=1&rp=1&ru=%2Fta%2Fcoding-interviews&qru=%2Fta%2Fcoding-interviews%2Fquestion-ranking

ac代码
class Solution {
public:
    double Power(double base, int exponent) {
        if(exponent == 0)
            return 1;

        bool flag = true;
        if(exponent < 0)
        {
            flag = false;
            exponent = - exponent;
        }
        bool flag2 = true;
        if(base < 0){
            flag2 = false;
            base = - base;
        }

        double sign = 1.0;
        if(flag2 == false && exponent % 2 == 1)
            sign = -1.0;

        double ans = 1; 
        double base2 = base;
        int n = exponent;
        while(n != 0)
        {
            int tmp = n & 1;
            n = n >> 1;
            if(tmp != 0)
                ans = ans * base2;  
            base2 *= base2;
        }
        if(flag == false)
            return sign * 1.0 / ans;
        return sign * ans;
    }
};

一般的规律（线性代数）

满足如下的线性关系表达式
这里写图片描述

都能写成 递推的矩阵相乘的情况，且矩阵的系数可以由前面的几项推出来，是确定的。

2阶情况
这里写图片描述

3阶情况
这里写图片描述

依次类推。。。

这样就可以根据前面几项 * 矩阵次幂 得到后面的项

例如
这里写图片描述

问题将转化为矩阵的乘幂运算，即快速幂问题，可以将时间复杂度优化为 O(logn)

372. Super Pow （leetcode）

题目地址
https://leetcode.com/problems/super-pow/

ac代码：

const int mod = 1337;

class Solution {
public:

    // a^n % 1337 其中 n < 10
    int npow(int a, int n)
    {
        a = a % mod;
        vector<int> apn(11);
        apn[0] = 1;
        apn[1] = a;
        for (int i = 2; i <= n; i++)
        {
            apn[i] = (apn[i - 1] * a) % mod;
        }
        return apn[n];
    }

    int superPow(int a, vector<int>& b) {
        a = a % mod;

        int len = b.size();
        if (len == 1)
        {
            return npow(a, b[0]);
        }

        int ans = 1;
        int n = a;
        for (int i = len - 1; i >= 0; i--){

            if (b[i] > 0)
            {
                int btmp = npow(n, b[i]);
                ans = (ans * btmp) % mod;
            }
            n = npow(n, 10); //n = n^10  即 a,a^10,a^100,a^1000 
        }

        return ans;
    }
};

50. Pow(x, n)

题目地址
https://leetcode.com/problems/powx-n/

ac代码如下,需要注意特殊数据，正负数等的处理

class Solution {
public:
    double myPow(double x, int n) {
        int INF = 0x7fffffff;
        int MIN = -INF - 1;
        if (n == 0)
            return 1;

        if (n == MIN) // 最大负数特殊处理
        {
            if (x < 0) // x是负数，也将变成正数
                x = -x;
            double ans = 1;
            n = INF;
            while (n != 0)
            {
                if ((n & 1) == 1)
                    ans *= x;
                x = x*x;
                n = n >> 1;// 带符号右移
            }
            ans *= x;
            return 1.0 / ans;
        }
        bool flag = true;
        if (n < 0)
        {
            n = -n;
            flag = false;
        }
        int sign = 1;
        if (x < 0)
        {
            x = -x;
            if ((n & 1) == 1) // 奇数
                sign = -1;
        }

        double ans = 1;
        while (n != 0)
        {
            if ((n & 1) == 1)
                ans *= x;
            x = x*x;
            n = n >> 1;// 带符号右移
        }
        if (flag)
            return sign * ans;
        else
            return sign * (1.0 / ans);
    }
};

斐波拉契加强版

http://www.nowcoder.com/practice/2393c500d43a4293aa7a662274aff4d1?tpId=49&tqId=29343&rp=4&ru=/ta/2016test&qru=/ta/2016test/question-ranking

对于斐波拉契经典问题，我们都非常熟悉，通过递推公式F(n) = F(n - 1) + F(n - 2)，我们可以在线性时间内求出第n项F(n)，现在考虑斐波拉契的加强版，我们要求的项数n的范围为int范围内的非负整数，请设计一个高效算法，计算第n项F(n)。第一个斐波拉契数为F(0) = 1。
给定一个非负整数，请返回斐波拉契数列的第n项，为了防止溢出，请将结果Mod 1000000007。
测试样例：
3
返回：2
12
返回 233

const int INF = 0x7fffffff;
const int MIN = -INF - 1;

const int N = 2;
const int Mod = 1000000007;

typedef struct matrix
{
    long long data[N][N];
    matrix()
    {
        for (int i = 0; i < N; i++)
        for (int j = 0; j < N; j++)
            data[i][j] = 0;
    }
}Ma;

// 矩阵相乘
Ma matrixMulti(Ma m1, Ma m2)
{
    Ma rs;
    for (int i = 0; i < N; i++)
    {
        for (int j = 0; j < N; j++)
        {
            rs.data[i][j] = 0;
            for (int k = 0; k < N; k++)
            {
                rs.data[i][j] = (rs.data[i][j] + m1.data[i][k] * m2.data[k][j]) % Mod;
            }
        }
    }
    return rs;
}

// 矩阵的n次幂
Ma matrixN(Ma m, int n)
{
    Ma rs;
    for (int i = 0; i < N; i++) //rs成为单位矩阵
        rs.data[i][i] = 1;

    while (n > 0){
        if (n & 1)
            rs = matrixMulti(rs, m);
        m = matrixMulti(m, m);
        n >>= 1;
    }

    return rs;
}

class Fibonacci {
public:

    int getNthNumber(int n) { //fabonacci加强版
        // write code here
        if (n == 0 || n == 1)
            return 1;

        Ma m; // fabonacci的二阶矩阵
        m.data[0][0] = 1;
        m.data[0][1] = 1;
        m.data[1][0] = 1;
        m.data[1][1] = 0;

        // [fn fn-1] = [f1 f0] * matrix^(n-1)
        Ma rs = matrixN(m, n - 1);

        long long ans = 0;
        for (int i = 0; i < N; i++)
            ans = (ans + rs.data[i][0]) % Mod;

        return ans;
    }
};

3阶问题，同上

http://www.nowcoder.com/practice/f26698ae586b46428c87a5dbcdae272c?tpId=49&tqId=29345&rp=4&ru=/ta/2016test&qru=/ta/2016test/question-ranking

在农场中，奶牛家族是一个非常庞大的家族，对于家族中的母牛，从它出生那年算起，第三年便能成熟，成熟的母牛每年可以生出一只小母牛。即母牛从出生开始的第三年便能做妈妈。最开始农场只有一只母牛，它从第二年开始生小母牛。请设计一个高效算法，返回第n年的母牛总数，已知n的范围为int范围内的正整数。
给定一个正整数n，请返回第n年的母牛总数，为了防止溢出，请将结果Mod 1000000007。
测试样例：
6
返回：9

const int INF = 0x7fffffff;
const int MIN = -INF - 1;

const int N = 3;
const int Mod = 1000000007;

typedef struct matrix
{
    long long data[N][N];
    matrix()
    {
        for (int i = 0; i < N; i++)
        for (int j = 0; j < N; j++)
            data[i][j] = 0;
    }
}Ma;

// 矩阵相乘
Ma matrixMulti(Ma m1, Ma m2)
{
    Ma rs;
    for (int i = 0; i < N; i++)
    {
        for (int j = 0; j < N; j++)
        {
            rs.data[i][j] = 0;
            for (int k = 0; k < N; k++)
            {
                rs.data[i][j] = (rs.data[i][j] + m1.data[i][k] * m2.data[k][j]) % Mod;
            }
        }
    }
    return rs;
}

// 矩阵的n次幂
Ma matrixN(Ma m, int n)
{
    Ma rs;
    for (int i = 0; i < N; i++) //rs成为单位矩阵
        rs.data[i][i] = 1;

    while (n > 0){
        if (n & 1)
            rs = matrixMulti(rs, m);
        m = matrixMulti(m, m);
        n >>= 1;
    }

    return rs;
}

class Cows {
public:

    int countSum(int n) { //fabonacci加强版
        // write code here
        if (n <= 3)// f1 = 1 f2 = 2 f3 = 3
            return n;

        Ma m; // fabonacci的二阶矩阵
        m.data[0][0] = 1;
        m.data[0][1] = 1;
        m.data[1][2] = 1;
        m.data[2][0] = 1;

        // fn = fn-1 + fn-3
        // [fn fn-1 fn-2] = [f3 f2 f1] * matrix^(n-3)
        Ma rs = matrixN(m, n - 3);


        long long ans = 0;
        for (int i = 0; i < N; i++)
            ans = (ans + (3 - i) * rs.data[i][0]) % Mod;

        return ans;
    }
};

参考学习：

矩阵 快速幂
http://www.cnblogs.com/yan-boy/archive/2012/11/29/2795294.html

POJ-3070Fibonacci（矩阵快速幂求Fibonacci数列）
http://www.cnblogs.com/dongsheng/archive/2013/06/02/3114073.html

master-dragon
关注
专栏目录
光速幂 - 加强版快速幂
青烟绕指柔的博客
02-03 957
考虑一个问题：a ^ b % p。 你可能会说，这不是快速幂板子吗？ 但是如果询问有1e7组呢？ 这时就需要用到光速幂了。 考虑指数b，我们其实可以把b根号分解，然后根号预处理。 模板： #pragma GCC optimize("-Ofast","-funroll-all-loops") #include<bits/stdc++.h> #define int long long u...
快速幂及其相关知识
08-01
因此，快速幂作为一种优化手段，能够显著降低时间复杂度至O(log n)，从而极大地提高了计算效率。 #### 二、公式求幂 公式求幂是最直观的方法，它利用C语言中的`pow()`函数进行计算。`pow()`函数位于`math.h`头文件...
快速幂的大数运算
m0_60343477的博客
04-21 647
快速幂运算1.什么是快速幂2.快速幂的“小数”运算3.高精度（大数）的快速幂 1.什么是快速幂 快速幂，是指在进行幂运算的时候，用一种快速方法得出答案。比如，要求2^100的值，那按照最简单的方式，就是一个一个2去相乘，然后最终得到答案，那么这样就要计算100次，非常浪费时间，那么快速幂就是使用一种技巧使得将其计算次数减少，快速得到答案。 2.快速幂的“小数”运算 对于系统内置类型的整型，暂且叫他“小数”，这个时候进行快速幂运算，代码如下: #include<cstdio> #include&l
幂函数的一次算法改进
技术何
12-08 1243
幂函数的一次算法改进
快速幂（递归）
abc1235454的博客
03-22 637
一个数的整数次幂，是我们在计算中经常用到的，但是怎么可以在 O(log(n)) 的时间内算出结果呢？ 代码框中的代码是一种实现，请分析并填写缺失的代码，求 x^y mod p 的结果。 #include <iostream> #include <cstdio> #include <queue> #include <cstring> #include ...
大数快速幂
林苏泽
10-14 670
快速幂 写一个大数的快速幂吧 就是数字大于long long的 原理其实很简单就是拆分一步步的拆分 ll qpow(ll a,ll b){//这个不能少吧 ll ans=1; while(b){ if(b&1) ans=ans*a%mod; b>>=1; a=a*a%mod; } return ...
分数快速幂计算.pptx
最新发布
05-20
- **大数据处理**: 在分布式系统中，利用快速幂可以有效减少数据处理时间。 - **高性能计算**: 在科学计算和工程模拟中，快速幂是提高计算效率的关键技术之一。 总结而言，分数快速幂不仅是一种高效的计算工具，也...
快速幂、快速乘算法.zip
03-24
快速幂和快速乘是两种高效的计算大整数乘法或幂运算的算法，它们在计算机科学，特别是算法设计和数学计算领域中具有重要的应用。这两种算法都基于分治策略，大大减少了运算次数，从而提高了计算效率。 快速幂算法...
分布式快速幂计算框架.pptx
05-20
#### 二、大数快速幂算法优化 **2.1 基于模数的快速幂计算** - **模运算**: 利用模运算的性质，将大数的幂次计算转化为较小模数下的小数幂次计算。 - **存储开销降低**: 减少中间结果的存储空间，提高计算效率。 -...
快速幂、快速乘算法(1).zip
03-27
快速幂和快速乘是计算机科学中两种高效计算大整数乘法和幂次的方法，尤其在处理大数据量和高精度计算时，它们的重要性不言而喻。这两种算法都是基于分治策略来降低时间复杂度，从传统的指数增长优化为对数时间复杂度...
快速幂取模，大数幂次求模，a^p%m
09-15
本函数输入a,p,m，结果输出为a的p次方对m求模的结果。
大数运算的奇妙技巧之“快速幂”
qq_40253497的博客
03-05 402
快速幂 当遇到题目要求输出返回值要%1E7的时候，基本上，该题最终可以写成ab的形式，可以用“快速幂”方法在O(logn)的时间复杂度求解ab。 核心思想： b可以展开为二进制，以a3次方为例，3 = 0b011。 a3 = (a1)1 * (a2)1 * (a4)0 发现：a1 → a2 → (a2)2→ …，可以通过分治法，不断迭代求得 当b的某位为1时, 它将会纳入最终结果的运算。 上述方法...
快速幂（java）
qq_51131591的博客
03-23 300
快速幂必备的数学知识问题引入1. 递归法2. 非递归实现（扩大底数法） 必备的数学知识 对于 加法有（a+b）%c=（(a%c)+(b%c)）%c 则乘法也有（ab）%c=（(a%c)(b%c)）%c //可以代入数去实验 问题引入 我们用程序去计算2的3次方显然等于8，但2的1000次方呢，显然结果是越界，编译器器会报错或输出一个0 对于大数运算，外面采用快速幂的方式进行降低时间复杂度的运算 例如计算2的n次方并对2019取模的值（n可以取大数） 1. 递归法 核心递归式 a的n次方对c取余 1.为奇数
超大数快速幂模板
青春苦瓜
08-03 261
#include <cstdio> #include <cstring> #include <cmath> #include <iostream> #include <algorithm> #include <string> #include <cstdlib> using namesp...
麦森数 大数加快速幂（Java版）
仰望星空的博客
11-22 603
校赛的时候我们有这道题，那时候不会做，后来百度了题解，做的时候我直接用Java大数暴力，结果超时了，当时也想过用快速幂。。但是我不会用大数的快速幂，最后参考了别人的题解写的。 AC代码： import java.math.BigInteger; import java.util.Scanner; public class Main { public static void main(Strin...
Java算法学习4——分治，记忆化递归与快速幂
Warddamn的博客
03-06 341
一、分治 分治其实就是将一个大问题分解为一个个小问题，然后再将每个子问题求出结果，最后将子问题结果合并为原问题的解。分治问题常常用递归方式来解决 二、递归 递归一定要记住递归体与临界点，下面为递归的两个例子 1.输出1…n的全排列 思路：可以分解成输出1为开头的全排列，2开头的全排列，3开头的全排列……的子问题来解决，然后利用一个hashtable数组，记录这个数字是否以及存在数组P中了，如果没有的话就加入，然后hashtable变为true，最后递归结束就把hashtable变成false，下面例子为生成
自然数的n次方-pow算法的优化。
weixin_43704107的博客
05-30 679
最简单的pow #include <iostream> using namespace std; int main(int argc, char const *argv[]) { int num=0; int n=0; int sum=1; cin >> num; cin >> n; for (int i = 0; i < n; ++...
大数的快速幂（模拟水题。）
永远鲜红の幼月的博客
02-20 1336
题目链接：http://noi.openjudge.cn/ch0204/2991/描述已知长度最大为200位的正整数n，请求出2011^n的后四位。输入第一行为一个正整数k，代表有k组数据，k&lt;=200接下来的k行，每行都有一个正整数n，n的位数&lt;=200输出每一个n的结果为一个整数占一行，若不足4位，去除高位多余的0样例输入3 5 28 792样例输出1051 81 5521这道题其...
快速幂取模算法详解及优化
"快速幂算法是一种高效的计算大数幂并取模的算法，常用于解决计算机科学中的数学问题，特别是在算法竞赛和数据结构的学习中。它解决了常规幂运算容易溢出且效率低下的问题。快速幂算法的核心是利用指数的二进制表示...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值