矩阵游戏

【题目描述】

婷婷是个喜欢矩阵的小朋友,有一天她想用电脑生成一个巨大的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;
}