BZOJ3240 NOI2013矩阵游戏（数论）

　　必修五题。

// luogu-judger-enable-o2
#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;
int read()
{
    int x=,f=;char c=getchar();
    while (c<''||c>'') {if (c=='-') f=-;c=getchar();}
    while (c>=''&&c<='') x=(x<<)+(x<<)+(c^),c=getchar();
    return x*f;
}
#define N 1000010
#define P 1000000007
int n,m,a[N],b[N],c,d,u,v,p,q,ans,x,y,k,tmp1,tmp2;
int ksm(int a,int k)
{
    if (k==) return ;
    int tmp=ksm(a,k>>);
    if (k&) return 1ll*tmp*tmp%P*a%P;
    else return 1ll*tmp*tmp%P;
}
int inv(int x){return ksm((x+P)%P,P-);}
int main()
{
    char ch=getchar();
    while (ch>=''&&ch<='') k=(10ll*k+(ch^))%(P-),tmp1=(10ll*tmp1+(ch^))%P,ch=getchar();
    while (ch<''||ch>'') ch=getchar();
    while (ch>=''&&ch<='') p=(10ll*p+(ch^))%(P-),tmp2=(10ll*tmp2+(ch^))%P,ch=getchar();
    u=read(),v=read(),c=read(),d=read();
    if (u>) p=(p+P-)%(P-),p=ksm(u,p),q=1ll*v*inv(u-)%P*(p+P-)%P;
    else tmp2=(tmp2+P-)%P,q=1ll*tmp2*v%P,p=;
    x=p,y=q;
    p=1ll*p*c%P;q=(d+1ll*c*q%P)%P;
    if (p>) k=(k+P-)%(P-),q=1ll*q*inv(p-)%P,p=ksm(p,k),q=1ll*q*(p+P-)%P,ans=(p+q)%P;
    else tmp1=(tmp1+P-)%P,ans=(+1ll*tmp1*q)%P;
    ans=(1ll*ans*x%P+y+P)%P;
    cout<<ans;
    return ;
}