HDU 2276 Kiki & Little Kiki 2

最新推荐文章于 2021-02-01 12:30:32 发布

原创最新推荐文章于 2021-02-01 12:30:32 发布 · 765 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#ini

ACM-矩阵专栏收录该内容

6 篇文章

订阅专栏

本文介绍了一种利用矩阵快速幂解决特定类型递推问题的方法。通过C++代码示例展示了如何进行矩阵乘法运算及矩阵快速幂计算，适用于求解大规模数值问题。

http://acm.hdu.edu.cn/showproblem.php?pid=2276

#include<iostream> using namespace std; #define N 101 #define FF(i,n) for(int i=0;i<n;i++) int ans[N]; int init[N][N]; void matrixMul(int a[][N],int b[][N],int n,int mod){ int i,j,k; int buf[N][N]={0}; FF(i,n)FF(k,n) if(a[i][k])FF(j,n)if(b[k][j]) buf[i][j]=(buf[i][j]+a[i][k]*b[k][j])%mod; FF(i,n)FF(j,n) a[i][j]=buf[i][j]; } void matrixMul(int a[],int b[][N],int n,int mod){ int i,j; int buf[N]={0}; FF(i,n)FF(j,n) buf[i]=(buf[i]+a[j]*b[i][j])%mod; FF(i,n) a[i]=buf[i]; } void matrixMul(int ans[],int init[][N],int n,int m,int mod){ for(;m;m>>=1){ if(m&1) matrixMul(ans,init,n,mod); matrixMul(init,init,n,mod); } } int main(void){ int m,n,i,j,mod; char ch[N]; while(~scanf("%d",&m)){ scanf("%s",ch); n=strlen(ch); mod=2; FF(i,n) ans[i]=ch[i]-'0'; memset(init,0,sizeof(init)); init[0][0]=init[0][n-1]=1; for(i=1;i<n;i++) init[i][i-1]=init[i][i]=1; matrixMul(ans,init,n,m,mod); FF(i,n) printf("%d",ans[i]); printf("/n"); } }