2019 牛客暑期多校 B generator 1 （矩阵快速幂+倍增）

最新推荐文章于 2019-08-17 20:58:00 发布

weixin_30649641

最新推荐文章于 2019-08-17 20:58:00 发布

阅读量55

点赞数

原文链接：http://www.cnblogs.com/Lis-/p/11285289.html

版权

题目：https://ac.nowcoder.com/acm/contest/885/B

题意：给你x0,x1，让你求出xn，递推式时xn=a*xn-1+b*xn-2

思路：这个n特别大，我自己没有摸清欧拉降幂的性质，瞎套了，然后其实因为底数是一个矩阵，并不能运用这一定理，但是这个n又这么大，我们就可以使用倍增

这里用2倍增有点麻烦，我们就直接用10倍增，然后这个递推式很明显就能看出是一个2*2的矩阵快速幂，然后求解即可

#include<bits/stdc++.h>
#define maxn 1000005
using namespace std;
typedef long long ll;
ll x0,x1,a,b,mod;
ll mod1;
char str[maxn];
char s[maxn];
struct jz//结构体写法的矩阵快速幂
{
    long long num[2][2];
    jz() { memset(num,0,sizeof(num)); }
    jz(ll a,ll b,ll c,ll d){
        num[0][0]=a;
        num[0][1]=b;
        num[1][0]=c;
        num[1][1]=d;
    };
    jz operator*(const jz &P)const {
        jz ans;
        for(int k=0;k<2;k++)
            for(int i=0;i<2;i++)
                for(int j=0;j<2;j++)
                    ans.num[i][j]=(ans.num[i][j]+num[i][k]*P.num[k][j]%mod)%mod;
        return ans;
    }
}COE,ans,unit;
int T_T;
jz F, A;
jz B, T;
jz pOw(jz X,long long m)//矩阵快速幂
{
    jz ans;
    ans.num[0][0]=1;
    ans.num[1][1]=1; 
    for(;m;m>>=1,X=X*X)
        if(m&1)
            ans=ans*X;
    return ans;
}
int main(){
    scanf("%lld%lld%lld%lld",&x0,&x1,&a,&b);
    scanf("%s%lld",str,&mod);
    unit=jz(x1,x0,0,0);
    COE=jz(a,1,b,0);
    for(int i=strlen(str)-1;i>=0;i--){//倍增
        unit=unit*pOw(COE,str[i]-'0');
        COE=pOw(COE,10); 
    }
    printf("%lld",unit.num[0][1]);
}