[bzoj3240][矩阵乘法]矩阵游戏

Description

婷婷是个喜欢矩阵的小朋友,有一天她想用电脑生成一个巨大的n行m列的矩阵(你不用担心她如何存储)。她生成的这个矩阵满足一个神奇的性质：若用F[i][j]来表示矩阵中第i行第j列的元素，则F[i][j]满足下面的递推式:

F[1][1]=1 F[i,j]=a*F[i][j-1]+b (j!=1) F[i,1]=c*F[i-1][m]+d (i!=1)
递推式中a,b,c,d都是给定的常数。

现在婷婷想知道F[n][m]的值是多少,请你帮助她。由于最终结果可能很大，你只需要输出F[n][m]除以1,000,000,007的余数。

Input

一行有六个整数n,m,a,b,c,d。意义如题所述

Output

包含一个整数，表示F[n][m]除以1,000,000,007的余数

Sample Input

3 4 1 3 2 6

Sample Output

85

HINT

样例中的矩阵为：

1 4 7 10

26 29 32 35

76 79 82 85

1<=N,M<=10^1000 000,a<=a,b,c,d<=10^9

题解

由于不想写高精度取模不想写费马小定理导致卡了一晚上常数超级难受
很容易看出，二进制快速幂是不可实现的
那么就写一个十进制快速幂
矩阵状态优化一下，矩阵乘法优化一下
优化了两小时终于卡过去了。。

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long LL;
const LL mod=1000000007;
int sum[2][1100000],length[2];
inline int read()
{
    int f=1,x=0;char ch=getchar();
    while(ch<'0' || ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0' && ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
inline void del(int op,int b)
{
    sum[op][1]-=b;
    for(int i=1;i<=length[op];i++)while(sum[op][i]<0)sum[op][i]+=10,sum[op][i+1]--;
}
struct matrix
{
    LL m[2][2];
    matrix(){memset(m,0,sizeof(m));}
}stx,sty,tmp;
inline matrix mul(matrix u,matrix v,int n,int m,int p)
{
    matrix ret;
    ret.m[0][0]=(u.m[0][0]*v.m[0][0]+u.m[0][1]*v.m[1][0])%mod;
    ret.m[1][0]=(u.m[1][0]*v.m[0][0]+u.m[1][1]*v.m[1][0])%mod;
    ret.m[1][1]=(u.m[1][0]*v.m[0][1]+u.m[1][1]*v.m[1][1])%mod;
    return ret;
}
inline matrix pow_mod_two(matrix u,int b)
{
    register matrix ret;
    ret.m[0][0]=ret.m[1][1]=1;
    while(b)
    {
        if(b&1)ret=mul(ret,u,2,2,2);
        u=mul(u,u,2,2,2);b/=2;
    }
    return ret;
}
char ch[2][1100000];
inline matrix pow_mod_ten(matrix u,int k,int op,int len)
{
    register matrix ret;
    ret.m[0][0]=ret.m[1][1]=1;
    while(k<=len)
    {   
        ret=mul(ret,pow_mod_two(u,sum[op][k]),2,2,2);
        u=pow_mod_two(u,10);
        k++;
    }
    return ret;
}
int a,b,c,d;
int main()
{
    scanf("%s%s",ch[0]+1,ch[1]+1);
    int lenx=strlen(ch[0]+1),leny=strlen(ch[1]+1);
    a=read();b=read();c=read();d=read();
    stx.m[0][0]=a;stx.m[1][0]=b;stx.m[1][1]=1;
    sty.m[0][0]=c;sty.m[1][0]=d;sty.m[1][1]=1;
    tmp.m[0][0]=1;tmp.m[0][1]=1;
    length[0]=lenx;
    for(int i=1;i<=lenx;i++)sum[0][lenx-i+1]=ch[0][i]-'0';
    length[1]=leny;
    for(int i=1;i<=leny;i++)sum[1][leny-i+1]=ch[1][i]-'0';
    del(1,1);
    register matrix P=pow_mod_ten(stx,1,1,length[1]);
    register matrix Q=mul(P,sty,2,2,2);
    del(0,1);
    Q=pow_mod_ten(Q,1,0,length[0]);
    Q=mul(Q,P,2,2,2);
    tmp=mul(tmp,Q,1,2,2);
    printf("%lld\n",tmp.m[0][0]);
    return 0;
}