fibonacci构造矩阵-总结

最新推荐文章于 2020-03-31 21:51:51 发布
最新推荐文章于 2020-03-31 21:51:51 发布
阅读量1.3k 收藏
点赞数
分类专栏： C/C++ 数据结构与算法 文章标签： c++ 算法

原文地址：

http://wenku.baidu.com/link?url=3PbnE95-IX6Sm26SQnkFkOb8xnqzzbcAqkAaPpiPZE-OXQGu3Aa6yrBR5Ot1a0a09-QYPQJ3i6maTlrQ0IX__j2AMA-tIUZ0QKw3LEYCTxi



Anotherkind of Fibonacci

Time Limit: 3000/1000 MS(Java/Others)    Memory Limit: 65536/65536 K (Java/Others)
Total Submission(s): 802    Accepted Submission(s): 326


Problem Description

As we all known , the Fibonacci series : F(0)= 1, F(1) = 1, F(N) = F(N - 1) + F(N - 2) (N >= 2).Now we define anotherkind of Fibonacci : A(0) = 1 , A(1) = 1 , A(N) = X * A(N - 1) + Y * A(N - 2) (N>= 2).And we want to Calculate S(N) , S(N) = A(0)2 +A(1)2+……+A(n)2.

 

 

Input

Thereare several test cases.
Each test case will contain three integers , N, X , Y .
N : 2<= N <= 231 – 1
X : 2<= X <= 231– 1
Y : 2<= Y <= 231 – 1

 

 

Output

Foreach test case , output the answer of S(n).If the answer is too big , divide itby 10007 and give me the reminder.

 

 

Sample Input

2 1 1

3 2 3

 

 

Sample Output

6

196

 

 

 

下为解题报告：林大陈宇

 

该题为矩阵连乘问题，主要考虑以下公式：

 

f(n)=x*f(n-1)+y*f(n-2)

 

 

先看 s(n)=s(n-1)+f(n)^2;

右边出现了n，不可以的，只能出现n-1  或n-2

所以：s(n-1)=s(n-2)+f(n-1)^2;  右边有2项，s(n-2)和f(n-1)^2必须要有了

  S(n-2)-----------às(n-1)

   F(n-1)^2------à f(n)^2=[xf(n-1)+yf(n-2)]^2=出现了f(n-2)^2和2xyf(n-1)f(n-2),所以共4项

 

 

 

最左边是1个空的矩阵(4*4)，也是我们要构造的矩阵（只要学过矩阵的人，就能把它算出来了吧）

 

然后代入矩阵连乘幂的模板算下就行了（该模板要会背着敲）

 

#include <iostream>

#include <stdio.h>

#include <cstring>

#define Mod 10007

using namespace std;

const int MAX = 4;

 

typedef  struct{

         int m[MAX][MAX];

}  Matrix;

MatrixP={1,1,0,0,

          0,0,0,0,

          0,1,0,0,

          0,0,0,0};

MatrixI={1,0,0,0,

          0,1,0,0,

          0,0,1,0,

          0,0,0,1};

Matrixmatrixmul(Matrix a,Matrix b) //矩阵乘法

{

       int i,j,k;

       Matrix c;

       for (i = 0 ; i < MAX; i++)

           for (j = 0;j < MAX;j++)

            {

                 c.m[i][j] = 0;

                 for (k = 0; k < MAX; k++)

                     c.m[i][j] += (a.m[i][k]*b.m[k][j])%Mod;

                 c.m[i][j] %= Mod;

            }

      return c;

}

 

Matrix quickpow(long long n)

{

      Matrix m = P, b = I;

              while (n >= 1)

       {

             if (n & 1)

                b = matrixmul(b,m);

             n = n >> 1;

             m =matrixmul(m,m);

      }

      return b;

}

 

 

int main()

{

    int n,x,y,sum;

    Matrix b;

/    while(cin>>n>>x>>y)

    {

      sum=0;

      x=x%Mod;

      y=y%Mod;

      P.m[1][1]=(x*x)%Mod;P.m[1][2]=(y*y)%Mod;P.m[1][3]=(2*x*y)%Mod;

     P.m[3][1]=x;P.m[3][3]=y;

       b=quickpow(n);

       for(int i=0;i<4;i++)

       sum+=b.m[0][i]%Mod;

       cout<<sum%Mod<<endl;

 

    }

        return 0;

}

 

注意： 注意每一步都要取模，别怕麻烦；

       对于x和y要是取负数的话（本题没有），则要(data%mod+mod)%mod;

       正的就直接data%mod就行；

       typedef struct{

         int m[MAX][MAX];

} Matrix;

该结构体的定义要注意，若将int  类型换成long long ，则运行时间会翻倍，可能会TLE,注意哦

 

----------------------------------------------------------------------------------------------------------------------

有兴趣的同学可以接着看，很有收获的，下面的方法和我的类似，我是把够造矩阵A放在左边，它是放在中间的，请用我的方法推1下，我给答案了。

 

 

构造常系数矩阵！

（一）Fibonacci数列f[n]=f[n-1]+f[n-2],f[1]=f[2]=1的第n项的快速求法（不考虑高精度）.

解法：

考虑1×2的矩阵【f[n-2],f[n-1]】。根据fibonacci数列的递推关系，我们希望通过乘以一个2×2的矩阵，得到矩阵【f[n-1],f[n]】=【f[n-1],f[n-1]+f[n-2]】

很容易构造出这个2×2矩阵A，即：

０１ 
１１

所以，有【f[1],f[2]】×A＝【f[2],f[3]】

又因为矩阵乘法满足结合律，故有：

【f[1],f[2]】×A n-1=【f[n],f[n+1]】

这个矩阵的第一个元素即为所求。

至于如何快速求出A n-1，相信大家都会，即递归地：n为偶数时，An=(A n/2)2；n为奇数时，An=(A n/2)2*A。

问题（一）解决。

 

我的方法：    （a）

还是填空：

     还没会吗？ (a)式中左边的空矩阵就是让你填的A

 

（二）数列f[n]=f[n-1]+f[n-2]+1,f[1]=f[2]=1的第n项的快速求法（不考虑高精度）.

解法：

仿照前例，考虑1×3的矩阵【f[n-2],f[n-1],1】，希望求得某3×3的矩阵A，使得此1×3的矩阵乘以A得到矩阵：【f[n-1],f[n],1】＝【f[n-1],f[n-1]+f[n-2]+1,1】

容易构造出这个3×3的矩阵A，即：

０１０ 
１１０ 
０１１

问题（二）解决。

 

我的推法：

 

 

A是3*3的矩阵，填空吧

 

 

（三）数列f[n]=f[n-1]+f[n-2]+n+1,f[1]=f[2]=1的第n项的快速求法（不考虑高精度）.

解法：

仿照前例，考虑1×4的矩阵【f[n-2],f[n-1],n,1】，希望求得某4×4的矩阵A，使得此1×4的矩阵乘以A得到矩阵：

【f[n-1],f[n],n+1,1】＝【f[n-1],f[n-1]+f[n-2]+n+1,n+1,1】

容易构造出这个4×4的矩阵A，即：

０１００ 
１１００ 
０１１０ 
０１１１

问题（三）解决……

 

我自己的推法

 

 

接着填空：

 

 

 

四）数列f[n]=f[n-1]+f[n-2],f[1]=f[2]=1的前n项和s[n]的快速求法（不考虑高精度）.

解法：

虽然我们有S[n]=F[n+2]-1，但本文不考虑此方法，我们想要得到更一般的方法。

考虑（一）的矩阵A，容易发现我们要求【f[1],f[2]】×（A+A2+A3+…+AN-1）。很多人使用一种很数学的方法构造一个2r*2r（r是A的阶数，这里为2）的矩阵来计算，这种方法比较麻烦且很慢，这里不再介绍。下面考虑一种新方法。

仿照之前的思路，考虑1×3的矩阵【f[n-2],f[n-1],s[n-2]】，我们希望通过乘以一个3×3的矩阵A，得到1×3的矩阵：

【f[n-1],f[n],s[n-1]】=【f[n-1],f[n-1]+f[n-2],s[n-2]+f[n-1]】

容易得到这个3×3的矩阵是：

０１０ 
１１１ 
００１

然后…………容易发现，这种方法的矩阵规模是(r+1)*(r+1)，比之前流行的方法好得多。

我用分块矩阵做过这题，但该方法的确应该推广，是主流的思想。

我的推法：

，所以矩阵见下，填空

为所求

 

 

 

（五）数列f[n]=f[n-1]+f[n-2]+n+1,f[1]=f[2]=1的前n项和s[n]的快速求法（不考虑高精度）.

解法：

结合（三）（四），容易想到……

考虑1×5的矩阵【f[n-2],f[n-1],s[n-2],n,1】,

我们需要找到一个5×5的矩阵A，使得它乘以A得到如下1×5的矩阵：

【f[n-1],f[n],s[n-1],n+1,1】

=【f[n-1], f[n-1]+f[n-2]+n+1,s[n-2]+f[n-1],n+1,1】

容易构造出A为：

０１０００ 
１１１００ 
００１００ 
０１０１０ 
０１０１１

然后……问题解决。

 

我的推法：

 

，填空

 

为构造的矩阵。

 

 

陈宇于2012年7月15日


duhaomin
关注
专栏目录
基于Fibonacci-PEG算法的LDPC码校验矩阵的构造
qq_37934722的博客
09-12 227
在LDPC码的构造中，校验矩阵起着至关重要的作用，它决定了编码和解码的性能。本文将介绍如何使用Matlab实现基于Fibonacci-PEG算法的LDPC码校验矩阵的构造，并提供相应的源代码。通过以上的代码和说明，我们可以使用Matlab实现基于Fibonacci-PEG算法的LDPC码校验矩阵的构造。校验矩阵的生成过程可以根据需要进行调整，例如调整校验矩阵的大小和停止条件，以满足具体应用需求。在上述示例中，我们设置校验矩阵的列数为100，行数为50，并指定最大迭代次数为200。最后，我们打印校验矩阵。
矩阵与斐波那契数列
SalamAbdu的博客
05-25 1269
矩阵快速幂 斐波那契数列 斐波那契数列（Fibonacci sequence），又称黄金分割数列。0、1、1、2、3、5、8、13、21、34、……在数学上，斐波那契数列以如下被以递推的方法定义： F(n)={n,0≤n≤1F(n−1)+F(n−2)（n≥2，n∈N∗）F(n) = \begin{cases} n , \quad 0 \le n \le 1 \\ F(n - 1)+F(n - 2)（n ≥ 2，n ∈ N*） \end{cases}F(n)={n,0≤n≤1F(n−1)+F(n−2)（n≥
Fibonacci（矩阵）
weixin_30514745的博客
03-13 154
Fibonacci Time Limit:1000MSMemory Limit:65536KB64bit IO Format:%I64d & %I64u Description In the Fibonacci integer sequence,F0= 0,F1= 1, andFn=Fn− 1+Fn− 2forn...
又见斐波那契(构造矩阵，矩阵快速幂)
obsorb_knowledge的博客
05-19 967
链接：https://www.nowcoder.com/acm/contest/105/G来源：牛客网题目描述这是一个加强版的斐波那契数列。 给定递推式 求F(n)的值，由于这个值可能太大，请对109+7取模。输入描述:第一行是一个整数T(1 ≤ T ≤ 1000)，表示样例的个数。 以后每个样例一行，是一个整数n(1 ≤ n ≤ 1018)。 输出描述:每个样例输出一行，一个整数，表示F(...
Fibonacci +矩阵
地衣的博客——qdu
06-09 388
#include<cstdio> #include<cstring> using namespace std; struct matrix { long long m[2][2]; }; const long long mod=10000; matrix mul(matrix a,matrix b) { matrix temp; memset(temp.m,0,siz...
Fibonacci --矩阵快速幂
jialileng的博客
08-01 172
菲波那契数列是指这样的数列: 数列的第一个是0和第二个数是1，接下来每个数都等于前面2个数之和。 给出一个正整数a，要求菲波那契数列中第a个数的后四位是多少。 Input 多组数据 -1结束 范围1~10^9 Output 第x项的后4位 Sample Input 0 9 999999999 1000000000 -1 Sample Output 0 34 626 6875 ac代码： #incl...
fibonacci构造矩阵-总结.doc
09-05
（一）Fibonacci数列f[n]=f[n-1]+f[n-2],f[1]=f[2]=1的第n项的快速求法（不考虑高精度）. 解法： 考虑1×2的矩阵【f[n-2],f[n-1]】。根据fibonacci数列的递推关系，我们希望通过乘以一个2×2的矩阵，得到矩阵【f[n-1],f[n]】=【f[n-1],f[n-1]+f[n-2]】
利用斐波那契数列构造QC-LDPC码的方法
02-23
在编码理论中，特别是低密度奇偶校验（LDPC）码的构造中，斐波那契数列也能找到其身影。 LDPC码是一类具有稀疏校验矩阵的线性分组码，这类码最早由Gallager在1963年提出。LDPC码因为其优异的纠错能力和相对简单的...
Fibonacci矩阵
11-21 2045
我们如下定义一个无限阶矩阵m 矩阵第0行正好是Fibonacci数列 也就是 m(0,0)=1,m(0,1)=2,m(0,2)=3, m(0,3)=5,m(0,4)=8,.... 矩阵第k行中第一个数字是前面k-1行中都没有出现的最小正整数, 所以 m(1,0)=4, 而m(k,1)=2*m(k,0)-k (这个关系对于第0行也成立) 所以m(1,1)=2*4-1=7 而第k行后面的任意一个数同第...
矩阵-斐波那契数列
热门推荐
许增强
06-08 1万+
利用矩阵来求解斐波那契数列的有关问题是ACM题中一个比较常见的题型。例：NYOJ 148(斐波那契数列2)。 有关斐波那契树列的规律详见这里。 (1)、对于n>1，都有f(n)与f(n-1)互质。 (2)、f(n)=f(i)*f(n-i-1)+f(i+1)*f(n-i)。 现在说说怎么利用矩阵来求解斐波那契数列。 我们可以先保存b=f(1),a=f(0)，然后每次设： b'=a+b a'
Fibonacci(矩阵快速幂 - 递推式构造 + 变种)
Sun
10-04 279
文章目录POJ 3070 Fibonacci2018年湘潭大学程序设计竞赛 G 又见斐波那契 POJ 3070 Fibonacci 题目链接：http://poj.org/problem?id=3070 题意：求Fn % 1000是多少？(公式题目已给出) 此题的公式就是Fn = Fn-1 + Fn-2，我们令:   Xn−1=[Fn−1Fn−2]X_{n-1} = \begin{bmatrix...
Fibonacci（矩阵快速幂）
oneplus123的博客
12-06 265
Fibonacci 时间限制: 1 Sec  内存限制: 128 MB 题目描述 In the Fibonacci integer sequence, F0 = 0, F1 = 1, and Fn = Fn − 1 + Fn − 2 for n ≥ 2. For example, the first ten terms of the Fibonacci sequence are: 0, 1,...
【牛客】斐波那契数列卷积（矩阵快速幂）（构造矩阵）
Nerdant
10-27 354
链接：https://ac.nowcoder.com/acm/contest/1087/C 题目描述 输入描述: 输出描述: 一行一个整数表示答案 解释： 我没做出来的原因是没推出来这个式子，矩阵的题目就是前后有联系，这样的式子我没推出来。 按照斐波那契的性质推了很久也没有发现什么公式，看题解也是直接说“很容易就得到”....不过我们可以通过直接对所求式子打...
【poj3070】Fibonacci（矩阵）
蜗牛
08-21 433
In the Fibonacci integer sequence, F0 = 0, F1 = 1, and Fn = Fn − 1 + Fn − 2 for n ≥ 2. For example, the first ten terms of the Fibonacci sequence are:0, 1, 1, 2, 3, 5, 8, 13, 21, 34, …An alternative fo
【全国互虐】Fibonacci矩阵
dee2784的博客
05-28 156
orz啊又被屠了 人生如此艰难 题意： 给定一个k维的n^k的超立方体 超立方体的元素Ai1,i2,...,ik 的值为f(i1+i2+...+ik-k+1) f为斐波那契数列 求该超立方体的所有元素和 1<=n,k<=10^9 题解： 其实看到数据范围 就大概猜到是矩阵乘法了 但是我考试的时候想了半天还是不知道矩阵乘法怎么搞 - - 其实矩阵乘法...
Fibonacci矩阵的探讨与应用
11-28
很不错的资料，写的比较详细，主要介绍的数列用矩阵实现的方法、原理等等
Fibonacci数列矩阵法
随笔
11-21 428
/*斐波那契数列矩阵解法*/ /* T(n)=T(n-1)+T(n-2) T(n) = T(n-1)+T(n-2) =1 1 T(n-1) = 1 1 1 1 T(n-2) T(n-1) = T(n-1) =1 0 T(n-2) = 1 0 1 0 T(n-3) */ #include<iostre...
fibonacci 利用矩阵乘法求解
fareast_mzh的博客
08-01 668
Fibonacci.hpp 假设一个整数是10，如何最快的求10的75次方 1. 75 =&gt; Binary 1001 1011 2. 10^75 = 10^64 + 10^8 + 10^2 + 10^1 10^1 =&gt; 10^2 =&gt; 10^4 =&gt; ... =&gt; 10^32 =&gt; 10^64 #pragma once #ifndef Fibona...
Fibonacci (矩阵快速幂)
王小平
10-02 581
题目链接 Description In the Fibonacci integer sequence, F0 = 0, F1 = 1, and Fn = Fn − 1 + Fn − 2 for n ≥ 2. For example, the first ten terms of the Fibonacci sequence are: 0, 1, 1, 2, 3, 5, 8, 13, 21, 34,...
矩阵K次幂（fibonacci）
ss123mlk的博客
03-31 6549
@[TOC](Fibanacci 的复杂度O(logn)方式) 快速幂 [快速幂运算方式] (https://blog.csdn.net/flyfish1986/article/details/48014523) 根据上面的链接 在求幂的时候 暴力算法是不停的乘本身A 需要很多次的调用自己本身 乘很多次 时间复杂度高 但是所采用“分治法”将 1.多次幂运算乘以自己这个大问题分成若干的小问题 2.然...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值