TZOJ——4954：矩阵游戏（快速幂+费马小定理+指数循环节+逆元）

描述：

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

输入：

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

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

输出：

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

样例输入：

3 4 1 3 2 6

样例输出

85

提示

样例中的矩阵为：

1 4 7 10
26 29 32 35
76 79 82 85

---

题目分析：题目中给出3个式子：式1表示矩阵左上角为1；式2表示矩阵的第2列至第m列都由前1列推出；式3表示矩阵的第2行至第n行的第1个数都由上1行的最后1个数推出。

刚开始我以为这是一道与矩阵相关的题，直到我看见n和m的数据范围。。。

数据范围这么大，肯定得优化，由费马小定理+指数循环节可以进一步的优化指数的大小，详细内容可以点击这里。指数问题优化掉了，接下来就是公式推导了：

最终可得：

令x=a^(m-1)c，y=a^(m-1)d，p=(m-1)b（a=1）或b[a^(m-1)-1]/(a-1)（a!=1）,q=(n-1)x（x=1）或[x^(n-1)-1]/(x-1)（x!=1），得F[n,m]=x^(n-1)[a^(m-1)+p]+(y+p)q。 

最后还剩一个问题：除法怎么办。一般情况下遇到除法，可以用快速幂将除以一个数转化为乘以这个数的逆元，这样就不用担心精度问题了。

代码示例：

#include<iostream>
#include<cstring>
using namespace std;
#define ll long long
const int N=1e6+6,M=1e9+7;
char sn[N],sm[N];
ll quick(ll a,ll b) //快速幂
{
    ll ans=1;
    a%=M;
    while(b)
    {
        if(b&1)
            ans=ans*a%M;
        a=a*a%M,b>>=1;
    }
    return ans;
}
void func1(ll &a1,ll &a2,char s[],int l)
{
    a1=a2=0;
    for(int i=0;i<l;i++)
    {
        a1=(a1*10+s[i]-'0')%M;
        a2=(a2*10+s[i]-'0')%(M-1);  //费马小定理
    }
}
int main()
{
    ll a,b,c,d,m1,m2,n1,n2,x,y,p,q,aa;
    int ln,lm;
    scanf("%s%s%lld%lld%lld%lld",sn,sm,&a,&b,&c,&d);
    ln=strlen(sn);
    lm=strlen(sm);
    func1(n1,n2,sn,ln);
    func1(m1,m2,sm,lm);
    aa=quick(a,m2-1);
    x=aa*c%M,y=aa*d%M;
    if(a==1)
        p=b*((m1-1+M)%M)%M;
    else
        p=b*(aa-1)%M*quick(a-1,M-2)%M;
    if(x==1)
        q=x*((n1-1+M)%M)%M;
    else
        q=(quick(x,(n2-1+M)%M)-1+M)%M*quick(x-1,M-2)%M;
    printf("%lld",(quick(x,(n2-1+M)%M)*(aa+p)%M+(y+p)%M*q%M)%M);
    return 0;
}