反正切函数的应用

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

(其中0 <= x <= 1) 公式(1)

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

PI=4arctan(1)=4(1-1/3+1/5-1/7+1/9-1/11+...) 公式(2)

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

tan(a+b)=[tan(a)+tan(b)]/[1-tan(a)*tan(b)] 公式(3)

通过简单的变换得到：

arctan(p)+arctan(q)=arctan[(p+q)/(1-pq)] 公式(4)

利用这个公式，令p=1/2,q=1/3，则(p+q)/(1-pq)=1，有

arctan(1/2)+arctan(1/3)=arctan[(1/2+1/3)/(1-1/2*1/3)]=arctan(1)

使用1/2和1/3的反正切来计算arctan(1)，速度就快多了。
我们将公式(4)写成如下形式

arctan(1/a)=arctan(1/b)+arctan(1/c)

其中a,b和c均为正整数。

我们的问题是：对于每一个给定的a（1 <= a <= 60000），求b＋c的值。我们保证对于任意的a都存在整数解。如果有多个解，要求你给出b+c最小的解。


Input
输入文件中只有一个正整数a，其中 1 <= a <= 60000。

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

Sample Input
1

Sample Output
5

Source
Noi 01

=================================================

// 首先用数学方法判断一下定义域和值域
//本题题意就是：[(b*c) - 1]/(b + c) = a
//a是常数，求b + c或其最大值
//程序保证有结果，其中a(1 <= a <= 60000)
//有正切函数的性质，我们可以轻松地判断出：b和c的取值范围是[1,+∞)
//题目保证有结果，也就是(ab + 1)/(b - a)一定有整除的时候
//人类不能做的，机器能做：无限判断，直到出现结果
//比较不好的方法，还是希望能从纯数学上解决比较快……

#include "stdafx.h" #include <iostream> using namespace std; int main(){ unsigned int a,b,c; unsigned int a_By_b,b_Less_a; cout<<"输入a ，a取值[1,600000]"<<endl; cin>>a; if (a >= 1 && a <= 60000) { for (b = a + 1;;b++) { b_Less_a = b - a; a_By_b = a*b; if(b_Less_a > a_By_b + 1) break; if (!((a_By_b + 1)%b_Less_a)) { c = (a_By_b + 1)/(b_Less_a); cout<<(b + c)<<endl; break; } } } system("pause"); return 0; }

==================================================

别人的解法：

#include "stdafx.h" #include <iostream> using namespace std; int main() { int a; cout<<"-----a-----"<<endl; cin>>a; unsigned int b, den, num, c, result=0; for(b=a+1; ; ++b){ num=1+a*b; den=b-a; if(b*den>num) break; if(0==num%den){ c=num/den; result=b+c; } } cout<<result<<endl; system("pause"); return 0; }

根据两个结果的比较，我觉得他是对的。但是不是很明白那个红色标出来的判断条件是怎么来的。估计他是根据那个判断条件找到的a + b的最小值。我没有那个判断条件，找出来的不是最小值。

还是得研究一下再……