一类常见的矩阵快速幂题型 总结

最新推荐文章于 2021-11-01 20:34:22 发布
最新推荐文章于 2021-11-01 20:34:22 发布
阅读量2k 收藏 5
点赞数 1
分类专栏： =====ACM===== 【数学/大数】 【矩阵计算/快速幂】 ACM算法编程 文章标签： 数论 矩阵快速幂
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/u013050857/article/details/48637557
版权

 

Googlecode jam 2008 Round 1A(c.Numbers)

【题意】计算的小数点前三位数，不足三位补0,正整数n的最大值为20亿。

【前提】：满足


的值在【0，1】范围

首先将


展开之后可以发现


的形式，同样的，有


因此，


是个整数，

其中



这正是解题的关键！

由于


所以的整数部分等于


根据以上的推导

只要高效的求出an就可以解决这个问题了

由于


为观察仔细，进一步展开得：


得出


的递推关系


因此，可以用矩阵表示这个递推关系，使用矩阵快速幂，在O（logn）的时间里求出和,由于此题只要取前3位，因此对1000取模就可以了

代码：

#include <stdio.h>
#include <string.h>
#include <iostream>
#include <algorithm>

using namespace std;
typedef long long LL;
const LL siz=2;          // max size of the matrix,
#define MODD(a,b) (((a%b)+b)%b)

LL A,B,N,M,ret;
struct mut{
    LL mat[siz][siz];
    mut(){
        memset(mat,0,sizeof(mat));
    }
    void init(LL a,LL b,LL c,LL d){
        mat[0][0]=a;
        mat[0][1]=b;
        mat[1][0]=c;
        mat[1][1]=d;
    }
    mut operator *(const mut &c){
        mut res;
        for(int i=0; i<siz; ++i){
            for(int k=0; k<siz; ++k){
                for(int j=0; j<siz; ++j){
                    res.mat[i][j]=MODD(res.mat[i][j]+mat[i][k]*c.mat[k][j],M);
                }
            }
        }
        return res;
    }
} AC;

mut poww(LL n){
    mut ans;
    ans.init(1,0,0,1);
    while(n){
        if(n&1) ans=ans*AC;
        n>>=1;
        AC=AC*AC;
    }
    return ans;
}
int main()
{
    int t;
    scanf("%d",&t);
    while(t--){
        scanf("%lld",&N);
        M=1000;
        AC.init(3,5,1,3);
        mut ans=poww(N);
        printf("%03lld\n",MODD(2*ans.mat[0][0]-1,M));
    }
    return 0;
}

HDU 2256    Problem of Precision

【题意】：


【思路】和上题思路一样：


代码：

/*  
* Problem: HDU No.2256
* Running time: 15MS  
* Complier: G++  
* Author: javaherongwei 
* Create Time: 23:56 2015/9/21 星期一
*/  

#include <stdio.h>
#include <string.h>
#include <iostream>
#include <algorithm>

using namespace std;
typedef long long LL;
const LL siz=2;          // max size of the matrix,
#define MODD(a,b) (((a%b)+b)%b)

LL A,B,N,M,ret;
struct mut{
    LL mat[siz][siz];
    mut(){
        memset(mat,0,sizeof(mat));
    }
    void init(LL a,LL b,LL c,LL d){
        mat[0][0]=a;
        mat[0][1]=b;
        mat[1][0]=c;
        mat[1][1]=d;
    }
    mut operator *(const mut &c){
        mut res;
        for(int i=0; i<siz; ++i){
            for(int k=0; k<siz; ++k){
              for(int j=0; j<siz; ++j){
                    res.mat[i][j]=MODD(res.mat[i][j]+mat[i][k]*c.mat[k][j],M);
                }
            }
        }
        return res;
    }
} AC;

mut poww(LL n){
    mut ans;
    ans.init(1,0,0,1);
    while(n){
        if(n&1) ans=ans*AC;
        n>>=1;
        AC=AC*AC;
    }
    return ans;
}
int main(){
    int t;
    scanf("%d",&t);
    while(t--){
        scanf("%lld",&N);
        M=1024;
        AC.init(5,12,2,5);
        mut ans=poww(N);
        printf("%lld\n",MODD(2*ans.mat[0][0]-1,M));
    }
    return 0;
}

HDU 1575 Tr A 

【题目链接】：click here~~

24K纯裸的矩阵计算：

【题意】：给定一个n*n的矩阵要求矩阵的k次幂之后的矩阵的对角线的和

代码：

/*
* Problem: HDU No.1575
* Running time: 0MS
* Complier: G++
* Author: javaherongwei
* Create Time: 8:10 2015/9/22 星期二
*/

#include <stdio.h>
#include <string.h>
#include <iostream>
#include <algorithm>

using namespace std;
typedef long long LL;
#define MODD(a,b) (((a%b)+b)%b)

LL A,B,K,M,ret;
LL N;
LL siz;          // max size of the matrix,
struct mut{
    LL mat[10][10];
    mut(){
        memset(mat,0,sizeof(mat));
    }
    void init(){
      for(int i=0; i<siz; ++i)
        for(int j=0; j<siz; ++j)
          if(i==j)mat[i][j]=1;
          else mat[i][j]=0;
    }
    mut operator *(const mut &c){
        mut res;
        for(int i=0; i<siz; ++i){
            for(int k=0; k<siz; ++k){
              for(int j=0; j<siz; ++j){
                    res.mat[i][j]=MODD(res.mat[i][j]+mat[i][k]*c.mat[k][j],M);
                }
            }
        }
        return res;
    }
} AC;

mut poww(LL n){
    mut ans;
    ans.init();
    while(n){
        if(n&1) ans=ans*AC;
        n>>=1;
        AC=AC*AC;
    }
    return ans;
}
int main(){
    int t;scanf("%d",&t);
    while(t--){
        scanf("%d %lld",&N,&K);M=9973;siz=N;
        for(int i=0; i<N; ++i){
            for(int j=0; j<N; ++j)
                scanf("%d",&AC.mat[i][j]);
        }
        mut ans=poww(K);
        LL sum=0;
        for(int i=0; i<N; ++i){
            sum+=ans.mat[i][i];
        }
        printf("%lld\n",MODD(sum,M));
    }
    return 0;
}

CodeForces 185A Plant

【题目链接】

点击打开cf 185A

【思路】: 递推+矩阵快速幂

分析:

1 题目要求找到在n年后向上三角形的个数

2 写出前面的几个数F(0) = 1 , F(1) = 3 , F(2) = 10 , F(3) = 36 ，我们可以发现每一个向上三角形在后面形成的时候必然是分成了3个向上+1个向下的，总共分成了4个，这是包含了向下的，所以要减去，然后发现F(1向下)=1，F(2想下)=3-1=2，因此向下的可以推出关系：


从而找到通项公式 ：



构造矩阵：


之后就是矩阵快速幂了

代码：

/*
* Problem: CodeForces 185A
* Running time: 8MS
* Complier: G++
* Author: javaherongwei
* Create Time: 9:10 2015/9/22 星期二
*/

#include <stdio.h>
#include <string.h>
#include <iostream>
#include <algorithm>

using namespace std;
typedef long long LL;
#define MODD(a,b) (((a%b)+b)%b)

LL A,B,K,M,ret;
LL N;
const int siz=2;          // max size of the matrix,
struct mut{
    LL mat[siz][siz];
    mut(){
        memset(mat,0,sizeof(mat));
    }
    void init(LL a,LL b,LL c,LL d){
      mat[0][0]=a;mat[0][1]=b;
      mat[1][0]=c;mat[1][1]=d;
    }
    mut operator *(const mut &c){
        mut res;
        for(int i=0; i<siz; ++i){
            for(int k=0; k<siz; ++k){
              for(int j=0; j<siz; ++j){
                    res.mat[i][j]=MODD(res.mat[i][j]+mat[i][k]*c.mat[k][j],M);
                }
            }
        }
        return res;
    }
} AC;

mut poww(LL n){
    mut ans;
    ans.init(1,0,0,1);
    while(n){
        if(n&1) ans=ans*AC;
        n>>=1;
        AC=AC*AC;
    }
    return ans;
}
int main(){
    //int t;scanf("%d",&t);
   // while(t--){
   while(~scanf("%lld",&N)){
        M=1e9+7;
        AC.init(4,-1,0,2);
        mut ans=poww(N);
        printf("%lld\n",MODD(ans.mat[0][0]*1+ans.mat[0][1]*1,M));
    }
    return 0;
}

HDU 2842 Chinese Rings

【题目链接】点击打开hdu2842

【思路】: 矩阵快速幂

【分析】:

1 题目的意思是给定n个环，和一些规则：

If the first k rings are all off and the (k + 1)th ring is on, then the (k + 2)th ring can be taken off or taken on with one step. (0 ≤ k ≤ 7)

要把所有的环全部拆下最少需要的步数

2 题目规定如果要拆第n个环，那么第n-1个要挂着，n-2环要被拆下。那么我们设f(n)表示拆下前n个环的最少的步骤

 那么考虑第n个环的情况，第n-1个环必须要挂着，n-2环要拆下，那么这一步就要f(n-2)，拆下第n个需要1步。然后只剩下第n-1个环，由于n-1环需要第n-2环挂着，所以我们需要把前n-2个环挂上去，所以需要f(n-2)，剩下n-1个需要拆下需要f(n-1)。那么总的需要:




3 接下来利用矩阵快速幂即可

代码：

/*
* Problem: HDU 2842
* Running time: 0MS
* Complier: G++
* Author: javaherongwei
* Create Time: 9:10 2015/9/22 星期二
*/

#include <stdio.h>
#include <string.h>
#include <iostream>
#include <algorithm>

using namespace std;
typedef long long LL;
#define MODD(a,b) (((a%b)+b)%b)

LL A,B,K,M,ret;
LL N;
const int siz=3;          // max size of the matrix,
struct mut{
    LL mat[siz][siz];
    mut(){
        memset(mat,0,sizeof(mat));
    }
    void init(LL a1,LL a2,LL a3,LL b1,LL b2,LL b3,LL c1,LL c2,LL c3){
      mat[0][0]=a1;mat[0][1]=a2;mat[0][2]=a3;
      mat[1][0]=b1;mat[1][1]=b2;mat[1][2]=b3;
      mat[2][0]=c1;mat[2][1]=c2;mat[2][2]=c3;
    }
    mut operator *(const mut &c){
        mut res;
        for(int i=0; i<siz; ++i){
            for(int k=0; k<siz; ++k){
              for(int j=0; j<siz; ++j){
                    res.mat[i][j]=MODD(res.mat[i][j]+mat[i][k]*c.mat[k][j],M);
                }
            }
        }
        return res;
    }
} AC;

mut poww(LL n){
    mut ans;
    ans.init(1,0,0,0,0,1,0,0,1);
    while(n){
        if(n&1) ans=ans*AC;
        n>>=1;
        AC=AC*AC;
    }
    return ans;
}
int main(){
    //int t;scanf("%d",&t);
   // while(t--){
   while(~scanf("%lld",&N),N){
         if(N==1){
            puts("1");
            continue;
        }
        M=200907;
        AC.init(1,2,1,1,0,0,0,0,1);
        mut ans=poww(N-2);
        printf("%lld\n",MODD(ans.mat[0][0]*2+ans.mat[0][1]*1+ans.mat[0][2]*1,M));
    }
    return 0;
}


herongweiV
关注
专栏目录
蓝桥杯Python组最后几天冲刺———吐血总结,练题总结,很管用我学会了
那人独钓寒江雪
04-07 2241
一、重要知识要点 1、穷举法 2、枚举法 3、动态规划 4、回溯法 5、图论 6、深度优先搜索（DFS） 7、广度优先搜索（BFS ） 8、二叉树 9、递归 10、分治法、矩阵法 11、排列组合 12、素数、质数、水仙花数 13、欧几里得定理gcd 14、求最大公约数、最小公倍数 15、海伦公式（求三角形面积） 16、博弈论 17、贪心 18、二分查找法 19、hash表 20、日期计算 21、矩形快速幂 22、树形DP 23、最短路径 24、最小生成
Codeforces - 450B 矩阵快速幂
TRDD的博客
05-06 647
【题目大意】 已知f1=x,f2=y,∀i(i≥2)fi=fi−1+fi+1.f1=x,f2=y,∀i(i≥2)fi=fi−1+fi+1.f_1=x,f_2=y,\forall i(i\geq 2)f_i=f_{i-1}+f_{i+1}.现给你n,x,yn,x,yn,x,y要求fnmod(109+7)fnmod(109+7)f_nmod(10^9+7)。 【数据范围】 1≤n≤2⋅109...
矩阵快速幂专题【完结】
caobihole
08-25 429
第一题 hdu 1757 A Simple Math Problem 点击打开链接 思路：矩阵快速幂 分析: 1 最简单的矩阵快速幂的题目，直接利用矩阵求解即可 点击打开查看代码 第二题 hdu 1575Tr A 点击打开hdu 1575 思路: 矩阵快速幂 分析: 1 题目给定一个n*n的矩阵要求矩阵的k次幂之后的矩阵的对角线的和 2 矩阵快速幂的裸题 点...
矩阵快速幂的一些题目
ypeijasd的专栏
07-16 644
矩阵快速幂求斐波那契数列第n项，递推公式 进而转换为矩阵幂次 矩阵快速幂代码 //矩阵快速幂计算 a^n #include <bits/stdc++.h> using namespace std; int mod=1e9+7; //模数 long long n; struct matrix //定义结构体来表示矩阵 { long long m[3][3]; }; matrix a, an, fib; //fib是最终斐波那契矩阵, a是矩阵, an为计算完a
矩阵快速幂模板+常用变式
black_hole6的博客
11-06 3252
矩阵乘法（矩阵乘矩阵） 在向量乘向量的运算中，是将每个元素与它对应的元素相乘，求所有乘积之和 那么矩阵乘矩阵是不是就是两个同型矩阵的对应元素相乘呢？ 图样图森破 两个矩阵相乘的前提是前一个矩阵的列数等于后一个矩阵的行数 举个栗子，AAA为n∗kn*kn∗k矩阵，BBB为k∗mk*mk∗m矩阵，CCC为m∗nm*nm∗n矩阵，那么AAA可以与BBB相乘，BBB可以与CCC相乘，CCC可以与...
矩阵快速幂题集
wangtao的博客
07-10 649
矩阵快速幂 矩阵必须全部初始化,注意矩阵大小 POJ - 3070 最简单的矩阵快速幂 f[n]=f[n-1]+f[n-2] ac: #include <iostream> #include <cstdio> #include <cstring> #define ll long long using namespace std; /* (f(n-...
基础算法题——斐波那契（快速求幂、斐波那契特性、矩阵）
Zhengyangxinn的博客
02-20 538
基础算法题——斐波那契 本人为一名普通二本学校自动化专业的大二学生，对编程有着少许兴趣。致力将算法写得更加通俗易懂。快速求幂、斐波那契数列特性、矩阵。
Leetcode各种题型题目+思路+代码（共176道题）
JJJJJJames的博客
03-09 9401
1、Leetcode-1047 题目描述： 题目思路： 使用stack的思路，先把字符串的第一个字母压栈，然后判断下一个待压栈的字母是否与stack[-1]的字母相同，如果相同，该字母不进入stack中，且stack.pop()将重复的字母剔除。如此遍历一遍字符串，最后返回stack中的元素即为消除重复字后的字符串。 举例：“abbaca” 建立栈 把a压栈 b与a不相同，压栈 第二个b和stack栈口元素相同，不要第二个b且把栈扣的b pop掉。 第四个元素：a，此时stack中还有一个元素a，相同，
一段话总结一道题：剑指offer75道题
qq_44793993的博客
08-07 1298
文章目录03：数组中重复的数字04：二维数组中的查找05:替换空格06：从头到尾打印链表07：重建二叉树（中序+前序）09:双栈实现队列10:斐波那契、青蛙跳台阶11：旋转数组的最小数字 剑指的地址 03：数组中重复的数字 前提：在一个长度为 n 的数组 nums 里的所有数字都在 0～n-1 的范围内 遍历数组中每一个数字num，并且将这个数字作为索引，将取到的数字nums[num]进行标记进行标记，当一个待标记的nums[num]已经是一个负数，说明这个num就是重复的数字 public int
刷题学习—数据结构(字符串、栈和队列、链表、数组与矩阵、哈希表、位运算)
np131131的博客
11-01 1195
参考了很多大佬的题解，仅作为自己学习笔记用。 数据结构相关 第一章 字符串 1.字符串循环移位包含 给定两个字符串 s1 和 s2，要求判定 s2 是否能够被 s1 做循环移位得到的字符串包含 解决办法：s1 进行循环移位的结果是 s1s1 的子字符串，因此只要判断 s2 是否是 s1s1 的子字符串即可。 public class Exer { public boolean strStr(String s1, String s2) { StringBuilder sb = new
【专题】—【数学】—【矩阵快速幂】
ACer-TomatoFish
08-16 4722
参考文章： http://blog.csdn.net/shiwei408/article/details/8818386
利用矩阵快速幂转换的题目
weixin_34326179的博客
09-01 159
T-shirt 题目地址：https://nanti.jisuanke.com/t/28873 JSZKC is going to spend his vacation! His vacation hasNNdays. Each day, he can choose a T-shirt to wear. Obviously, he doesn't want to wear a si...
矩阵快速幂 题目整理 HDU5015、CF392 C、HDU4686、CF385 E
qcwlmqy的博客
07-30 207
矩阵快速幂 矩阵 在数学中，矩阵（Matrix）是一个按照长方阵列排列的复数或实数集合。
十个利用矩阵乘法解决的经典题目
weixin_30243533的博客
03-24 247
好像目前还没有这方面题目的总结。这几天连续看到四个问这类题目的人，今天在这里简单写一下。这里我们不介绍其它有关矩阵的知识，只介绍矩阵乘法和相关性质。不要以为数学中的矩阵也是黑色屏幕上不断变化的绿色字符。在数学中，一个矩阵说穿了就是一个二维数组。一个n行m列的矩阵可以乘以一个m行p列的矩阵，得到的结果是一个n行p列的矩阵，其中的第i行第j列位置上的数等于前一个矩阵第i行上的m个数与后一个矩阵...
矩阵快速幂模板+例题
yangzijiangac的博客
02-06 365
矩阵快速幂这个算法，理解起来很容易，但是我之前自己写的代码有bug，也是因为上课不听课，对形参和实参没理解，平常用的都是全局变量，是不是实参影响不大，这次定义一个结构体的矩阵，矩阵需要初始化为0，然后，因为形参和实参没怎么理解，导致输出的答案差异很大，前提是矩阵快速幂，矩阵需要初始化数组，并不会默认为0 重要的事情讲三遍： 矩阵需要初始化，并不会默认为0！ 矩阵需要初始化，并不会默认为0！ ...
矩阵快速幂（快速幂）模板题目Decoding Genome
名字不能相同
05-22 1020
这类题目的难点在于，数据量大，需要用到高精度，也就是快速幂 https://codeforces.com/group/vFwRVj9WjO/contest/328371/problem/H #include <iostream> #include <cstdio> #include <cstdlib> #include <cstring> #include <cmath> #include <algorithm> #define M
矩阵快速幂模板+例题（快速入门）
serendipity的博客
03-23 1084
矩阵快速幂模板+例题 本来早在一个月前就应该把这给掌握的，硬是得等到卡题了，才想到来补知识 文章目录矩阵快速幂模板+例题模板例题一、HDU 1575：Tr A二、POJ 3070：Fibonacci三、HNUCM 1620: Fy's dota2 模板 传送门：矩阵快速幂基础讲解 单纯的求某一矩阵的快速幂是没有什么意义的，主要是由于所求的递推式的数比较大时，如果使用暴力递推，必定会超时，这时矩阵快...
HDU 5950(Recursive sequence)
大宝哥spring
10-30 2947
Recursive sequence Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 29    Accepted Submission(s): 19 Problem Description Farmer John
【矩阵快速幂】经典题 hdu2157 how many ways、woj642 Lost In WHU
Hi_jiaxinwei的博客
04-24 367
How many ways?? Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 3151    Accepted Submission(s): 1202 Problem Description 春天到了, HDU校
矩阵论题型解析与求解策略
"该资源是关于矩阵论的题型总结，主要针对北邮程云鹏版教材，涵盖了矩阵论中的重要概念和解题方法。" 矩阵论是线性代数的一个分支，它深入研究矩阵的性质及其在数学、物理学、工程学等领域的应用。这个题型总结主要...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值