正题
题目链接:https://ac.nowcoder.com/acm/contest/1111/A
题目大意
2∗22*22∗2的矩阵AAA,求AnA^nAn
解题思路
我们将其十进制的转化过来,如
A2016=(A103)2∗(A10)∗(A)6A^{2016}=(A^{10^3})^2*(A^{10})*(A)^6A2016=(A103)2∗(A10)∗(A)6
这样时间复杂度为O(log10(n)∗10)O(log_{10}(n)*10)O(log10(n)∗10)
codecodecode
#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
const int Size=2,XJQ=7;
struct matrix{int a[Size][Size];
}f,ans;
int m,l;
char n[100010];
matrix operator *(matrix a,matrix b)
{matrix c;memset(c.a,0,sizeof(c.a));for(int i=0;i<2;i++)for(int j=0;j<2;j++)for(int k=0;k<2;k++)(c.a[i][j]+=a.a[i][k]*b.a[k][j]%XJQ)%=XJQ;return c;
}
int main()
{scanf("%s",n);l=strlen(n)-1;for(int i=0;i<=l/2;i++)swap(n[i],n[l-i]);for(int i=0;i<2;i++)for(int j=0;j<2;j++)scanf("%d",&f.a[i][j]);for(int i=0;i<2;i++)ans.a[i][i]=1;for(int i=0;i<=l;i++){for(int j=1;j<=n[i]-'0';j++)ans=ans*f;f=f*f*f*f*f*f*f*f*f*f;}for(int i=0;i<2;i++){for(int j=0;j<2;j++)printf("%d ",ans.a[i][j]);putchar('\n');}
}