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

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

相关文章

使用EasyPoi快速导出Word文档功能的实现步骤

《使用EasyPoi快速导出Word文档功能的实现步骤》EasyPoi是一个基于ApachePOI的开源Java工具库,旨在简化Excel和Word文档的操作,本文将详细介绍如何使用EasyPoi快速... 目录一、准备工作1、引入依赖二、准备好一个word模版文件三、编写导出方法的工具类四、在Export

Python绘制TSP、VRP问题求解结果图全过程

《Python绘制TSP、VRP问题求解结果图全过程》本文介绍用Python绘制TSP和VRP问题的静态与动态结果图,静态图展示路径,动态图通过matplotlib.animation模块实现动画效果... 目录一、静态图二、动态图总结【代码】python绘制TSP、VRP问题求解结果图(包含静态图与动态图

Python多线程实现大文件快速下载的代码实现

《Python多线程实现大文件快速下载的代码实现》在互联网时代,文件下载是日常操作之一,尤其是大文件,然而,网络条件不稳定或带宽有限时,下载速度会变得很慢,本文将介绍如何使用Python实现多线程下载... 目录引言一、多线程下载原理二、python实现多线程下载代码说明:三、实战案例四、注意事项五、总结引

C#使用Spire.XLS快速生成多表格Excel文件

《C#使用Spire.XLS快速生成多表格Excel文件》在日常开发中,我们经常需要将业务数据导出为结构清晰的Excel文件,本文将手把手教你使用Spire.XLS这个强大的.NET组件,只需几行C#... 目录一、Spire.XLS核心优势清单1.1 性能碾压:从3秒到0.5秒的质变1.2 批量操作的优雅

Mybatis-Plus 3.5.12 分页拦截器消失的问题及快速解决方法

《Mybatis-Plus3.5.12分页拦截器消失的问题及快速解决方法》作为Java开发者,我们都爱用Mybatis-Plus简化CRUD操作,尤其是它的分页功能,几行代码就能搞定复杂的分页查询... 目录一、问题场景:分页拦截器突然 “失踪”二、问题根源:依赖拆分惹的祸三、解决办法:添加扩展依赖四、分页

c++日志库log4cplus快速入门小结

《c++日志库log4cplus快速入门小结》文章浏览阅读1.1w次,点赞9次,收藏44次。本文介绍Log4cplus,一种适用于C++的线程安全日志记录API,提供灵活的日志管理和配置控制。文章涵盖... 目录简介日志等级配置文件使用关于初始化使用示例总结参考资料简介log4j 用于Java,log4c

使用Redis快速实现共享Session登录的详细步骤

《使用Redis快速实现共享Session登录的详细步骤》在Web开发中,Session通常用于存储用户的会话信息,允许用户在多个页面之间保持登录状态,Redis是一个开源的高性能键值数据库,广泛用于... 目录前言实现原理:步骤:使用Redis实现共享Session登录1. 引入Redis依赖2. 配置R

Linux如何快速检查服务器的硬件配置和性能指标

《Linux如何快速检查服务器的硬件配置和性能指标》在运维和开发工作中,我们经常需要快速检查Linux服务器的硬件配置和性能指标,本文将以CentOS为例,介绍如何通过命令行快速获取这些关键信息,... 目录引言一、查询CPU核心数编程(几C?)1. 使用 nproc(最简单)2. 使用 lscpu(详细信

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

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

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

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