HDU 1757 A Simple Math Problem

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

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

0 ·

CC 4.0 BY-SA版权

文章标签：

#math #ini

ACM-矩阵专栏收录该内容

6 篇文章

订阅专栏

本文探讨了在ACM HDU平台上的矩阵乘法问题，通过使用高效的矩阵乘法算法，结合模运算，实现了对大规模矩阵的快速处理。重点介绍了三种不同情况下的矩阵乘法优化方法，并通过实例展示了算法的应用。

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

#include<iostream> using namespace std; #define N 10 #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 n=10; int m,mod; int i,j; while(~scanf("%d%d",&m,&mod)){ memset(init,0,sizeof(init)); FF(i,n) ans[i]=9-i; FF(j,n) scanf("%d",&init[0][j]); for(int i=1;i<n;i++) init[i][i-1]=1; if(m<10){ printf("%d/n",ans[9-m]); continue; } m-=9; matrixMul(ans,init,n,m,mod); printf("%d/n",ans[0]); } }