【题目描述】
婷婷是个喜欢矩阵的小朋友,有一天她想用电脑生成一个巨大的n行m列的矩阵(你不用担心她如何存储)。她生成的这个矩阵满足一个神奇的性质:若用F[i][j]来表示矩阵中第i行第j列的元素,则F[i][j]满足下面的递推式:
递推式中a,b,c,d都是给定的常数。
现在婷婷想知道F[n][m]的值是多少,请你帮助她。由于最终结果可能很大,你只需要输出F[n][m]除以1,000,000,007的余数。
【输入格式】
包含一行有六个整数n,m,a,b,c,d。意义如题所述。
【输出格式】
包含一个整数,表示F[n][m]除以1,000,000,007的余数。
【样例输入】
3 4 1 3 2 6
【样例输出】
85
【提示】
样例中的矩阵为:
【数据范围】
矩阵快速幂
#include<fstream>
#include<string>
using namespace std;
const long long p=1000000000+7;
ifstream cin("matrixb.in");
ofstream cout("matrixb.out");
string s1,s2;
int a,b,c,d;
long long C[2][2],A[2][2],x[2][2],y[2][2],z[2][2],t[2][2];
long long D[20][2][2];
struct bigint
{
int d[1000010];
int len;
}
n,m;
void st(long long a[2][2])
{
a[0][0]=1;
a[0][1]=0;
a[1][0]=0;
a[1][1]=1;
}
void multiple(long long a[2][2],long long b[2][2],long long c[2][2])
{
long long x1=(b[0][0]*c[0][0]+b[0][1]*c[1][0])%p;
long long x2=(b[0][0]*c[0][1]+b[0][1]*c[1][1])%p;
long long x3=(b[1][0]*c[0][0]+b[1][1]*c[1][0])%p;
long long x4=(b[1][0]*c[0][1]+b[1][1]*c[1][1])%p;
a[0][0]=x1;
a[0][1]=x2;
a[1][0]=x3;
a[1][1]=x4;
}
void power(long long x[2][2],long long y[2][2],bigint&n)
{
st(x);
st(D[0]);
long long t[2][2];
for(int i=1;i<10;i++)
multiple(D[i],D[i-1],y);
for(int i=0;i<n.len;i++)
{
multiple(x,x,D[n.d[i]]);
multiple(t,D[1],D[9]);
for(int i=1;i<10;i++)
multiple(D[i],D[i-1],t);
}
}
void minus1(bigint&x)
{
int temp=0;
x.d[0]--;
for(int i=0;i<x.len;i++)
{
x.d[i]-=temp;
if(x.d[i]<0)
{
temp=1;
x.d[i]+=10;
}
else return;
}
if(!x.d[x.len-1])
x.len--;
}
int main()
{
// scanf("%s %s %d %d %d %d", s1, s2, &a, &b, &c, &d);
//freopen()
ios::sync_with_stdio(false);
cin>>s1>>s2>>a>>b>>c>>d;
n.len=s1.length();
for(int i=n.len-1;i>=0;i--)
n.d[n.len-i-1]=s1[i]-'0';
m.len=s2.length();
for(int i=m.len-1;i>=0;i--)
m.d[m.len-i-1]=s2[i]-'0';
A[0][0]=a;
A[0][1]=b;
A[1][0]=0;
A[1][1]=1;
C[0][0]=c;
C[0][1]=d;
C[1][0]=0;
C[1][1]=1;
minus1(m);
power(x,A,m);
multiple(y,C,x);
minus1(n);
power(z,y,n);
multiple(y,x,z);
cout<<(y[0][0]+y[0][1])%p;
return 0;
}