【HDOJ5950】Recursive sequence(矩阵乘法,快速幂)
题意:f[1]=a,f[2]=b,f[i]=2f[i-2]+f[i-1]+i^4(i>=3),多组询问求f[n]对2147493647取模
N,a,b < 2^31
思路:重点在于i^4的处理
对于i转移矩阵中可以记录下它的0,1,2,3,4次项
i的幂又可以由i-1的幂运算得出,最后推出的系数是二项式展开的系数
试试新的矩乘模板
- #include<cstdio>
- #include<cstring>
- #include<string>
- #include<cmath>
- #include<iostream>
- #include<algorithm>
- #include<map>
- #include<set>
- #include<queue>
- #include<vector>
- using namespace std;
- typedef long long ll;
- typedef unsigned int uint;
- typedef unsigned long long ull;
- typedef pair<int,int> PII;
- typedef vector<int> VI;
- #define fi first
- #define se second
- #define MP make_pair
- #define N 2100000
- #define MOD 2147493647
- #define eps 1e-8
- #define pi acos(-1)
- const int MAXN=;
- int read()
- {
- int v=,f=;
- char c=getchar();
- while(c<||<c) {if(c=='-') f=-; c=getchar();}
- while(<=c&&c<=) v=(v<<)+v+v+c-,c=getchar();
- return v*f;
- }
- struct matrix //矩阵类
- {
- int n,m;
- ll data[MAXN][MAXN];
- };
- matrix ma,mb;
- ll a,b,c,d,p,n;
- matrix matrixMul(matrix a, matrix b)
- {
- matrix re;
- if(a.m!=b.n)
- {
- printf("error\n");
- return re;
- }
- memset(re.data,,sizeof(re.data));
- re.n = a.n; re.m = b.m;
- for(int i = ; i <= a.n; i++)
- {
- for(int j = ; j <= a.m; j++)
- {
- if(a.data[i][j] == ) continue;
- for(int k = ; k <= b.m; k++)
- {
- re.data[i][k] += (a.data[i][j] % MOD * b.data[j][k] % MOD) % MOD;
- re.data[i][k] %= MOD;
- }
- }
- }
- return re;
- }
- matrix matrixPow(matrix a,int b)
- {
- matrix re;
- if(a.n!=a.m)
- {
- printf("error2\n");
- return re;
- }
- re.n=re.m=a.n;
- memset(re.data,,sizeof(re.data));
- for(int i=;i<=re.n;i++) re.data[i][i]=;
- while(b)
- {
- if(b&) re=matrixMul(re,a);
- a=matrixMul(a,a);
- b>>=;
- }
- return re;
- }
- void inputMat(int n,int m,matrix &a,ll *b)
- {
- a.n = n; a.m = m;
- for(int i = ; i <= n; i++)
- for(int j = ; j <= m; j++)
- a.data[i][j] = *(b + (i - ) * m + (j - ));
- }
- void init(){
- ll pt[][] = {b,a,,,,,};
- inputMat(,,ma,*pt);
- ll pt2[][] = {,,,,,,,
- ,,,,,,,
- ,,,,,,,
- ,,,,,,,
- ,,,,,,,
- ,,,,,,,
- ,,,,,,,};
- inputMat(,,mb,*pt2);
- }
- int main()
- {
- int cas;
- scanf("%d",&cas);
- while(cas--)
- {
- scanf("%I64d%I64d%I64d",&n,&a,&b);
- int i=;
- a%=MOD;
- b%=MOD;
- if(n == )
- printf("%I64d\n",a);
- else if(n == )
- printf("%I64d\n",b);
- else
- {
- init();
- ma=matrixMul(ma,matrixPow(mb,n-));
- }
- printf("%I64d\n",ma.data[][]);
- }
- return ;
- }
【HDOJ5950】Recursive sequence(矩阵乘法,快速幂)的更多相关文章
- Qbxt 模拟赛 Day4 T2 gcd(矩阵乘法快速幂)
/* 矩阵乘法+快速幂. 一开始迷之题意.. 这个gcd有个规律. a b b c=a*x+b(x为常数). 然后要使b+c最小的话. 那x就等于1咯. 那么问题转化为求 a b b a+b 就是斐波 ...
- 洛谷 P4910 帕秋莉的手环 矩阵乘法+快速幂详解
矩阵快速幂解法: 这是一个类似斐波那契数列的矩乘快速幂,所以推荐大家先做一下下列题目:(会了,差不多就是多倍经验题了) 注:如果你不会矩阵乘法,可以了解一下P3390的题解 P1939 [模板]矩阵加 ...
- ACM学习历程—HDU5667 Sequence(数论 && 矩阵乘法 && 快速幂)
http://acm.hdu.edu.cn/showproblem.php?pid=5667 这题的关键是处理指数,因为最后结果是a^t这种的,主要是如何计算t. 发现t是一个递推式,t(n) = c ...
- 矩阵乘法快速幂 codevs 1732 Fibonacci数列 2
1732 Fibonacci数列 2 时间限制: 1 s 空间限制: 128000 KB 题目等级 : 钻石 Diamond 题解 查看运行结果 题目描述 Description 在“ ...
- codevs1281 矩阵乘法 快速幂 !!!手写乘法取模!!! 练习struct的构造函数和成员函数
对于这道题目以及我的快速幂以及我的一节半晚自习我表示无力吐槽,, 首先矩阵乘法和快速幂没必要太多说吧,,嗯没必要,,我相信没必要,,实在做不出来写两个矩阵手推一下也就能理解矩阵的顺序了,要格外注意一些 ...
- [vijos1725&bzoj2875]随机数生成器<矩阵乘法&快速幂&快速乘>
题目链接:https://vijos.org/p/1725 http://www.lydsy.com/JudgeOnline/problem.php?id=2875 这题是前几年的noi的题,时间比较 ...
- [codevs]1250斐波那契数列<矩阵乘法&快速幂>
题目描述 Description 定义:f0=f1=1, fn=fn-1+fn-2(n>=2).{fi}称为Fibonacci数列. 输入n,求fn mod q.其中1<=q<=30 ...
- 【BZOJ-1009】GT考试 KMP+DP+矩阵乘法+快速幂
1009: [HNOI2008]GT考试 Time Limit: 1 Sec Memory Limit: 162 MBSubmit: 2745 Solved: 1694[Submit][Statu ...
- 矩阵乘法快速幂 codevs 1574 广义斐波那契数列
codevs 1574 广义斐波那契数列 时间限制: 1 s 空间限制: 256000 KB 题目等级 : 钻石 Diamond 题目描述 Description 广义的斐波那契数列是指形如 ...
- BZOJ-1875 HH去散步 DP+矩阵乘法快速幂
1875: [SDOI2009]HH去散步 Time Limit: 20 Sec Memory Limit: 64 MB Submit: 1196 Solved: 553 [Submit][Statu ...
随机推荐
- dfs染色法判定二分图
#include<iostream> #include<cstring> using namespace std; ][],color[],n; int dfs(int x,i ...
- ValueError: option names {'--alluredir'} already added 报错
运行测试用例 import pytest from WXP2P_2.test_data2.login_case import logindata_error1,logindata_error2,log ...
- Robot Framework(十) 执行测试用例——测试执行
3.2测试执行 本节描述如何执行从解析的测试数据创建的测试套件结构,如何在失败后继续执行测试用例,以及如何正常停止整个测试执行. 3.2.1执行流程 执行套房和测试 设置和拆卸 执行顺序 3.2.2继 ...
- 一个小笔记(5):A*算法
A-Star算法是一种静态路网中求解最短路径最有效的直接搜索方法其实百科有 http://baike.baidu.com/link?url=CvmkWQIAmztYgMq3Nk1WyWkDiC0koV ...
- HDU-1217-Arbitrage(SPFA)
这题和以往的求最短路的题目略微有点不一样,以往求的都是最小的,这题求的是大的,而且还是乘法. 我们求的时候初始化的时候就要进行相反的初始化了,把它们初始化为0,然后比较大的就更新. 因为这题的点少边多 ...
- mysqldump导出备份数据库报Table ‘performance_schema.session_variables‘ doesn‘t exist
今天在bash进行本地数据库往云端数据库导数据的时候,在本地导出.sql文件这第一步就出现了错误问题,导出sql文件的命令: mysqldump -u 用户名 -p 数据库名 > xxx.sql ...
- this version of SLF4J requires log4j version 1.2.12 or later.
org.springframework.beans.factory.BeanCreationException: Error creating bean with name 'sessionFacto ...
- Python中的socket网络编程(TCP/IP,UDP)讲解
在网络编程中的一个基本组件就是套接字(socket).套接字基本上是两个端点的程序之间的"信息通道".程序可能分布在不同的计算机上,通过套接字互相发送信息.套接字包括两个:服务器套 ...
- Android 8.0 adb shell dumpsys activity activities | findstr mFocusedActivity 获取当前的 activity 显示空的
adb shell dumpsys activity activities | findstr mFocusedActivity Android 7.0 现象: Android 8.0 现象: 改用: ...
- Oracle审计相关对象的迁移
目录 创建审计用的表空间 在线迁移 查询结果 在日常的数据库维护中,经常出现因为数据库登录审计的功能启动,导致system表空间被用满.从而出现异常,一般建议把aud$相关对象迁移到其他表空间,从而避 ...