标签：right log OI dfrac ne 数学 2k operatorname 技巧
在OI比赛中，如果能够灵活地运用一些数学小技巧，是能够很好地优化计算的时间和正确性的。
既然说了是小技巧，那么这些指的都是一些技巧，一般是不会单独成题的。
光速幂
有的时候，我们要去求解一个数或者一个矩阵的若干次幂，而这个指数在一般情况下是暴力无法接受的，这个时候我们会想到使用快速幂，即：
\(a^n=\begin{cases}\left(a^2\right)^{\dfrac{n}{2}}&2\mid n\\ a\left(a^2\right)^{\dfrac{n-1}{2}}&2\nmid n\end{cases}\)
这种做法可以做到优秀的 \(O(logn)\) 的复杂度。
但是，如果我们继续提高 \(n\) 的上限，long long存不下了呢？
栗子：Luogu P1397 [NOI2013] 矩阵游戏
这道题目有一个特别和善的数据范围： \(n \leqslant 10^{1000000}\)
似乎要写高精了，如果还是仿照快速幂的方式写，不难发现，由于高精度是按位处理，我们会发现，我们的复杂度会被卡到 \(O(log^2n)\) ，这个复杂度是我们无法接受的。
但我们想想，我们存的高精度是十进制的，如果我们用十进制来做快速幂，会怎么样？
我们考虑这样的形式：
\(a^n=\left(a^{10}\right)^{\left\lfloor\dfrac{n}{10}\right\rfloor}*a^{n \bmod 10}\)
我们只需要先预处理出 \(a^2\) 到 \(a^9\) 即可。
这样，我们规避了高精度除单精的问题，可以将时间优化到 \(O(logn)\) 了。
不难发现，我们平时写的快速幂就是光速幂的一种特殊形式，也就是对于二进制的处理方法
那么，是不是对于任意k进制，都可以用光速幂呢？
其实在形式上都是一样的，就是把上面的 \(10\) 换成 \(k\) 就可以了。
闲话（不知道正不正确，欢迎大佬指正）：那么对于不同的 \(n\) ，最优的 \(k\) 会是多少呢？
由于我们需要进行预处理，首先会有一个 \(k\) 的运算，如果不考虑做除法的时间，后面的运算是 \(log_kn\) 的复杂度
那么它的总复杂度是 \(k+log_kn\) ，这个函数是先减后增的。
我们来找它的极小点，对 \(k\) 求导
\((k+\operatorname{log}_kn)'=(k)'+(\operatorname{log}_kn)'=1+\left(\dfrac{1}{\operatorname{log}_nk}\right)'\)
\(\left(\dfrac{1}{\operatorname{log}_nk}\right)'=-\dfrac{1}{\operatorname{log}_n^2k}*(\operatorname{log}_nk)'=-\dfrac{1}{\operatorname{log}_n^2k}*\operatorname{log}_ne(\operatorname{ln}k)'=-\dfrac{\operatorname{log}_ne}{k\operatorname{log}_n^2k}\)
\((k+\operatorname{log}_kn)'=1-\dfrac{\operatorname{log}_ne}{k\operatorname{log}_n^2k}\)
当导数取 \(0\) 时：\(0=1-\dfrac{\operatorname{log}_ne}{k\operatorname{log}_n^2k}\)
则 \(\dfrac{\operatorname{log}_ne}{k\operatorname{log}_n^2k}=1\)
所以 \(\operatorname{log}_ne=\operatorname{log}_n^2k\)
所以 \(e=k^{\operatorname{log}_nk}\)
\(e^{\operatorname{ln}n}=k^{\operatorname{ln}nlog_nk}\)
\(n=k^{\operatorname{ln}k}\)
对于上面的题目，大概 \(k\) 是在 \(10^5\) 左右
标签：right,log,OI,dfrac,ne,数学,2k,operatorname,技巧 来源： https://www.cnblogs.com/Xun-Xiaoyao/p/16624066.html
本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享； 2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关； 3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关； 4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除； 5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。