通过这个题真的学到了不少东西,最起码矩阵快速幂算是入门了,普通快速幂也彻底明白了(以前都是打模板),了解了费马小定理 关键点 求(a^fib[b])%p 的值其中p是素数,0p-1 那么上式 (a^fib[b])%p = (a^(p-1)*a^(p-1)*....*a^(p-1)*a^m)%...
分类:
其他好文 时间:
2015-04-06 06:16:37
阅读次数:
132
裸题,最简单fib的应用模板,算是新技能get 吧。 其实和快速幂差不多了,只是矩阵代替的递推式。 1 #include 2 #include 3 #include 4 using namespace std; 5 const int maxn = 1005; 6 struct node...
分类:
其他好文 时间:
2015-04-06 00:53:32
阅读次数:
201
传送门:Gift题意:由n(n#include #include #include #include #include #define LL long long#define N 25#define mod 1000000007using namespace std;/***************...
分类:
其他好文 时间:
2015-04-06 00:51:39
阅读次数:
243
题目1081:递推数列时间限制:1 秒内存限制:32 兆特殊判题:否提交:5885解决:800题目描述:给定a0,a1,以及an=p*a(n-1) + q*a(n-2)中的p,q。这里n >= 2。 求第k个数对10000的模。输入:输入包括5个整数:a0、a1、p、q、k。输出:第k个数a(k)对...
分类:
其他好文 时间:
2015-04-05 18:41:37
阅读次数:
128
大致题意:简单题意就不解释了。
不过可以建议大家可以先做POJ 3070,先学会快速幂的基本思想。
没有做过的可以查看我的博客:点击打开链接
然后我们已经会使用矩阵快速幂求解A^k,则如何求解 A+A^2+A^3……
这里还是运用了二分的思想 ,例如:
令S(N)=A+A^2+……+A^N;
则S(6)=A+A^2+A^3+A^4+A^5+A^6=(1+A^3)...
分类:
其他好文 时间:
2015-04-04 12:17:44
阅读次数:
168
C - Matrix Power Series Time Limit:3000MS Memory Limit:131072KB 64bit IO Format:%I64d & %I64uSubmit Status Practice POJ 3233DescriptionGiven a n × n m...
分类:
其他好文 时间:
2015-04-04 12:01:26
阅读次数:
128
题意:给一个n*n的矩阵A,求k次幂之和S=A+A2+A3+ … +Ak思路:矩阵快速幂。#include#include#includeusing namespace std;typedef struct node{ int matrix[55][55];}Matrix;Matrix a,s...
分类:
其他好文 时间:
2015-04-04 11:59:15
阅读次数:
134
其实矩阵构造就是对公式的化简,最后运用矩阵快速幂求值下面来看一题Everybody knows Fibonacci numbers, now we are talking about the Tribonacci numbers: T[0] = T[1] = T[2] = 1; T[n] = T[n...
分类:
其他好文 时间:
2015-04-04 10:36:09
阅读次数:
178
题目链接:hdu 4965,题目大意:给你一个 n*k 的矩阵 A 和一个 k*n 的矩阵 B,定义矩阵 C= A*B,然后矩阵 M= C^(n*n),矩阵中一切元素皆 mod 6,最后求出 M 中所有元素的和。题意很明确了,便赶紧敲了个矩阵快速幂的模板(因为编程的基本功不够还是调试了很久),然.....
分类:
其他好文 时间:
2015-04-03 21:02:14
阅读次数:
144
矩阵快速幂,1001. Fibonacci 2求斐波那契第n项!毕竟数据量太大!http://soj.sysu.edu.cn/show_problem.php?pid=1001&cid=1740 1 #include 2 #include 3 4 using namespace std; 5 ...
分类:
其他好文 时间:
2015-04-02 15:01:25
阅读次数:
148