HDU 2256 Problem of Precision 数论矩阵快速幂
2018-06-17 23:38:53来源:未知 阅读 ()
题目要求求出(√2+√3)2n的整数部分再mod 1024。
(√2+√3)2n=(5+2√6)n
如果直接计算,用double存值,当n很大的时候,精度损失会变大,无法得到想要的结果。
我们发现(5+2√6)n+(5-2√6)n是一个整数(2√6的奇数次幂总会正负抵消掉),并且(5-2√6)n是小于1的。所以我们就只需要求出S(n)-1即可。令
An=(5+2√6)n; Bn=(5-2√6)n.
Sn=An+Bn Sn为整数。
Sn*((5+2√6)+(5-2√6))=Sn*10
Sn*10=(5+2√6)n+1+(5-2√6)n+1+(5+2√6)n-1+(5-2√6)n-1
Sn*10=Sn+1+Sn-1
递推式:Sn=10*Sn-1-Sn-2
然后转化为矩阵快速幂求Sn
#include <iostream> #include <cstdio> #include <cstring> using namespace std; const int Mod=1024; const int N=2; struct Mat { int mat[N][N]; }a; Mat Multiply(Mat a, Mat b) { Mat c; memset(c.mat, 0, sizeof(c.mat)); for(int k = 0; k < 2; ++k) for(int i = 0; i < 2; ++i) if(a.mat[i][k]) for(int j = 0; j < 2; ++j) if(b.mat[k][j]) c.mat[i][j] = (c.mat[i][j] +a.mat[i][k] * b.mat[k][j])%Mod; return c; } Mat QuickPower(Mat a, int k) { Mat c; memset(c.mat,0,sizeof(c.mat)); for(int i = 0; i < 2; ++i) c.mat[i][i]=1; for(; k; k >>= 1) { if(k&1) c = Multiply(c,a); a = Multiply(a,a); } return c; } void InitMat(Mat &A) { A.mat[0][0]=10; A.mat[0][1]=-1; A.mat[1][0]=1; A.mat[1][1]=0; } int main() { //freopen("in.txt","r",stdin); int t; scanf("%d",&t); while(t--) { int n; scanf("%d",&n); if(n==1) printf("9\n"); else if(n==2) printf("97\n"); else { InitMat(a); a=QuickPower(a,n-2); int ans=(a.mat[0][0]*98+a.mat[0][1]*10-1)%1024; //我们求的是S[n]-1 while(ans<0) ans+=1024; printf("%d\n",ans); } } return 0; }
标签:
版权申明:本站文章部分自网络,如有侵权,请联系:west999com@outlook.com
特别注意:本站所有转载文章言论不代表本站观点,本站所提供的摄影照片,插画,设计作品,如需使用,请与原作者联系,版权归原作者所有
- HDU-2955-Robberies(0-1背包) 2020-03-30
- hdu1455 拼木棍(经典dfs) 2020-02-29
- anniversary party_hdu1520 2020-02-16
- hdu1062 text reverse 2020-01-27
- hdu4841 2020-01-26
IDC资讯: 主机资讯 注册资讯 托管资讯 vps资讯 网站建设
网站运营: 建站经验 策划盈利 搜索优化 网站推广 免费资源
网络编程: Asp.Net编程 Asp编程 Php编程 Xml编程 Access Mssql Mysql 其它
服务器技术: Web服务器 Ftp服务器 Mail服务器 Dns服务器 安全防护
软件技巧: 其它软件 Word Excel Powerpoint Ghost Vista QQ空间 QQ FlashGet 迅雷
网页制作: FrontPages Dreamweaver Javascript css photoshop fireworks Flash