zoukankan
html css js c++ java
HDU1588 Gauss Fibonacci 矩阵应用
/* 由f(n)=f(n-1)+f(n-2)构造矩阵 令 Fn(1,2) = |f(n) f(n-1)| A(2 2) = |1 1| |1 0| F1(1 2) = |1 0| Fn = F(n-1) * A = F1 * A ^ (n-1) 代入n = g(i) = k * i + b 得:Fn = F1 * A ^ (k * i + b - 1) 1)b > 1 Fn = F1 * A^(b-1) * (A^k)^i Sum(Fn) = F1 * A^(b-1) * Sum((A^k)^i) 2)b = 1 Fn = F1 * (A^k)^i Sum(Fn) = F1 * Sum((A^k)^i) 3)b = 0 Fn = F1 * A^(k-1) * (A^k)^i------------------n = n - 1 Sum(Fn) = F1 * A^(k-1) * Sum((A^k)^i)--------n = n - 1 其中Sum((A^k)^i)用solve进行二分求和 */ #include "Mat.h" #include <iostream> using namespace std; int main() { Mat A(2, 2), F1(1, 2), e(2, 2), ans(1, 2), temp(2, 2); int b, n, k; while(cin>>k>>b>>n>>mod) { //全1矩阵 A.clear(2);A.s[1][1] = 0; F1.clear(1); e.clear(1); ans.clear(); temp = A; if(b > 1) { temp.Er_work(b-1); F1.Multiply(temp); } else if(b == 0) { n--; temp.Er_work(k-1); F1.Multiply(temp); } A.Er_work(k); A.solve(n-1); A.Add(e); F1.Multiply(A); printf("%I64d\n",F1.s[0][0]); } return 0; }
查看全文
相关阅读:
测试开发面试集锦_数据库
测试开发面试集锦_linux
测试开发面试题集锦_java
Java中equals 和==的区别
定时清理文件shell脚本
java文件上传,upload使用
python 获取错误日志,并发送邮件
c语言代码审计规范
渗透测试之nmap
渗透测试之GoogleHack
原文地址:https://www.cnblogs.com/windmissing/p/2559894.html
最新文章
Sql快速查看表结构
JS事件类型及介绍
JS中数组去重的几种方法
keep the bar green to keep the code clean——Junit详解(二)
keep the bar green to keep the code clean——Junit详解(一)
使用Word发表博客
结对项目——学生管理系统
JUnit 4 单元测试
Git ——分布式版本控制系统
取消文本文档中每一行前缀(如编号,行号等)
热门文章
Java中8种基本数据类型及取值范围
android和ios的系统特性区别
bindService和startService的区别
测试不同阶段的困扰何规划
谷歌开发者工具使用解读--新人
我们web前端常用的一些Array对象及应用
IE浏览器兼容的常见问题及解决方案
#define的用法
关于Failed to check the status of the service com.mytest.service.ItemService. No provider available for the service【已解决】
测试开发面试集锦_测试方面
Copyright © 2011-2022 走看看