8.9 暑假集训——矩阵快速幂

最新推荐文章于 2023-01-09 14:43:26 发布
最新推荐文章于 2023-01-09 14:43:26 发布
阅读量237 收藏
点赞数
分类专栏： 17年暑假集训 矩阵快速幂
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_37280023/article/details/77075338
版权

暂时还有题目没做出来，把已经写出来的题目先整理出来吧，顺便理一下矩阵快速幂的常见用法和思路，为做下面的题目做一个缓冲

由等式求递推矩阵，加速

F-M斐波那契数列

题意： 类似斐波那契数列给出一种数列的新定义：
F[0] = a
F[1] = b
F[n] = F[n-1] * F[n-2] ( n > 1 )
给出a, b, n，（ 0 <= a, b, n <= 10^9），要求输出F[n]的值%(1e9+7)

思路： 写出几个F[n]，a,b, ab, ab2, a2b3, a3b5, a5b8……
可以 观察出 这个序列的次数是斐波那契数列，由于n比较大，用矩阵快速幂加速一下~ 然后因为斐波那契数做次数可能比较大，所以用到费马小定理，取余（1e9+6）
这里写图片描述

之后问题的 关键只要写出斐波那契数列的递推矩阵 就可以啦~ 所以关键如图
这里写图片描述

话说这个名字前的M是代表matrix嘛？ 正式比赛的时候会给这种提示咩？？？ 

#include <iostream>
#include <cstdio>
#include <cstring>
#include <string>
#include <algorithm>
#include <cmath>
using namespace std;
typedef long long LL;
struct Squ{
    LL s[2][2];  //一开始WA了还找不出是哪里，后来看到这里设成了int觉得都不想对自己说什么了……
};
Squ e,t;
const LL mod=1e9+7;
LL a,b,n;

LL quickpow(LL a,LL k){
    a%=mod;
    LL ans=1;
    while(k){
        if(k&1)
            ans=(ans*a)%mod;
        a=(a*a)%mod;
        k/=2;
    }
    return ans;
}

Squ mul(Squ a,Squ b){
    Squ ans;
    for(int i=0;i<2;i++)
    for(int j=0;j<2;j++){
        ans.s[i][j]=0;
        for(int k=0;k<2;k++){
            ans.s[i][j]+=a.s[i][k]*b.s[k][j];
            ans.s[i][j]%=(mod-1);
        }
    }
    return ans;
}

Squ squickpow(Squ a,LL k){
    Squ ans=e;
    while(k){
        if(k&1) ans=mul(ans,a);
        a=mul(a,a);
        k/=2;
    }
    return ans;
}

int main(){
    for(int i=0;i<2;i++)    //初始化单位矩阵
        for(int j=0;j<2;j++)
            if(i==j) e.s[i][j]=1;
            else e.s[i][j]=0;
    t.s[0][0]=t.s[0][1]=t.s[1][0]=1;
    t.s[1][1]=0;
    while(~scanf("%lld%lld%lld",&a,&b,&n)){
        a%=mod,b%=mod;
        if(n==0){  //两种情况单独处理
            printf("%lld\n",a%mod); continue;
        }
        else if(n==1){
            printf("%lld\n",b%mod); continue;
        }
        Squ ans=squickpow(t,n-2);  //对应图中的n-2次
        LL fn=(ans.s[0][0]+ans.s[0][1])%(mod-1);
        LL fn1=(ans.s[1][0]+ans.s[1][1])%(mod-1);
        long long res=(quickpow(a,fn1)*quickpow(b,fn))%mod;
        printf("%lld\n",res);
    }
    return 0;
}

Recursive sequence

题意：　这里写图片描述
告诉F[1],F[2]和n,（F[1],F[2],n<2^31）, 求F[n]

思路： 知等式求递推矩阵，看到递推关系式中的i4，第一要想办法把它化成（i-1）4，所以可得
这里写图片描述

#include <iostream>
#include <cstdio>
#include <cstring>
#include <string>
#include <cstdio>
#include <algorithm>
#include <cmath>
using namespace std;
typedef long long LL;
const LL mod=2147493647;
struct Squ{
    LL s[7][7];
};
Squ squ={{1,2,1,4,6,4,1,1,0,0,0,0,0,0,0,0,1,4,6,4,1,0,0,0,1,3,3,1,0,0,0,0,1,2,1,0,0,0,0,0,1,1,0,0,0,0,0,0,1}};  //注意像这样的结构体初始化，只能在它声明的时候
Squ e;

Squ mul(Squ a,Squ b){  //模板
    Squ ans;
    for(int i=0;i<7;i++)
    for(int j=0;j<7;j++){
        ans.s[i][j]=0;
        for(int k=0;k<7;k++){
            ans.s[i][j]+=a.s[i][k]*b.s[k][j];
            ans.s[i][j]%=mod;
        }
    }
    return ans;
}
Squ squickpow(Squ a,LL k){ //模板
    Squ ans=e;
    while(k){
        if(k&1) ans=mul(ans,a);
        a=mul(a,a);
        k/=2;
    }
    return ans;
}

int main(){
    for(int i=0;i<7;i++)
        for(int j=0;j<7;j++)
            if(i==j) e.s[i][j]=1;
            else e.s[i][j]=0;
    LL a,b,n;
    int t;
    scanf("%d",&t);
    while(t--){
        scanf("%lld%lld%lld",&n,&a,&b);
        a%=mod;
        b%=mod;
        if(n==1) printf("%lld\n",a);
        else if(n==2) printf("%lld\n",b);
        else{
            Squ aa=squickpow(squ,n-2); //想一想，为什么是n-2
            LL f0[7]={b,a,16,8,4,2,1};
            LL fn[7];
            for(int i=0;i<7;i++){
                fn[i]=0;
                for(int j=0;j<7;j++){
                    fn[i]+=aa.s[i][j]*f0[j];
                    fn[i]%=mod;
                }
            }
            printf("%lld\n",fn[0]%mod);
        }
    }
    return 0;
}

找思路，套模板

B-How many ways

题意： 知道n个点，输入m个有向边，给定s,e,k（0 < n <= 20, m <= 100，k < 20），求从s到e且经过k条边的路线 有多少种方案
思路： 就是离散数学中学的邻接矩阵， 设给的初始有向边的矩阵为 C，则C的n次方中，s,e对应矩阵的点c’[s][e]表示从s到e长度为k的路线个数。注意一条边可以走多次。
所以整理一下，只要知道初始矩阵的k次幂后，对应c’[s][e]的大小就是所求的解了~

#include <iostream>
#include <cstdio>
#include <cstring>
#include <string>
#include <cmath>
#include <algorithm>
using namespace std;
struct Squ{
    int s[25][25];
};
Squ e,squ;
int n,m;
const int mod=1000;

Squ mul(Squ a,Squ b){  //模板，矩阵相乘
    Squ squ;
    for(int i=0;i<n;i++)
    for(int j=0;j<n;j++){
        squ.s[i][j]=0;
        for(int k=0;k<n;k++)
            squ.s[i][j]+=a.s[i][k]*b.s[k][j];
            squ.s[i][j]%=mod;
    }
    return squ;
}

Squ squickpow(Squ squ,int k){ //模板，矩阵快速幂 精髓
    Squ ans=e;
    while(k){
        if(k&1)
            ans=mul(ans,squ);
        k/=2; squ=mul(squ,squ);
    }
    return ans;
}

int main(){
    for(int i=0;i<25;i++)
        for(int j=0;j<25;j++)
            if(i==j) e.s[i][j]=1;
            else e.s[i][j]=0;   //初始化单位矩阵

    while(~scanf("%d%d",&n,&m)&&n){
        memset(squ.s,0,sizeof(squ.s));
        for(int i=0;i<m;i++){
            int a,b;
            scanf("%d%d",&a,&b);
            squ.s[a][b]=1;
        }
        int t;
        scanf("%d",&t);
        while(t--){
            int s,e,k;
            scanf("%d%d%d",&s,&e,&k);  //所以正式解题就这三步啊orz
            Squ a=squickpow(squ,k); //
            int res=a.s[s][e]%mod; //
            printf("%d\n",res);
        }
    }
    return 0;
}

Matrix Power Series

题意： 给一个n*n，(n ≤ 30) 的矩阵A，和次数 k (k ≤ 109) 模数m (m < 104)，最后求如图的矩阵S，矩阵中每个元素取余m
这里写图片描述
思路：

思路一：二分

这里写图片描述

思路二：矩阵套矩阵

这里写图片描述

#include <iostream>
#include <cstdio>
#include <cstring>
#include <string>
#include <algorithm>
#include <cmath>
using namespace std;
typedef long long LL;
struct Squ{
    LL s[35][35];
};
struct mSqu{
    Squ s[2][2];
}msqu;

Squ squ,e,t;
LL n,mod,k;

Squ mul(Squ a,Squ b){  //矩阵乘法
    Squ ans;
    for(int i=0;i<n;i++)
    for(int j=0;j<n;j++){
        ans.s[i][j]=0;
        for(int k=0;k<n;k++){
            ans.s[i][j]+=a.s[i][k]*b.s[k][j];
            ans.s[i][j]%=mod;
        }
    }
    return ans;
}

mSqu mmul(mSqu a,mSqu b){  //元素为矩阵的矩阵乘法，类似矩阵
    mSqu ans;
    for(int i=0;i<2;i++)
    for(int j=0;j<2;j++){
        ans.s[i][j]=t;
        for(int k=0;k<2;k++){
            Squ temp=mul(a.s[i][k], b.s[k][j]);
            for(int ii=0;ii<n;ii++){
                for(int jj=0;jj<n;jj++){
                    ans.s[i][j].s[ii][jj]+=temp.s[ii][jj];
                    ans.s[i][j].s[ii][jj]%=mod;
                }
            }
        }
    }
    return ans;
}

mSqu msquick(mSqu a,int k){  //元素为矩阵的矩阵快速乘法
    mSqu ans;

    for(int i=0;i<2;i++)
        for(int j=0;j<2;j++)
            if(i==j) ans.s[i][j]=e;
            else ans.s[i][j]=t;
    while(k){
        if(k&1)
            ans=mmul(ans,a);  
            a=mmul(a,a);
            k/=2;
        }
    return ans;
}

int main(){  //只有一组数据
    scanf("%lld%lld%lld",&n,&k,&mod);
    for(int i=0;i<n;i++)
    for(int j=0;j<n;j++){
        scanf("%lld",&squ.s[i][j]);
        t.s[i][j]=0;
        if(i==j) e.s[i][j]=1;
        else e.s[i][j]=0;
    }
    msqu.s[0][0]=squ;
    msqu.s[0][1]=msqu.s[1][1]=e;
    msqu.s[1][0]=t;
    mSqu ans=msquick(msqu,k);
    Squ res=t;

    for(int i=0;i<n;i++){
        for(int j=0;j<n;j++){
        res.s[i][j]=ans.s[0][0].s[i][j]+ans.s[0][1].s[i][j];
        if(i==j) res.s[i][j]--;  //注意这里-1，因为这种算法多加了一个单位矩阵，需要减去
        res.s[i][j]%=mod;
        if(j) printf(" ");
        printf("%lld",res.s[i][j]);
        }
        printf("\n");
    }
    return 0;
}

用矩阵快速幂打辅助的数学分析题

So Easy!

题意： 给出a,m,b,n, （0< a, m < 2^15, (a-1)^2< b < a^2, 0 < b, n < 2^31)，求这里写图片描述
思路：
这里写图片描述

#include <iostream>
#include <cstring>
#include <cstdio>
#include <string>
#include <algorithm>
#include <cmath>
using namespace std;
typedef long long LL;
struct Squ{
    LL s[2][2];
};
Squ e={{1,0,0,1}};
Squ squ;
LL a,b,n,mod;

Squ mul(Squ a,Squ b){
    Squ ans;
    for(int i=0;i<2;i++)
    for(int j=0;j<2;j++){
        ans.s[i][j]=0;
        for(int k=0;k<2;k++){
            ans.s[i][j]+=a.s[i][k]*b.s[k][j];
            ans.s[i][j]%=mod;
        }
    }
    return ans;
}

Squ squickpow(Squ squ,int k){
    Squ ans=e;
    while(k){
        if(k&1) ans=mul(ans,squ);
        k/=2;
        squ=mul(squ,squ);
    }
    return ans;
}

int main(){
    while(~scanf("%lld%lld%lld%lld",&a,&b,&n,&mod)){
        a%=mod,b%=mod;
        LL ans=0;
        if(n==1)    ans=2*a%mod;
        else if(n==2) ans=(2*a*a+2*b)%mod;
        else{
            squ.s[0][0]=(2*a)%mod;
            squ.s[0][1]=(b-a*a+mod*mod)%mod; //注意这里的取模  第一次交WA是因为只加了一个mod，只有加mod*mod才可以
            squ.s[1][0]=1;
            squ.s[1][1]=0;
            Squ res=squickpow(squ,n-2);
            LL f1=2*a%mod;
            LL f2=(2*a*a+2*b)%mod;
            ans=res.s[0][0]%mod*f2%mod+res.s[0][1]*f1%mod;
            ans%=mod;
        }
        printf("%lld\n",ans%mod);

    }
    return 0;
}
「已注销」
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
【从零开始学习 SystemVerilog】8.9、SystemVerilog 约束—— solve before
努力，加油每一天鸭~
09-12 529
【从零开始学习 SystemVerilog】8.9、SystemVerilog 约束—— solve before
openssh8.9版本
03-08
《OpenSSH 8.9版本在CentOS 7上的应用与配置详解》 OpenSSH是Secure Shell的开源实现，它提供了安全的网络服务，包括加密的远程登录、文件传输等功能。OpenSSH 8.9p1是其最新的稳定版本，带来了多项改进和安全修复...
快速幂与矩阵快速幂
RDJ_Widow的博客
01-21 167
快速幂 ll qpow(ll x,ll k) { ll res=1; x=x%mod; while (k)//把次数k看成二进制数 { if (k&1) res=res*x%mod;//遇到二进制位为1，就乘进答案里 x=x*x%mod;//加大权值 k>>=1; } retur...
矩阵快速幂（个人心得+白书板子）
DDDDDDDDDD4444444444的博客
05-31 117
矩阵快速幂的大体思想与快速幂想法一致，只是普通的数乘，变成了矩阵相乘，所以问题的关键在于矩阵乘法上，即（看代码） ll pow(ll x,ll n,ll mod){ ll res=1;//此处的1变为单位阵 while(n>0){ if(n&1)res=res*x%mod;//res*x变成矩阵乘法 x=x*x%mod;//同上 n&g...
8.9zju集训日记
ysk850874665的博客
08-10 154
       和新队员的第一次比赛，前期开题方向基本正确，签到的速度比较快，中期读了旋转卡壳，矩阵和km的三道题目，都有一定的想法，但三个人意见不同没有往一个方向想，但其实旋转卡壳和km的题目思路几乎都对了，但是旋转卡壳少一点点思路，km又缺少比较优秀的模板，一个题目给了n=1000，O（n^3）的复杂度，虽然时间给了七秒，但还是以为过不了，事实上赛后补题，发现自己的做法只跑了七百多毫秒。   ...
0019——8.9——面向对象概述
weixin_43332626的博客
08-09 70
1.对象本身就是一级指针，类是产生对象的模板
【矩阵论】8. 常用矩阵总结——单阵，正规阵，幂0阵，幂等阵，循环阵
努力，奋斗！
01-09 2838
单阵： 单阵谱公式：单阵逆公式，单阵函数公式 单阵条件： 充分条件：正规阵是单阵，方阵有n个不同的特向，A有n个无关特向，k重根必须有k个无关特向 充要条件：0化式判定 正规阵定义 正规阵特点 常见正规阵：对角阵，H阵，U阵 正规阵构造方法：倍数法则，平移法则，U相似，多项式正规 幂0阵条件，幂0阵特根特向，幂0阵不是单阵 幂等阵一定是单阵A幂等I-A幂等r(A)+r(I-A)=n 循环阵：基本循环阵，循环阵可由基本循环阵表示，傅里叶阵
暑假集训-8.9总结
weixin_30469895的博客
08-09 69
学习内容：数学知识 今日完成题数：4 今日看书情况：20页 学习算法的总结： 看书笔记：http://note.youdao.com/noteshare?id=ab09364903107dea47bc0a60fbf9c9d6 今日做题总结： https://www.acwing.com/activity/content/code/content/8...
暑假集训 - 8.9 总结
weixin_30782331的博客
08-09 77
学习内容：set 完成题数：3题 看书情况：6页 做题总结： AcWing 146：https://www.cnblogs.com/buhuiflydepig/p/11326141.html codeforces 854c：https://www.cnblogs.com/buhuiflydepig/p/11329858.html 心得： 不会用一些STL容器，T...
C#入门8.9——Main方法
乔培宸的博客
08-16 1262
Main方法是程序的入口点。程序从这里开始，也是从这里结束。 执行过程：程序在执行编写的源代码时，是先找Main方法，然后开始执行Main方法中“{”开始后的第一句代码，依次执行，如果遇到Main方法中有调用其他的方法时，便会根据方法名称找到定义方法的代码，然后执行这个方法内的代码，执行完这个方法后，再返回到Main方法继续执行，直到遇到Main方法的结束符“}”，执行程序结束。 C#的Mai
Proteus 8.9SP2仿真软件
perseverance51博客
04-22 1924
Proteus 8.9 SP2仿真软件
CentOS7升级OPENSSH至8.9P1
04-25
本文将详细讲解如何在`CentOS7`系统中将`OpenSSH`升级到8.9P1版本，这是一个重要的安全更新，无需进行复杂的编译过程。 首先，我们理解`OpenSSH`是开放源代码的实现，它提供了一个安全的网络通信协议，用于替代不...
顶尖电子称条码秤LS系列称固件，最新固件，固件
最新发布
08-28
顶尖电子称条码秤LS系列称固件，最新固件，固件
爬取淘宝热销(热门)咖啡商品信息数据集
08-28
本资源旨在收集淘宝平台上热销咖啡商品的相关信息，包括但不限于：店铺所在省份、城市位置、商品名称、销售价格、销量数据、单价（以人民币计）、付款人数、是否提供包邮服务、是否来自天猫平台以及满减优惠的详细信息。这些信息均来源于淘宝平台的公开透明数据，可供学习参考之用。请注意，本资源仅用于学术研究和个人学习目的，不得用于商业用途或其他非法活动。
爬取的淘宝新品推荐男装商品信息公开透明数据集
08-28
本资源专注于收集淘宝平台最新推荐的男装商品信息，内容涵盖商品的店铺所在省份、城市位置、商品的名称、销售价格、累积销量、单价（以人民币计价）、付款的顾客人数、是否提供包邮服务、是否为天猫平台的商品，以及相关的满减优惠情况。这些详细的数据点均来源于淘宝平台的公开透明信息，经过精确抓取和整理，旨在为分析电商平台上的新品推荐策略和消费者购买行为提供实用数据。 这些数据严格遵循淘宝平台的公开政策和隐私保护原则获取，确保了信息的合法性与合规性。然而，本资源仅作为学习参考之用，意在帮助研究人员、市场分析师或学生等理解电商领域的商品推荐机制、销售动态及市场趋势。 任何将此数据用于商业目的或其他未授权的活动都是不恰当的，甚至可能触犯相关法律条款。 在使用这些数据进行学术研究或个人学习时，用户应自觉遵守相关法律法规，尊重数据来源和版权，正确引用数据源，并不得用于任何形式的商业盈利。
STM32F4的环境配置，可直接运行STM32F4系列的芯片
08-28
STM32F4的环境配置，可直接运行STM32F4系列的芯片
技术资料PCB资料模拟电子技术资料03j.zip
08-28
技术资料PCB资料模拟电子技术资料03j.zip
RoseHA 8.9 for Windows快速安装及配置指南
"RoseHA 8.9 for Windows快速安装说明v1.8" 本文档是RoseHA 8.9 for Windows的快速安装指南，旨在帮助用户顺利部署和管理高可用性集群。RoseHA是一款旨在确保业务连续性的软件，通过在Windows环境下提供故障转移和...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值