kuangbin专题十九 UVA10689 矩阵快速幂

最新推荐文章于 2020-11-24 12:54:48 发布

Start_to_crazy

最新推荐文章于 2020-11-24 12:54:48 发布

阅读量304

点赞数

CC 4.0 BY-SA版权

分类专栏：矩阵

本文链接：https://blog.csdn.net/Start_to_crazy/article/details/79082664

矩阵专栏收录该内容

9 篇文章

订阅专栏

本文介绍了一种解决类似斐波那契数列问题的方法，通过矩阵快速幂算法求解特定位置的数值，并给出完整的C++代码实现。

题意:
给定一个类似斐波那契数列的数列的头两项，求的n项取模10^m。
题解：
斐波那契数列变形，直接堆上去模板就可以过了，这道题我WA了两次都是因为装逼，认为一看到模板题就不用动手写一下，结果头两项乘错了，晕。

#include<stdio.h>
#include<string.h>
#include<math.h>
#include<algorithm>
using namespace std;
#define LL long long int
struct node
{
    int m[2][2];
    node()
    {
        memset(m,0,sizeof(m));
    }   
}; 
int MOD,A,B,n,m;
node cla(node a,node b)
{
    node c;
    for(int i=0;i<2;i++)
        for(int j=0;j<2;j++)
            for(int k=0;k<2;k++)
            if(a.m[i][k]&&b.m[k][j])
            {
                c.m[i][j]+=a.m[i][k]*b.m[k][j];
                c.m[i][j]%=MOD;
            }
    return c;
}
void POW(int k)
{
    node a,c;
    for(int i=0;i<2;i++) c.m[i][i]=1;
    a.m[0][0]=1,a.m[0][1]=1;
    a.m[1][0]=1,a.m[1][1]=0;
    while(k)
    {
        if(k&1) c=cla(c,a);
        a=cla(a,a);
        k>>=1;
    }
    printf("%d\n",(B*c.m[0][0]%MOD+A*c.m[0][1]%MOD)%MOD);
}
int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d%d%d%d",&A,&B,&n,&m);
        MOD=pow(10,m);
        POW(n-1);
    }
}