快速幂模板
来源:互联网 发布:ug软件的教学大纲 编辑:程序博客网 时间:2024/05/22 03:05
下面是 m^n % k 的快速幂:
// m^n % k
下面是矩阵快速幂:
//HOJ 3493
/*===================================*/
|| 快速幂(quickpow)模板
|| P 为等比,I 为单位矩阵
|| MAX 要初始化!!!!
||
/*===================================*/
/*****************************************************/
// m^n % k
int quickpow(int m,int n,int k){ int b = 1; while (n > 0) { if (n & 1) b = (b*m)%k; n = n >> 1 ; m = (m*m)%k; } return b;}
下面是矩阵快速幂:
//HOJ 3493
/*===================================*/
|| 快速幂(quickpow)模板
|| P 为等比,I 为单位矩阵
|| MAX 要初始化!!!!
||
/*===================================*/
/*****************************************************/
#include <cstdio>const int MAX = 3;typedef struct{ int m[MAX][MAX];} Matrix;Matrix P = {5,-7,4, 1,0,0, 0,1,0, };Matrix I = {1,0,0, 0,1,0, 0,0,1, }; Matrix matrixmul(Matrix a,Matrix b) //矩阵乘法{ int i,j,k; Matrix c; for (i = 0 ; i < MAX; i++) for (j = 0; j < MAX;j++) { c.m[i][j] = 0; for (k = 0; k < MAX; k++) c.m[i][j] += (a.m[i][k] * b.m[k][j])%9997; c.m[i][j] %= 9997; } return c;} Matrix quickpow(long long n){ Matrix m = P, b = I; while (n >= 1) { if (n & 1) b = matrixmul(b,m); n = n >> 1; m = matrixmul(m,m); } return b;} /*************************************/int main(){ Matrix re; int f[3] = {2,6,19}; long long n; while (scanf("%I64d",&n) && n != 0) { if (n == 1) printf("1\n"); else if (n <= 4) printf("%d\n",f[n-2]); else { re = quickpow(n - 4); printf("%d\n",(((re.m[0][0]*f[2]) + (re.m[0][1]*f[1]) + (re.m[0][2]*f[0])) %9997 + 9997) % 9997); } } return 0;}
0 0
- 快速幂模板
- 矩阵快速幂模板
- 快速幂 ---模板
- 快速幂模板
- 快速幂模板
- 矩阵快速幂--模板
- 矩阵快速幂 模板
- 矩阵快速幂模板
- 矩阵快速幂模板
- 快速幂模板
- 快速幂小模板
- 【矩阵快速幂模板】
- 矩阵快速幂模板
- 快速幂模板
- 快速幂算法模板
- 矩阵快速幂模板
- 矩阵快速幂模板
- 整数快速幂【模板】
- 墨水屏
- 11 装饰
- (转)Linux中gcc g++常用编译选项以及makefile的编写
- 取消Eclipse等号、分号、空格代码自动补全
- 异质链表保存学生信息
- 快速幂模板
- P91
- java基础篇——抽象类的特点简述
- SQLiteDatabase的CRUD操作实现
- iOS 的枚举
- 黑马程序员—C重点--总结一些琐碎的知识点
- 曲面屏幕手机
- Mac10.10重新安装PHP问题
- ICE(一)——安装