快速幂/快速幂取模/矩阵求解快速幂

2024-05-29 19:18
文章标签 快速 矩阵 求解 取模

本文主要是介绍快速幂/快速幂取模/矩阵求解快速幂,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

因为做了一个Fibonacci的题需要用矩阵+快速幂求解所以引发了下面一系列问题!!!

快速幂:

            快速幂时间复杂度为 O(log2N), 与朴素的O(N)相比效率有了极大的提高。

            例如:

3 ^ 999 = 3 * 3 * 3 * … * 3 
直接乘要做998次乘法。但事实上可以这样做: 
3 ^ 2 = 3 * 3 
3 ^ 4 = (3 ^ 2) * (3 ^ 2) 
…………
3 ^ 256 = (3 ^ 128) * (3 ^ 128) 
3 ^ 512 = (3 ^ 256) * (3 ^ 256) 
再相乘: 
3 ^ 999 
= 3 ^ (512 + 256 + 128 + 64 + 32 + 4 + 2 + 1) 
= (3 ^ 512) * (3 ^ 256) * (3 ^ 128) * (3 ^ 64) * (3 ^ 32) * (3 ^ 4) * (3 ^ 2) * 3 
这样只要做16次乘法。
代码中k%2的意思是如果当前这一位的二进制是 1 还是 0 ,如果是 1 ,则表示将其乘 ans (即加入最终结果),否则不乘 ans (即不加入最终结果)。

 代码:

#include <stdio.h>
int power(int n, int k)
{int ans = 1;while( k ){if(k&1)ans *= n ;k= k/2;n *= n;}return ans ;
}
int main ( ) 
{int n, k;while(scanf("%d %d",&n, &k)!=EOF){printf("%d\n",power(n,k));}
}

快速幂取模:

       快速幂取模还不算太懂,只知道需要对每一个值取模,公式 a*b%c=((a%c)*b)%c

代码:

#include<stdio.h>
int main()
{long long a,b,c,d;int T;scanf("%d",&T);while(T--){scanf("%lld%lld%lld",&a,&b,&c);long long  ans=1;while( b ){if( b%2 )ans = ( ans * a ) % c;b/=2;a = ( a * a ) % c;}printf("%lld\n",ans);}return 0;
}


矩阵+快速幂 求解Fibonacci数列:

           矩阵:矩阵可以看成一个n×m的数表,用二维数组表示
           矩阵乘法:定义矩阵A,B。A和B可以乘法操作当且仅当A的大小是a×b,B的大小是b×c,设矩阵C=AB,则C的大小是a×c,

        且有:

   

            最普通的矩阵乘法是直接三个for循环直接计算而已,所以复杂度是O(n3).

        

            
  运用矩阵乘法快速幂,可以快速计算出矩阵B^(n-1),这样实现将时间复杂度降低到O(log n).

代码(nyoj 698题 A Coin Problem 别人优代码 ):

#include<stdio.h>
#include<iostream>
#include<algorithm>
#include<string>
#include<cstring>
#include<map>
#include<vector>
using namespace std;
const int Max = 10000;
struct M
{long long p[2][2];
};
M mult(M a,M b)
{int i,j,k;M c;for(i=0;i<2;i++)for(j=0;j<2;j++){c.p[i][j]=0;for(k=0;k<2;k++){c.p[i][j]=(c.p[i][j]+a.p[i][k]*b.p[k][j]%Max)%Max;}}return c;
}
M pow(M a,long long k)
{M b={1,0,0,1};while(k){if(k&1)b=mult(b,a);a=mult(a,a);k>>=1;}return b;
}
int main()
{int t;long long n;scanf("%d",&t);while(t--){scanf("%lld",&n);M a={1,1,1,0};M c;c=pow(a,n-1);printf("%lld\n",(c.p[0][0]*2+c.p[0][1]*1)%Max);}return 0;
}


还不明白请点击~>
 

这篇关于快速幂/快速幂取模/矩阵求解快速幂的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1014411

相关文章

一文详解如何在idea中快速搭建一个Spring Boot项目

《一文详解如何在idea中快速搭建一个SpringBoot项目》IntelliJIDEA作为Java开发者的‌首选IDE‌,深度集成SpringBoot支持,可一键生成项目骨架、智能配置依赖,这篇文... 目录前言1、创建项目名称2、勾选需要的依赖3、在setting中检查maven4、编写数据源5、开启热

C/C++中OpenCV 矩阵运算的实现

《C/C++中OpenCV矩阵运算的实现》本文主要介绍了C/C++中OpenCV矩阵运算的实现,包括基本算术运算(标量与矩阵)、矩阵乘法、转置、逆矩阵、行列式、迹、范数等操作,感兴趣的可以了解一下... 目录矩阵的创建与初始化创建矩阵访问矩阵元素基本的算术运算 ➕➖✖️➗矩阵与标量运算矩阵与矩阵运算 (逐元

MybatisX快速生成增删改查的方法示例

《MybatisX快速生成增删改查的方法示例》MybatisX是基于IDEA的MyBatis/MyBatis-Plus开发插件,本文主要介绍了MybatisX快速生成增删改查的方法示例,文中通过示例代... 目录1 安装2 基本功能2.1 XML跳转2.2 代码生成2.2.1 生成.xml中的sql语句头2

8种快速易用的Python Matplotlib数据可视化方法汇总(附源码)

《8种快速易用的PythonMatplotlib数据可视化方法汇总(附源码)》你是否曾经面对一堆复杂的数据,却不知道如何让它们变得直观易懂?别慌,Python的Matplotlib库是你数据可视化的... 目录引言1. 折线图(Line Plot)——趋势分析2. 柱状图(Bar Chart)——对比分析3

一文教你Java如何快速构建项目骨架

《一文教你Java如何快速构建项目骨架》在Java项目开发过程中,构建项目骨架是一项繁琐但又基础重要的工作,Java领域有许多代码生成工具可以帮助我们快速完成这一任务,下面就跟随小编一起来了解下... 目录一、代码生成工具概述常用 Java 代码生成工具简介代码生成工具的优势二、使用 MyBATis Gen

使用animation.css库快速实现CSS3旋转动画效果

《使用animation.css库快速实现CSS3旋转动画效果》随着Web技术的不断发展,动画效果已经成为了网页设计中不可或缺的一部分,本文将深入探讨animation.css的工作原理,如何使用以及... 目录1. css3动画技术简介2. animation.css库介绍2.1 animation.cs

SpringBoot快速搭建TCP服务端和客户端全过程

《SpringBoot快速搭建TCP服务端和客户端全过程》:本文主要介绍SpringBoot快速搭建TCP服务端和客户端全过程,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,... 目录TCPServerTCPClient总结由于工作需要,研究了SpringBoot搭建TCP通信的过程

一文教你Python如何快速精准抓取网页数据

《一文教你Python如何快速精准抓取网页数据》这篇文章主要为大家详细介绍了如何利用Python实现快速精准抓取网页数据,文中的示例代码简洁易懂,具有一定的借鉴价值,有需要的小伙伴可以了解下... 目录1. 准备工作2. 基础爬虫实现3. 高级功能扩展3.1 抓取文章详情3.2 保存数据到文件4. 完整示例

快速修复一个Panic的Linux内核的技巧

《快速修复一个Panic的Linux内核的技巧》Linux系统中运行了不当的mkinitcpio操作导致内核文件不能正常工作,重启的时候,内核启动中止于Panic状态,该怎么解决这个问题呢?下面我们就... 感谢China编程(www.chinasem.cn)网友 鸢一雨音 的投稿写这篇文章是有原因的。为了配置完

Python利用ElementTree实现快速解析XML文件

《Python利用ElementTree实现快速解析XML文件》ElementTree是Python标准库的一部分,而且是Python标准库中用于解析和操作XML数据的模块,下面小编就来和大家详细讲讲... 目录一、XML文件解析到底有多重要二、ElementTree快速入门1. 加载XML的两种方式2.