hdu1757---A Simple Math Problem(矩阵)
来源:互联网 发布:独角仙 数码兽数据库 编辑:程序博客网 时间:2024/05/16 10:07
problem Description
Lele now is thinking about a simple function f(x).
If x < 10 f(x) = x.
If x >= 10 f(x) = a0 * f(x-1) + a1 * f(x-2) + a2 * f(x-3) + …… + a9 * f(x-10);
And ai(0<=i<=9) can only be 0 or 1 .
Now, I will give a0 ~ a9 and two positive integers k and m ,and could you help Lele to caculate f(k)%m.
Input
The problem contains mutiple test cases.Please process to the end of file.
In each case, there will be two lines.
In the first line , there are two positive integers k and m. ( k<2*10^9 , m < 10^5 )
In the second line , there are ten integers represent a0 ~ a9.
Output
For each case, output f(k) % m in one line.
Sample Input
10 9999 1 1 1 1 1 1 1 1 1 1 20 500 1 0 1 0 1 0 1 0 1 0
Sample Output
45 104
Author
linle
Source
2007省赛集训队练习赛(6)_linle专场
Recommend
lcy | We have carefully selected several similar problems for you: 1588 3117 2276 2256 2254
很明显的矩阵递推题
/************************************************************************* > File Name: hdu1757.cpp > Author: ALex > Mail: zchao1995@gmail.com > Created Time: 2015年03月11日 星期三 19时57分20秒 ************************************************************************/#include <map>#include <set>#include <queue>#include <stack>#include <vector>#include <cmath>#include <cstdio>#include <cstdlib>#include <cstring>#include <iostream>#include <algorithm>using namespace std;const double pi = acos(-1);const int inf = 0x3f3f3f3f;const double eps = 1e-15;typedef long long LL;typedef pair <int, int> PLL;int mod;class MARTIX{ public: int mat[105][105]; MARTIX(); MARTIX operator + (const MARTIX &b)const; MARTIX operator * (const MARTIX &b)const; MARTIX & operator = (const MARTIX &b);};MARTIX :: MARTIX(){ memset (mat, 0, sizeof(mat));}MARTIX MARTIX :: operator + (const MARTIX &b)const{ MARTIX ret; for (int i = 0; i < 10; ++i) { for (int j = 0; j < 10; ++j) { ret.mat[i][j] = mat[i][j] + b.mat[i][j]; ret.mat[i][j] %= mod; } } return ret;}MARTIX MARTIX :: operator * (const MARTIX &b)const{ MARTIX ret; for (int i = 0; i < 10; ++i) { for (int j = 0; j < 10; ++j) { ret.mat[i][j] = 0; for (int k = 0; k < 10; ++k) { ret.mat[i][j] += mat[i][k] * b.mat[k][j]; ret.mat[i][j] %= mod; } } } return ret;}MARTIX & MARTIX :: operator = (const MARTIX &b){ for (int i = 0; i < 10; ++i) { for (int j = 0; j < 10; ++j) { this -> mat[i][j] = b.mat[i][j]; } } return *this;}MARTIX fastpow(MARTIX A, LL n){ MARTIX ans; for (int i = 0; i < 10; ++i) { for (int j = 0; j < 10; ++j) { ans.mat[i][j] = (i == j); } } while (n) { if (n & 1) { ans = ans * A; } n >>= 1; A = A * A; } return ans;}void Debug(MARTIX A){ for (int i = 0; i < 10; ++i) { for (int j = 0; j < 10; ++j) { printf("%d ", A.mat[i][j]); } printf("\n"); }}int main (){ LL k; while (~scanf("%lld%d", &k, &mod)) { if (k < 10) { printf("%lld\n", k); continue; } MARTIX A; for (int i = 0; i < 10; ++i) { scanf("%d", &A.mat[i][0]); } for (int i = 0; i < 9; ++i) { A.mat[i][i + 1] = 1; } MARTIX ans = fastpow(A, k - 9); MARTIX F; for (int i = 0; i < 10; ++i) { F.mat[0][i] = 10 - i - 1; }// Debug(F); ans = F * ans; printf("%d\n", ans.mat[0][0]); } return 0;}
- HDU1757 A Simple Math Problem 矩阵应用
- hdu1757 A Simple Math Problem(矩阵乘法)
- hdu1757 A Simple Math Problem 矩阵乘法
- hdu1757---A Simple Math Problem(矩阵)
- A Simple Math Problem hdu1757
- hdu1757 A Simple Math Problem
- HDU1757 A Simple Math Problem
- HDU1757 A Simple Math Problem
- hdu1757 - A Simple Math Problem
- A Simple Math Problem HDU1757
- Hdu1757 - A Simple Math Problem - 矩阵快速幂
- HDU1757-A Simple Math Problem(矩阵快速幂)
- hdu1757-- A Simple Math Problem(矩阵快速幂优化)
- hdu1757 - A Simple Math Problem 矩阵快速幂
- hdu1757 A Simple Math Problem(矩阵快速幂)
- hdu1757 A Simple Math Problem (矩阵快速幂)
- hdu1757 A Simple Math Problem(矩阵快速幂)
- [HDU1757]A Simple Math Problem(矩阵快速幂)
- 初尝2D骨骼动画编辑工具SPINE
- 全排列生成算法之字典序
- 嵌入式C语言中printf函数的使用(二)——SWO引脚的巧妙应用
- poj 1287||zoj 1372 Networking 最小生成树 kruskal 克鲁斯卡尔算法
- JAVA2核心技术中的GBC类
- hdu1757---A Simple Math Problem(矩阵)
- IOS 调整内存中的图片大小
- HDU 5011-Game(尼姆博弈)
- CSDN-markdown编辑器简洁的语法与呈现效果对照2
- 编程好的规则
- Spinner入门1
- Configure,Makefile.am, Makefile.in, Makefile文件之...
- 彩色图像--色彩空间 综述
- PAT 1044. Shopping in Mars (25)