记录编号 |
577949 |
评测结果 |
AAAAAAAAAAAAAAAAAAAA |
题目名称 |
[NOI 2013]矩阵游戏 |
最终得分 |
100 |
用户昵称 |
op_组撒头屯 |
是否通过 |
通过 |
代码语言 |
C++ |
运行时间 |
1.674 s |
提交时间 |
2022-12-13 15:14:35 |
内存使用 |
1.83 MiB |
显示代码纯文本
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll mod=1000000007;
struct sdf{
ll p[2][2];
sdf(){
for (int i=0;i<=1;i++){
for (int j=0;j<=1;j++)p[i][j]=0;
}
}
sdf operator*(const sdf&x){
sdf tmp;
for (int i=0;i<=1;i++){
for (int j=0;j<=1;j++){
for (int k=0;k<=1;k++){
tmp.p[i][j]=(tmp.p[i][j]+p[i][k]*x.p[k][j]%mod)%mod;
}
}
}
return tmp;
}
}A[12],B,C,D[12],E,bs,tmp;
ll a,b,c,d;
string n,m;
int nl,ml;
sdf pw(sdf x,int y){
if (!y)return bs;
if (y&1)return pw(x,y-1)*x;
sdf P=pw(x,y/2);return P*P;
}
void fst1(int pt){
if (pt>=nl)return ;
C=pw(C,10);
C=C*A[n[pt]-'0'];
fst1(pt+1);
}
void fst2(int pt){
if (pt>=ml)return ;
E=pw(E,10);
E=E*D[m[pt]-'0'];
fst2(pt+1);
}
int main(){
freopen ("matrixb.in","r",stdin);
freopen ("matrixb.out","w",stdout);
cin>>m>>n;
scanf("%lld%lld%lld%lld",&a,&b,&c,&d);
ml=m.length(),nl=n.length();
for (int i=ml-1;i>=0;i--){
if (m[i]=='0')m[i]='9';
else{
m[i]--;break;
}
}
for (int i=nl-1;i>=0;i--){
if (n[i]=='0')n[i]='9';
else{
n[i]--;break;
}
}
A[1].p[0][0]=a;A[1].p[0][1]=b;A[1].p[1][0]=0;A[1].p[1][1]=1;
B.p[0][0]=c;B.p[0][1]=d;B.p[1][0]=0;B.p[1][1]=1;
bs.p[0][0]=bs.p[1][1]=1;bs.p[0][1]=bs.p[1][0]=0;
A[0]=D[0]=C=E=bs;
for (int i=2;i<=9;i++)A[i]=A[i-1]*A[1];
if (n[0]!='0')fst1(0);
else fst1(1);
D[1]=C*B;
for (int i=2;i<=9;i++)D[i]=D[i-1]*D[1];
if (m[0]!='0')fst2(0);
else fst2(1);
E=E*C;
ll ans=(E.p[0][0]+E.p[0][1])%mod;
printf("%lld\n",ans);
return 0;
}