Another kind of Fibonacce(矩阵快速幂,HDU3306)
Another kind of Fibonacci
Time Limit: 3000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others)
Total Submission(s): 2470 Accepted Submission(s): 987
Each test case will contain three integers , N, X , Y .
N : 2<= N <= 231 – 1
X : 2<= X <= 231– 1
Y : 2<= Y <= 231 – 1
3 2 3
196
#include <iostream>
#include <cstdio>
using namespace std;
typedef long long ll;
#define mod 10007
ll n;
ll x;
ll y;
struct matrix
{
ll v[4][4];
};
matrix mul(matrix m1, matrix m2)
{
matrix c;
for (int i = 0; i<4; i++)
for (int j = 0; j<4; j++)
{
c.v[i][j] = 0;
for (int k = 0; k<4; k++)
c.v[i][j] += (m1.v[i][k] * m2.v[k][j]) % mod;
c.v[i][j] %= mod;
}
return c;
}
matrix qupow(matrix A)
{
n--;
matrix res = A, m = A;
while (n)
{
if (n & 1)
res = mul(res, m);
n >>= 1;
m = mul(m, m);
}
return res;
}
void solve(matrix A)
{
matrix temp = qupow(A);
ll ans = 0;
for (int i = 0; i<4; i++)
ans += temp.v[i][0]%mod;
ans%=mod;
printf("%lld\n",ans);
//cout << ans << endl;
}
int main()
{
matrix A;
ll a, b, c;
while (~scanf("%lld%lld%lld",&n,&x,&y))
{
a = x*x%mod;
b = y*y%mod;
c = x*y * 2%mod;
//matrix A = { 1,0,0,0 ,1,a,1,x,0,b,0,0 0,c,0,y };
A.v[0][0]=1;
for(int i=1;i<4;i++)
A.v[0][i]=0;
A.v[1][0]=1;
A.v[1][1]=a%mod;
A.v[1][2]=1;
A.v[1][3]=x%mod;
A.v[2][0]=0;
A.v[2][1]=b%mod;
A.v[2][2]=0;
A.v[2][3]=0;
A.v[3][0]=0;
A.v[3][1]=c%mod;
A.v[3][2]=0;
A.v[3][3]=y%mod;
/*for(int i=0;i<4;i++)
for(int j=0;j<4;j++)
cin>>A.v[i][j];*/
solve(A);
}
return 0;
}
Another kind of Fibonacce(矩阵快速幂,HDU3306)的更多相关文章
- hdu3306 Another kind of Fibonacci【矩阵快速幂】
转载请注明出处:http://www.cnblogs.com/KirisameMarisa/p/4187670.html 题目链接:http://acm.hdu.edu.cn/showproblem. ...
- 矩阵快速幂 HDU 4565 So Easy!(简单?才怪!)
题目链接 题意: 思路: 直接拿别人的图,自己写太麻烦了~ 然后就可以用矩阵快速幂套模板求递推式啦~ 另外: 这题想不到或者不会矩阵快速幂,根本没法做,还是2013年长沙邀请赛水题,也是2008年Go ...
- 51nod 算法马拉松18 B 非010串 矩阵快速幂
非010串 基准时间限制:1 秒 空间限制:131072 KB 分值: 80 如果一个01字符串满足不存在010这样的子串,那么称它为非010串. 求长度为n的非010串的个数.(对1e9+7取模) ...
- 51nod 1113 矩阵快速幂
题目链接:51nod 1113 矩阵快速幂 模板题,学习下. #include<cstdio> #include<cmath> #include<cstring> ...
- 【66测试20161115】【树】【DP_LIS】【SPFA】【同余最短路】【递推】【矩阵快速幂】
还有3天,今天考试又崩了.状态还没有调整过来... 第一题:小L的二叉树 勤奋又善于思考的小L接触了信息学竞赛,开始的学习十分顺利.但是,小L对数据结构的掌握实在十分渣渣.所以,小L当时卡在了二叉树. ...
- HDU5950(矩阵快速幂)
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=5950 题意:f(n) = f(n-1) + 2*f(n-2) + n^4,f(1) = a , f(2 ...
- 51nod 1126 矩阵快速幂 水
有一个序列是这样定义的:f(1) = 1, f(2) = 1, f(n) = (A * f(n - 1) + B * f(n - 2)) mod 7. 给出A,B和N,求f(n)的值. Input 输 ...
- hdu2604(递推,矩阵快速幂)
题目链接:hdu2604 这题重要的递推公式,找到公式就很easy了(这道题和hdu1757(题解)类似,只是这道题需要自己推公式) 可以直接找规律,推出递推公式,也有另一种找递推公式的方法:(PS: ...
- 矩阵乘法&矩阵快速幂&矩阵快速幂解决线性递推式
矩阵乘法,顾名思义矩阵与矩阵相乘, 两矩阵可相乘的前提:第一个矩阵的行与第二个矩阵的列相等 相乘原则: a b * A B = a*A+b*C a*c+b*D c d ...
随机推荐
- C语言使用clock进行计时
#include "stdio.h" #include "stdlib.h" #include "time.h" int main( voi ...
- WCF消息压缩
对于WCF应用来说,传输前压缩请求消息和回复消息,不但可以降低网络流量,也可以提高网络传输的性能 一.消息压缩方案 二.用于数据压缩与解压缩组件 三.用于消息压缩与解压的组件 四.用于对请求/回复消息 ...
- hive 锁表问题
报错如下: Unable to acquire IMPLICIT, EXCLUSIVE lock dms@pc_user_msg@month=201611 after 100 attempts. 显示 ...
- sublime Text3快捷键使用大全
Ctrl+D 选中光标所占的文本,继续操作则会选中下一个相同的文本.Alt+F3 选中文本按下快捷键,即可一次性选择全部的相同文本进行同时编辑.举个栗子:快速选中并更改所有相同的变量名.函数名等.Ct ...
- Android中的socket本地通讯框架
一.先分析Native层: 1.C++基类SocketListener: class SocketListener { int mSock; const char *mSocketNa ...
- shell脚本编写步骤及其常用命令和符号
1,什么是Shell Shell 是kernel的一个外壳,是一个命令解析器,负责用户与内核的交互.2,Shell脚本 Shell脚本类似于批处理,可以方便的执行大量命令.3,编写sh ...
- SQL使用视图的优缺点
视图是为了查询方便!也就是多个表的总结!但是不能对视图增删改! 在做数据库开发中使用视图的优点有: 1.视图的好处就是在你做复杂的查询逻辑时可以简化你的思考过程. 2.用视图可以隐藏一定的信息,用过滤 ...
- 安卓 ADB常见问题整理
以下都是ADB连接问题,可以通过尝试如下步骤,由简单度排序 1. 插拔下USB连接线 2. 关闭USB模式再打开 3. 执行以下命令 adb kill-server adb start-server ...
- AFNetworking 用法详解
之前一直使用ASIHttpRequest 做网络请求 ,后来新公司用AFNetWorking ,经过一段时间学习总结一下二者的优缺点: 1.AFNetWorking的优缺点 优点: 1.维护和使用者比 ...
- 更改oracle字符集
修改oracle字符集 方法/步骤 oracle数据库的字符集更改 A.oracle server 端 字符集查询 select userenv('language') from dual 其中 ...