hdu 4565 So Easy!(矩阵快速幂)

2024-06-05 03:48
文章标签 快速 矩阵 hdu easy 4565

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

题目链接:hdu 4565 So Easy!


题目大意:计算题目中所给式子的值。


解题思路:矩阵快速幂。

An = (a+√b)^n, Bn = (a-√b)^n, Cn = An + Bn = (a+√b)^n+(a-√b)^n;

因为说An和Bn共轭,展开式可以互相抵消,所以保证说Cn一定是整数。


(a-1)^2 < b < a^2

-> a-1 < √b < a

-> 0 < a - √b < 1

-> 0 < (a - √b) < 1

-> [ Bn ] = 1

所以有Cn = [ An ]

Sn = Cn % m.


Sn * 2a = Sn * ((a + √b) + (a - √b)) = ((a + √b)^n + (a - √b)^n)*((a + √b) + (a - √b)) 

= ((a + √b)^(n+1) + (a - √b)^(n+1)) + (a^2-b)((a + √b)^(n-1) + (a - √b)^(n-1))

= Sn+1 + (a^2-b)Sn-1


Sn+1 = 2aSn + (b-a^2)Sn-1.


#include <stdio.h>
#include <string.h>
#include <math.h>#define MOD(a,b) (((a%b)+b)%b)
typedef long long ll;
const int N = 2;int m;
struct mat {ll s[N][N];mat () {memset(s, 0, sizeof(s));}void set (ll a, ll b, ll c, ll d) {s[0][0] = a; s[0][1] = b;s[1][0] = c; s[1][1] = d;}mat operator * (const mat& c) {mat ans;for (int i = 0; i < N; i++) {for (int j = 0; j < N; j++) {for (int k = 0; k < N; k++)ans.s[i][j] = MOD(ans.s[i][j] + s[i][k] * c.s[k][j], m);}}return ans;}
}A;mat powMat (int n) {mat ans;ans.set (1, 0, 0, 1);while (n) {if (n&1)ans = ans * A;n /= 2;A = A * A;}return ans;
}int main () {int a, b, n;while (scanf("%d%d%d%d", &a, &b, &n, &m) == 4) {double t = pow (a + sqrt(b), 2);ll s1 = 2 * a % m;ll s2 = (ll)ceil(t) % m;if (n <= 1) {printf("%lld\n", s1);} else {A.set (2 * a, b - a * a, 1, 0);mat ans = powMat (n-2);printf("%lld\n", MOD(s2 * ans.s[0][0] + s1 * ans.s[0][1], m));}}return 0;
}

    



这篇关于hdu 4565 So Easy!(矩阵快速幂)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



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

相关文章

springboot security快速使用示例详解

《springbootsecurity快速使用示例详解》:本文主要介绍springbootsecurity快速使用示例,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝... 目录创www.chinasem.cn建spring boot项目生成脚手架配置依赖接口示例代码项目结构启用s

C++快速排序超详细讲解

《C++快速排序超详细讲解》快速排序是一种高效的排序算法,通过分治法将数组划分为两部分,递归排序,直到整个数组有序,通过代码解析和示例,详细解释了快速排序的工作原理和实现过程,需要的朋友可以参考下... 目录一、快速排序原理二、快速排序标准代码三、代码解析四、使用while循环的快速排序1.代码代码1.由快

Win32下C++实现快速获取硬盘分区信息

《Win32下C++实现快速获取硬盘分区信息》这篇文章主要为大家详细介绍了Win32下C++如何实现快速获取硬盘分区信息,文中的示例代码讲解详细,感兴趣的小伙伴可以跟随小编一起学习一下... 实现代码CDiskDriveUtils.h#pragma once #include <wtypesbase

Spring AI与DeepSeek实战一之快速打造智能对话应用

《SpringAI与DeepSeek实战一之快速打造智能对话应用》本文详细介绍了如何通过SpringAI框架集成DeepSeek大模型,实现普通对话和流式对话功能,步骤包括申请API-KEY、项目搭... 目录一、概述二、申请DeepSeek的API-KEY三、项目搭建3.1. 开发环境要求3.2. mav

Python如何快速下载依赖

《Python如何快速下载依赖》本文介绍了四种在Python中快速下载依赖的方法,包括使用国内镜像源、开启pip并发下载功能、使用pipreqs批量下载项目依赖以及使用conda管理依赖,通过这些方法... 目录python快速下载依赖1. 使用国内镜像源临时使用镜像源永久配置镜像源2. 使用 pip 的并

SpringBoot集成图片验证码框架easy-captcha的详细过程

《SpringBoot集成图片验证码框架easy-captcha的详细过程》本文介绍了如何将Easy-Captcha框架集成到SpringBoot项目中,实现图片验证码功能,Easy-Captcha是... 目录SpringBoot集成图片验证码框架easy-captcha一、引言二、依赖三、代码1. Ea

SpringBoot快速接入OpenAI大模型的方法(JDK8)

《SpringBoot快速接入OpenAI大模型的方法(JDK8)》本文介绍了如何使用AI4J快速接入OpenAI大模型,并展示了如何实现流式与非流式的输出,以及对函数调用的使用,AI4J支持JDK8... 目录使用AI4J快速接入OpenAI大模型介绍AI4J-github快速使用创建SpringBoot

SpringBoot整合easy-es的详细过程

《SpringBoot整合easy-es的详细过程》本文介绍了EasyES,一个基于Elasticsearch的ORM框架,旨在简化开发流程并提高效率,EasyES支持SpringBoot框架,并提供... 目录一、easy-es简介二、实现基于Spring Boot框架的应用程序代码1.添加相关依赖2.添

使用Python快速实现链接转word文档

《使用Python快速实现链接转word文档》这篇文章主要为大家详细介绍了如何使用Python快速实现链接转word文档功能,文中的示例代码讲解详细,感兴趣的小伙伴可以跟随小编一起学习一下... 演示代码展示from newspaper import Articlefrom docx import

shell脚本快速检查192.168.1网段ip是否在用的方法

《shell脚本快速检查192.168.1网段ip是否在用的方法》该Shell脚本通过并发ping命令检查192.168.1网段中哪些IP地址正在使用,脚本定义了网络段、超时时间和并行扫描数量,并使用... 目录脚本:检查 192.168.1 网段 IP 是否在用脚本说明使用方法示例输出优化建议总结检查 1