矩阵乘法-洛谷P2233 [HNOI2002] 公交车路线

最新推荐文章于 2024-04-09 09:22:59 发布

原创最新推荐文章于 2024-04-09 09:22:59 发布 · 589 阅读

0 ·

CC 4.0 BY-SA版权

矩阵乘法专栏收录该内容

13 篇文章

订阅专栏

本文介绍了一种利用矩阵快速幂解决特定问题的方法，并通过一个示例详细展示了如何实现矩阵的定义、乘法运算及快速幂的过程。适用于理解矩阵运算基本原理及提高算法效率。

https://daniu.luogu.org/problem/show?pid=2233
矩阵这个很显然啊；
然后直接快速幂就好了；
至于为什么，这个就是矩阵的基本性质；
可以看我相关的博客；
然后到了E点就不懂，直接在乘法的时候处理一下就好了

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cstring>
#include<map>
#define Ll long long
using namespace std;
int a[9][9]
{//A B C D E F G H
{0,0,0,0,0,0,0,0,0},
{0,0,1,0,0,0,0,0,1},//A
{0,1,0,1,0,0,0,0,0},//B
{0,0,1,0,1,0,0,0,0},//C
{0,0,0,1,0,1,0,0,0},//D
{0,0,0,0,1,0,1,0,0},//E
{0,0,0,0,0,1,0,1,0},//F
{0,0,0,0,0,0,1,0,1},//G
{0,1,0,0,0,0,0,1,0},//H
};
struct jv{
    int m[9][9];
    jv(){memset(m,0,sizeof m);}
}ans,c;
int n,mo=1e3;
jv cheng(jv x,jv y){
    jv z;
    for(int i=1;i<=8;i++)
    for(int j=1;j<=8;j++)
    for(int k=1;k<=8;k++)if(k!=5)
    z.m[i][j]=(z.m[i][j]+x.m[i][k]*y.m[k][j])%mo;
    return z;
}
int main()
{
    scanf("%d",&n);
    n--;
    for(int i=1;i<=8;i++)
    for(int j=1;j<=8;j++)
        c.m[i][j]=a[i][j];
    ans=c;
    while(n){
        if(n&1)ans=cheng(ans,c);
        n>>=1;
        c=cheng(c,c);
    }
    printf("%d",ans.m[1][5]);
}