题目链接:https://vjudge.net/problem/HDU-4565 So Easy! Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 5525    Accepted Submission(s): 1841 Problem Description A sequence Sn is defined as:Where a, b,…

/* HDU4565 So Easy! http://acm.hdu.edu.cn/showproblem.php?pid=4565 数论 快速幂 矩阵快速幂 题意:求[(a+sqrt(b))^n ]%m [ ]代表向上取证 联想到斐波那契数列,所以第一感觉是矩阵快速幂 后来发现根本不用这么麻烦,我们认为a+b*sqrt(c)中的a为实部 b为虚部,就能直接用二元乘法模拟矩阵快速幂了,因此这两种方法 是等价的.于是乎,我们就解决了精度问题. 然而即使我们得出了结果中的a+b*sqrt(c)也不能…

A sequence Sn is defined as: Where a, b, n, m are positive integers.┌x┐is the ceil of x. For example, ┌3.14┐=4. You are to calculate Sn. You, a top coder, say: So easy! 矩阵快速幂 #include<stdio.h> #include<string.h> #include<math.h> typedef…

easy 个屁啊,一点都不easy,题目就是要求公式的值,但是要求公式在最后的取模前的值向上取整.再取模,无脑的先试了高速幂 double  fmod来做,结果发现是有问题的.这题要做肯定得凑整数,凑整  题目给 a+√b 那么加上a-√b就能够了.但是这样加上后面怎么处理还得减去.想了半年也想不出来. 原来用了负数的共轭思想.还有就是题目给的b的范围 是 ((a-1)*(a-1),a*a).所以 a-√b的值的 不管多少次方 的值都是小于1的,所以对于原式子 改装成 ((a + √b) ^n+…

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4565 题解:(a+√b)^n=xn+yn*√b,(a-√b)^n=xn-yn*√b, (a+√b)^n=2*xn-(a-√b)^n,(0<=a-√b<=1),所以整数部分就是2*xn 然后再利用两个公式 (a+√b)^(n+1)=(a+√b)*(xn+yn*√b) (a-√b)^(n+1)=(a-√b)*(xn-yn*√b) 联立得到 x(n+1)=a*xn+b*yn y(n+1)=xn+a*yn…

[HDU4565]So Easy! 题面 要你求 \[ \lceil (a+\sqrt b)^n \rceil \% m \] 其中\(0<a,m<2^{15},(a-1)^2<b<a^2,0<b,n<2^{31}\) 题解 这个向上取整放在这里很丑对吧,我们化一下柿子. \[ \because (a-1)^2<b<a^2\\ \therefore a-\sqrt b \] 因为\((a-\sqrt b)^n\)是个很小的小数且\(a-\sqrt b\)与\…

转载:http://www.klogk.com/posts/hdu4565/ 这里写的非常好,看看就知道了啊. 题意很easy.a,b,n都是正整数.求 Sn=⌈(a+b√)n⌉%m,(a−1)2<b<a2 这个题目也是2008年Google Codejam Round 1A的C题. 做法事实上很easy.记(a+b√)n为An,配项 Cn=An+Bn=(a+b√)n+(a−b√)n 两项恰好共轭,所以Cn是整数. 又依据限制条件 (a−1)2<b<a2⇒0<a−b√<…

Problem Description A sequence Sn is defined as:Where a, b, n, m are positive integers.┌x┐is the ceil of x. For example, ┌3.14┐=4. You are to calculate Sn. You, a top coder, say: So easy!  Input There are several test cases, each test case in one lin…

Problem Description A sequence Sn is defined as: Where a, b, n, m are positive integers.┌x┐is the ceil of x. For example, ┌3.14┐=4. You are to calculate Sn. You, a top coder, say: So easy! Input There are several test cases, each test case in one lin…

libsvm中有进行参数调优的工具grid.py和easy.py可以使用,这些工具可以帮助我们选择更好的参数,减少自己参数选优带来的烦扰. 所需工具:libsvm.gnuplot 本机环境:Windows7(64 bit) ,Python3.5 1.相关程序的下载和安装: 1.1.下载libsvm,我用的是libsvm-3.18.zip,下载后直接解压到任意位置,我解压到C:\libsvm-3.18下. 1.2.下载python,我下的是python-3.5.msi,双击该文件安装到默认位置,我…