技术标签: DP
传送门 POJ 3070
裸的矩阵快速幂,若题中没有给出矩阵形式,可以根据递推式写出。
一般的,对于 m 项递推式,如果记递推式为
a n + m = ∑ i = 0 m − 1 b i a n + i a_{n+m}=\sum\limits_{i=0}^{m-1}b_{i} a_{n+i} an+m=i=0∑m−1bian+i
则可以把递推式写成如下矩阵形式
( a n + m a n + m − 1 ⋮ a n + 1 ) = ( b m − 1 ⋯ b 1 b 0 1 ⋯ 0 0 ⋮ ⋱ ⋮ ⋮ 0 ⋯ 1 0 ) ( a n + m − 1 a n + m − 2 ⋮ a n ) \begin{pmatrix} a_{n+m} \\ a_{n+m-1}\\ \vdots \\ a_{n+1} \\ \end{pmatrix} = \begin{pmatrix} b_{m-1} & \cdots & b_{1} & b_{0} \\ 1 & \cdots &0 & 0 \\ \vdots & \ddots & \vdots & \vdots \\ 0 & \cdots & 1 & 0 \\ \end{pmatrix}\begin{pmatrix} a_{n+m-1} \\ a_{n+m-2} \\ \vdots \\ a_{n} \\ \end{pmatrix} ⎝⎜⎜⎜⎛an+man+m−1⋮an+1⎠⎟⎟⎟⎞=⎝⎜⎜⎜⎛bm−11⋮0⋯⋯⋱⋯b10⋮1b00⋮0⎠⎟⎟⎟⎞⎝⎜⎜⎜⎛an+m−1an+m−2⋮an⎠⎟⎟⎟⎞
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <vector>
using namespace std;
#define inf 0x3f3f3f3f
#define mod 10000
typedef vector<int> vec;
typedef vector<vec> mat;
mat mul(mat &A, mat &B)
{
mat C(A.size(), vec(B[0].size()));
for (int i = 0; i < A.size(); i++)
{
for (int k = 0; k < B.size(); k++)
{
for (int j = 0; j < B[0].size(); j++)
{
C[i][j] = (C[i][j] + A[i][k] * B[k][j]) % mod;
}
}
}
return C;
}
mat pow(mat A, int n)
{
mat B(A.size(), vec(A.size()));
for (int i = 0; i < A.size(); i++) B[i][i] = 1;
while (n)
{
if (n & 1) B = mul(B, A);
A = mul(A, A);
n >>= 1;
}
return B;
}
int main()
{
int n;
while (~scanf("%d", &n) && n != -1)
{
mat A(2, vec(2));
A[0][0] = 1, A[0][1] = 1;
A[1][0] = 1, A[1][1] = 0;
A = pow(A, n);
printf("%d\n", A[0][1]);
}
return 0;
}
文章浏览阅读547次,点赞21次,收藏27次。据统计,中国软件外包市场的潜力和机会已远远超过软件王国印度,不过由于软件人才的严重不足致使我国软件发展遭遇“瓶颈”。国家为了大力培养软件人才,不断采取积极有效的措施。我国对软件测试人才的需求数量还将持续增加,因此软件测试工程师也就成为了IT职业的亮点。那么一般究竟需要哪些主流的软件测试工具呢?_软件测试所需软件
文章浏览阅读134次。一、总结与计划1.1 总结 SpringBoot自动配置回顾 基于SpringBoot整合SSM开发 一个SpringBoot项目本身就是依赖Spring:spring-boot-starter 整合SpringMVC:spring-boot-starter-web 整合MyBatis Druid:druid-spring-boot-starter MyBatis:mybatis-spring-boot-sta_thymeleafautoconfig
文章浏览阅读340次,点赞6次,收藏6次。.NET framework 不能正常安装解决办法_net framework 3.5无法勾选
文章浏览阅读2.2k次。在学习springcloud的时候跟着老师创建项目,在maven中导入springcloud依赖和springcloudAlibaba依赖的时候,maven坐标一直爆红,刚开始以为是springcloud和springcloudAlibaba的版本有问题,这里我说一下我jdk使用的是11,maven是3.3.9.springboot是2.7.15,springcloud和springcloudAlibaba使用的是2021.0.5,我换了很多个版本都没有用。编译之后控制台报错为。_com.alibaba.cloud:spring-cloud-alibaba-dependencies:jar:2021.0.5.0 was not f
文章浏览阅读6.5w次,点赞14次,收藏103次。ping 的原理 ping 程序是用来探测主机到主机之间是否可通信,如果不能ping到某台主机,表明不能和这台主机建立连接。ping 使用的是ICMP协议,它发送icmp回送请求消息给目的主机。ICMP协议规定:目的主机必须返回ICMP回送应答消息给源主机。如果源主机在一定时间内收到应答,则认为主机可达。 ICMP协议通过IP协议发送的,IP协议是一种无连接的,不可_ping192.168.1.1mac
文章浏览阅读1.1k次。可以用while(n /= 10)的方法。_c++求位数
文章浏览阅读2.7k次,点赞6次,收藏4次。最近在项目开发和学习的过程中,涉及到了Vant UI(简称Vant)的使用,主要还是涉及到了Card卡片和Cell单元格的使用,这个系列的文章就简单记录一下这个实战过程和学习经历,这篇文章先介绍 Vant UI入门与使用。_vant ui
文章浏览阅读1.3k次,点赞20次,收藏22次。php+mysql项目如何运行_phpstudy怎么运行项目
文章浏览阅读1.2k次。参考链接: Android NFC 开发实例 NfcA/NfcB/NfcF/NfcV/IsoDep/Ndef/Mifare/Felica/Pboc/ISOxxxx 都是些什么鸟玩意?android NFC学习笔记(一) NFC课件1.NFC基础知识具体测试过程一、NFC基础知识1、NFC是什么NFC,即Near Field Communication,近距离无线通讯技术,是一种短距离的(通常<_android nfc +dep
文章浏览阅读783次。因double有效位数为16位这就会出现存储小数位数不够的情况,小数的二进制有时并不能精确表示。bigDecimal是java.math包中提供的API类,用来对超过16位有效位的数进行精确的运算。双精度浮点数使用8个字节(64位)存储空间来存储一个浮点数,包括符号位1位,阶码11位,尾数52位。单精度浮点数占用4个字节(32位)存储空间来存储一个浮点数,包括符号位1位,阶码8位,尾数23位。BigDecimal(double)创建一个具有参数所指定双精度值的对象。..._java金额计算float和bigdecimal对比
文章浏览阅读1.8k次,点赞125次,收藏107次。《嵌入式虚拟化技术与应用》:深入浅出阐述嵌入式虚拟机原理,实现“小而能”嵌入式虚拟机!
文章浏览阅读114次。【代码】02.fgetc()和fputc()---->cat。_通过fgetc与fputc函数实现cat功能