【NOIOL】【数论 / 数值】正反切函数的应用

题目描述

反正切函数可展开成无穷级数，有如下公式

(其中)

使用反正切函数计算是一种常用的方法。例如，最简单的计算的方法：



然而，这种方法的效率很低，但我们可以根据角度和的正切函数公式：



通过简单的变换得到：



利用这个公式，令，则，有



使用和的反正切来计算，速度就快多了。
我们将该公式写成如下形式



其中和均为正整数。
我们的问题是：对于每一个给定的()，求的值。我们保证对于任意的都存在整数解。如果有多个解，要求你给出最小的解。
　

输入描述

输入文件中只有一个正整数，其中。

输出描述

输出文件中只有一个整数，为 的值。  

样例输入

1

样例输出

5

题目来源

NOI 2001

---

将公式化简，极为求a=(bc-1)/(b+c)，可是还是没有什么头绪。

再次变换==> c=(ab+1)/(b-a)==>c=a+(a^2+1)/(b-a)

                        没头绪……          perfect！

这样我们只须枚举(b-a)，sqrt(60000^2+1)≈60000，不会超时。

#include<stdio.h> #include<math.h> #include<limits.h> long long a,ans=INT_MAX,t,i,b,c; int main() { scanf("%lld",&a); t=a*a+1; for (i=1;i<=sqrt(t);++i) if (t%i==0) { b=i+a; c=a+t/i; if (b+c<ans) ans=b+c; } printf("%lld/n",ans); return 0; }