luogu1397 [NOI2013]矩阵游戏 (等比数列求和)

一个比较显然的等比数列求和，但有一点问题就是n和m巨大..

考虑到他们是在幂次上出现，所以可以模上P-1（费马小定理）

但是a或c等于1的时候，不能用等比数列求和公式，这时候就要乘n和m，又要变成模P

所以我们一开始就模P*(P-1)好了...

很大，要用龟速乘

#include<bits/stdc++.h>
#define CLR(a,x) memset(a,x,sizeof(a))
#define MP make_pair
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> pa;
const int maxn=1;
const ll P=1e9+7,PP=P*(P-1);

inline ll fmul(ll x,ll y){
    ll re=0;
    while(y){
        if(y&1) re=(re+x)%PP;
        x=(x+x)%PP,y>>=1;
    }return re;
}

inline ll rd(){
    ll x=0;char c=getchar();int neg=1;
    while(c<'0'||c>'9'){if(c=='-') neg=-1;c=getchar();}
    while(c>='0'&&c<='9') x=(fmul(x,10)+c-'0')%PP,c=getchar();
    return x*neg;
}

inline ll fpow(ll x,ll y){
    ll re=1;
    while(y){
        if(y&1) re=re*x%P;
        x=x*x%P,y>>=1;
    }
    return re;
}

int main(){
    // freopen("testdata.in","r",stdin);
    ll n=(rd()-1+PP)%PP,m=(rd()-1+PP)%PP,a=rd(),b=rd(),c=rd(),d=rd();
    ll e=fpow(a,m),f;
    if(a!=1) f=(e-1)*fpow(a-1,P-2)%P*b%P;
    else f=m%P*b%P;
    ll g=e*c%P,h=(f*c%P+d)%P;
    ll i=fpow(g,n),j;
    if(g!=1) j=(i-1)*fpow(g-1,P-2)%P*h%P;
    else j=n%P*h%P;
    printf("%lld\n",((e*(i+j)+f)%P+P)%P);
    return 0;
}