HDU 4565 -- So Easy! (矩阵幂模板)
来源:互联网 发布:小鲜肉演技知乎 编辑:程序博客网 时间:2024/06/06 03:36
题目链接:hdu4565
题目要求这个得值
但是取模前有更号,所以无法直接计算,我们发现
0< a, m < 215, (a-1)2< b < a2, 0 < b, n < 231
所以 0 <| a+sqrt( b ) | < 1
可得表达式:,由二项式展开可知等号右边一坨是整数并且加的数小于一,所以等式成立
q=(a+sqrt(b))*(a-sqrt(b))=a*a-b
接下来就与loj1070相同:博客
另Sn=(a+sqrt(b))^n+(a-sqrt(b))^n。
就可以求递推式 :Sn=pS(n-1)-q*S(n-2)
即 (S2 S1)*(p 1)^(n-2)=(Sn S(n-2) )
-q 0
#include <cstdio>#include <cstring>#include <cmath>#include <iostream>using namespace std;typedef __int64 ll;#define N 3 // 这里开大了超时,原来开的32int K,mod;struct Matrix{ int r,c; ll m[N][N]; Matrix(){} Matrix(int r,int c):r(r),c(c){} Matrix operator *(const Matrix& B)//乘法 { Matrix T(r,B.c); for(int i=1;i<=T.r;i++) { for(int j=1;j<=T.c;j++) { ll tt = 0; for(int k=1;k<=c;k++) tt += m[i][k]*B.m[k][j] % mod; T.m[i][j] = tt % mod; } } return T; } Matrix Unit(int h) // 对角线矩阵 { Matrix T(h,h); memset(T.m,0,sizeof(T.m)); for(int i=1;i<=h;i++) T.m[i][i] = 1; return T; } Matrix Pow(int n) //矩阵幂 { Matrix P = *this,Res = Unit(r); while(n) { if(n&1) Res = Res*P; P = P*P; n >>= 1; } return Res; } void Print()//输出 { for(int i=1;i<=r;i++) { for(int j=1;j<=c;j++) printf("%d ",m[i][j]); printf("\n"); } }}Single;int main(){ll a,b,n;while(cin >> a >> b >> n >> mod){ll p=2*a%mod;//ll q=a*a%mod-b;Matrix cnt(2,2),ans(2,1);if(n==1){cout << p << endl;continue;}ans.m[1][1] = p*p%mod-2*q;ans.m[2][1] = p;// 操作矩阵cnt.m[1][1] = ans.m[2][1];cnt.m[1][2] = (-q +mod)%mod;cnt.m[2][1] = 1;cnt.m[2][2] = 0;ans = cnt.Pow(n-2) * ans;cout << ans.m[1][1] << endl;}return 0;}
阅读全文
0 0
- HDU 4565 -- So Easy! (矩阵幂模板)
- HDU 4565 So Easy! (矩阵快速幂)
- hdu 4565 - So Easy!(矩阵快速幂)
- hdu 4565 So Easy! 矩阵快速幂
- HDU 4565 So Easy!(矩阵快速幂)
- HDU 4565 So Easy!(矩阵快速幂)
- HDU 4565So Easy! 矩阵快速幂
- hdu 4565 So Easy!(矩阵快速幂)
- HDU - 4565 So Easy! 矩阵快速幂
- HDU 4565 so easy ! 矩阵快速幂
- hnu13150(hdu 4565)SO EASY!(矩阵快速幂)
- hdu 4565 So Easy!(矩阵+快速幂)
- hdu - 4565 - So Easy!(矩阵快速幂)
- HDU 4565 So Easy!(数学+矩阵快速幂)
- hdu 4565 So Easy!(推公式&矩阵快速幂)
- HDU 4565 So Easy!(矩阵快速幂)
- hdu 4565 So Easy!(矩阵快速幂)
- HDU 4565 so easy(矩阵快速幂)
- Java HashMap工作原理及实现
- 透明视频播放
- eclipse中配置Pydev插件出现的问题
- 百度之星资格赛1003 度度熊与邪恶大魔王
- lvm逻辑卷和设备的快照
- HDU 4565 -- So Easy! (矩阵幂模板)
- 其他对象
- Unity 中 启动页 滚动 代码实现
- 水波纹特效—Ripple
- bzoj 3749: [POI2015]Łasuchy
- Struts2知识点回顾(六)文件上传与下载
- null和“”,==和equals()的混用和区别
- 树莓派raspberry之samba
- HDOJ 1869 六度分离(floyd最短路~)