【POJ 3233】矩阵快速幂+二分求等比数列前N项和
来源:互联网 发布:java compareto方法 编辑:程序博客网 时间:2024/05/16 19:33
http://acm.hust.edu.cn/vjudge/contest/view.action?cid=20851#overview
开了一套矩阵和高斯消元方面的题。
POJ3233这题给一个矩阵A,然后求(A+A^2+A^3+……+A^k)%m
一看到这个首先想到的是等比数列求和,但是矩阵不能除,所以要考虑另一种方案。
把这个矩阵拆成两个部分,就是(A+A^2+A^3+……+A^(k/2))+A^(k/2)*(A+A^2+A^3+……+A^(k/2)),这样就可以对(A+A^2+A^3+……+A^(k/2))进行递归计算。如果k是奇数就再+A^k。
这个方法也可以用于等比数列前N项和取模问题,如果要用等比数列前N项和公式的话要涉及求逆元或者卢卡斯定理等比较麻烦的问题,不如采用二分法来得简单。
matrix expmod(matrix a,int n)//快速幂{e.init();while (n){if(n & 1) e=e*a;a=a*a;n>>=1;}return e;}matrix solve(int k)//二分{if(k==1) return a;c=solve(k>>1);matrix ans=c+c*expmod(a,(k>>1));if(k & 1) ans=ans+expmod(a,k);return ans;}
- 【POJ 3233】矩阵快速幂+二分求等比数列前N项和
- poj 1845 Sumdiv(数论:欧拉函数+二分求等比数列前n项和+快速幂取模)
- POJ 3233 二分求等比数列 矩阵快速幂
- poj 1845(等比数列前n项和及快速幂)
- POJ 3233 Matrix Power Series --二分求矩阵等比数列和
- 求等比数列的前n项和
- POJ 3233 矩阵快速幂 等比数列
- HDU 1588 Gauss Fibonacci(矩阵快速幂+二分求等比数列和)
- HDU - 1588 Gauss Fibonacci (矩阵快速幂+二分求等比数列和)
- POJ 1845 Sumdiv (唯一分解定理+求等比数列前n项和)
- POJ 3233 Matrix Power Series (矩阵快速幂+等比数列二分求和)
- 【POJ 1845】 Sumdiv (整数唯分+约数和公式+二分等比数列前n项和+同余)
- 快速矩阵幂+二分等比数列求和
- 矩阵快速二分求n次幂
- POJ 3233-Matrix Power Series(矩阵快速幂+二分求矩阵和)
- POJ 3233 Matrix Power (矩阵快速幂+等比数列求和)
- 等比数列的前n项和
- POJ 3233 矩阵运算,等比数列二分求和,矩阵
- 移植的一种思路
- Flex 4里的fx、mx以及s命名空间
- os Crash Log 分析汇总
- Solairs分区
- Eclipse导入源码工程
- 【POJ 3233】矩阵快速幂+二分求等比数列前N项和
- GNU内嵌汇编 ARM 版
- 第四周上机程序报告一
- java 加载类和反射机制的总结
- 两个JavaBean之间的复制,但是bean的属性名不一致。
- 在 ASP.NET 中虚拟目录指向远程共享时发生"System.Security.SecurityException: 安全性错误"的解决办法
- 2013-3-22.2
- 第四周任务——三角形类
- linux 驱动开发