Hdu 2256 Problem of Precision[矩阵快速幂 + 数学]

2024-06-07 03:38

本文主要是介绍Hdu 2256 Problem of Precision[矩阵快速幂 + 数学],希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=2256

题目的意思很简单粗暴。

简单思考一下,可以推出递推公式fi = fi-1 * a^2, a = sqrt(2) + sqrt(3)。

递推公式出来了,直接快速幂就好。。。但是, 我们发现10^9是很消耗精度, 可以肯定是很多的正确的做法,都是因为精度的问题而挂的。。。

思考无果后, 看了题解。。

表示还需要努力学习。。!!!

表示很吊。。。。构造5 - 2*sqrt(6),是为了避免sqrt(6)所带来的精度。。

Code:

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;const int N = 4;
const int mod = 1023;
struct Matrix
{int n, m;int a[N][N];Matrix(){memset(a, 0, sizeof(a));}
}unit, A;Matrix operator *(Matrix a, Matrix b)
{Matrix ans;ans.n = a.n; ans.m = b.m;for(int i = 1; i <= a.n; i ++){for(int j = 1; j <= b.m; j ++){for(int k = 1; k <= b.n; k ++)ans.a[i][j] = (ans.a[i][j] + a.a[i][k] * b.a[k][j]) & mod;}}return ans ;
}int power(int k)
{Matrix ans = unit, B = A;while(k){if(k & 1) ans = ans * B;B = B * B;k = k >> 1;}return 2 * ans.a[1][1] - 1;
}int main()
{int T;scanf("%d", &T);while(T --){int n;scanf("%d", &n);unit.n = 1; unit.m = 2;unit.a[1][1] = 5; unit.a[1][2] = 2;A.n = 2; A.m = 2;A.a[1][1] = 5; A.a[1][2] = 2; A.a[2][1] = 12; A.a[2][2] = 5;printf("%d\n", power(n - 1) & mod);}return 0;
}


期间,见到了一个 &1023的位运算。。

其实 x & 1023 = x % 1024。。这样的形式是第一次见,表示很好很强大。。

这种形式只使用于 mod = 2 ^ n - 1。为什么??我想不用解释。。。

这篇关于Hdu 2256 Problem of Precision[矩阵快速幂 + 数学]的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



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

相关文章

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快速接入OpenAI大模型的方法(JDK8)

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

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

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

使用C#代码计算数学表达式实例

《使用C#代码计算数学表达式实例》这段文字主要讲述了如何使用C#语言来计算数学表达式,该程序通过使用Dictionary保存变量,定义了运算符优先级,并实现了EvaluateExpression方法来... 目录C#代码计算数学表达式该方法很长,因此我将分段描述下面的代码片段显示了下一步以下代码显示该方法如

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

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

Rust中的Option枚举快速入门教程

《Rust中的Option枚举快速入门教程》Rust中的Option枚举用于表示可能不存在的值,提供了多种方法来处理这些值,避免了空指针异常,文章介绍了Option的定义、常见方法、使用场景以及注意事... 目录引言Option介绍Option的常见方法Option使用场景场景一:函数返回可能不存在的值场景