codevs 1281 Xn数列 (矩阵乘法)

Posted 一入OI深似海

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了codevs 1281 Xn数列 (矩阵乘法)相关的知识,希望对你有一定的参考价值。

/*
再来个题练练手 
scanf longlong 有bug.......
 
*/
#include<cstdio>
#include<iostream>
#include<cstring>
#define ll long long 
using namespace std;
ll n,m,x,y,x0,g;
ll f[3][3],a[3][3];
ll slow_mul(ll a,ll b,ll c)
{
    ll ans=0;
    a=a%c;b=b%c;
    while(b)
      {
          if(b&1)
            {
                b--;
                ans+=a;
                ans%=c;
          }
        a<<=1;a%=c;b>>=1;
      }
    return ans;
}
void mul(ll a[3][3],ll b[3][3])
{
    ll c[3][3];
    memset(c,0,sizeof(c));
    for(int i=1;i<=2;i++)
      for(int j=1;j<=2;j++)
        for(int k=1;k<=2;k++)
          c[i][j]=(c[i][j]+slow_mul(a[i][k],b[k][j],m))%m;
    for(int i=1;i<=2;i++)
      for(int j=1;j<=2;j++)
        a[i][j]=c[i][j];
}
int main()
{
    cin>>m>>x>>y>>x0>>n>>g;
    //scanf("%ld%ld%ld%ld%ld%ld",&m,&x,&y,&x0,&n,&g);
    f[1][1]=x0;f[1][2]=1;a[1][1]=x;a[1][2]=0;a[2][1]=y;a[2][2]=1;
    while(n)
      {
          if(n&1)mul(f,a);
          mul(a,a);
          n>>=1;
      }
    printf("%ld\n",f[1][1]%g);
    return 0;
}

 

以上是关于codevs 1281 Xn数列 (矩阵乘法)的主要内容,如果未能解决你的问题,请参考以下文章

CODEVS1281Xn数列

codevs 1281 Xn数列

Xn数列(codevs 1281)

C++之路进阶——矩阵乘法(Xn数列)

codevs1281 Xn数列

codevs1281 矩阵乘法 快速幂 !!!手写乘法取模!!! 练习struct的构造函数和成员函数