HDOJ 233 Matrix 5015【矩阵快速幂】
2015-08-15 19:19
330 查看
233 MatrixTime Limit: 10000/5000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others)Total Submission(s): 1355 Accepted Submission(s): 806 Problem Description In our daily life we often use 233 to express our feelings. Actually, we may say 2333, 23333, or 233333 ... in the same meaning. And here is the question: Suppose we have a matrix called 233 matrix. In the first line, it would be 233, 2333, 23333... (it means a0,1 = 233,a0,2 = 2333,a0,3 = 23333...) Besides, in 233 matrix, we got ai,j = ai-1,j +ai,j-1( i,j ≠ 0). Now you have known a1,0,a2,0,...,an,0, could you tell me an,m in the 233 matrix? Input There are multiple test cases. Please process till EOF. For each case, the first line contains two postive integers n,m(n ≤ 10,m ≤ 109). The second line contains n integers, a1,0,a2,0,...,an,0(0 ≤ ai,0 < 231). Output For each case, output an,m mod 10000007. Sample Input 1 1 1 2 2 0 0 3 7 23 47 16 Sample Output 234 2799 72937 Hint Source 2014 ACM/ICPC Asia Regional Xi'an Online Recommend hujie | We have carefully selected several similar problems for you: 5390 5389 5388 5387 5386 |
= ai-1,j+ai,j-1 ....就根据这个公式算出所有元素。然后输出anm
。
解题分析:
此题一看就是裸快速幂~因为m比较大。开始的时候是想直接推算出整个矩阵,后来发现推不出来。然后就得想到推下一列。就是有第一列推第二列,第二列推第三列。。。。等等。。
我们先写几个元素
a[1]
a[2]
a[3]
想到有233~得出来这个数啊。
23
a[1]
a[2]
a[3]
3
推下一列有:
23*10+3
a[1]+23*10+3
a[2]+a[1]+23*10+3
a[3]+a[2]+a[1]+23*10+3
3
然后转移矩阵A就可以推出来了:
10 0 0 0 1
10 1 0 0 1
10 1 1 0 1
10 1 1 1 1
0 0 0 0 1
然后对A求快速幂,再乘以a数列就可以了
#include <stdio.h> #include <math.h> #include <vector> #include <queue> #include <string> #include <string.h> #include <stdlib.h> #include <iostream> #include <algorithm> using namespace std; typedef long long LL; const int Modn = 10000007; const int MAXN = 15; LL res[MAXN][MAXN]; int N; void Matmul(LL X[MAXN][MAXN],LL Y[MAXN][MAXN]) { LL t[MAXN][MAXN]={0}; for(int i=0;i<N;i++) for(int k=0;k<N;k++) if(X[i][k]) for(int j=0;j<N;j++) t[i][j]=(t[i][j]+X[i][k]*Y[k][j]%Modn)%Modn; for(int i=0;i<N;i++) for(int j=0;j<N;j++) X[i][j]=t[i][j]; } void Matrix(LL X[MAXN][MAXN],LL n) { for(int i=0;i<N;i++) for(int j=0;j<N;j++) res[i][j]=(i==j); while(n){ if(n&1) Matmul(res,X); Matmul(X,X); n>>=1; } } int main() { LL n; LL m; while(scanf("%lld%lld",&n,&m)!=EOF){ LL a[MAXN]; LL A[MAXN][MAXN]; memset(A,0,sizeof(A)); N=n+2; for(int i=1;i<=n;i++) scanf("%lld",&a[i]); a[0]=23;a[n+1]=3; for(int i=0;i<N;i++){ if(i!=N-1)A[i][0]=10; for(int j=1;j<=i;j++){ if(i!=N-1)A[i][j]=1; } A[i][N-1]=1; } Matrix(A,m); LL ans=0; //for(int i=0;i<N;i++) for(int j=0;j<N;j++){ ans=(ans+a[j]*res[N-2][j]%Modn)%Modn; } printf("%lld\n",ans); } return 0; }
相关文章推荐
- final, finally, finalize
- VS2013常用快捷键
- [mahout in action]mahout中用于推荐引擎的数据结构(二)
- wamp配置
- 使用easyBCD 安装 Win8 + Ubuntu 13.10 双系统
- VS2012 快捷键
- 十道海量数据处理面试题与十个方法大总结
- zoj3344 第一类斯特林数+java大数
- 文件工具类
- Linux C语言程序设计(十三)——Linux上的系统与文件命令
- poj 3417 Network 【LCA】【树中增新边后 求每条树边被环所覆盖的次数】
- LeetCode(215) Kth Largest Element in an Array
- Android之Intent
- 【C#】TreeView应用
- 不是技术牛人,如何拿到国内IT巨头的Offer
- 【J2SE】这才是java虚引用(PhantomReference)
- 算法学习之排序算法(四)(希尔排序)
- 萧墙HTML5手机发展之路(51)——jquerymobile在提高页面访问速度
- C++之友元
- 二叉树俩个子节点的公共父节点