反切函数的应用基础详解（核心代码两行）

反正切函数可展开成无穷级数，有如下公式
(其中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§+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/21/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
题意描述：
给出a求b+c的最小值，其中a与b和c的关系是：arctan(1/a)=arctan(1/b)+arctan(1/c)
基本思路是将公式进行转化
转化过程：令A=arctan(1/a);B=arctan(1/b);C=arctan(1/c);
tanA=1/a;tanB=1/b;tanC=1/c;
所以将公式转化为A=B+C；
tan(A)=tan(B+C) a=(bc-1)/b+c;
所以a,b,c关系又为ab+ac=bc-1；令
b=a+m;c=a+n;求出nm的值为a^2+1;b+c=2a+n+m;
所以最终问题转化为求出n+m的最小值
n+m最小值解法
因为nm=a^2+1，所以n和m之间肯定会有一个值<=a,(如果都大于两者之积肯定会大于aa+1;)所以假设m是那个小于a的数，让它从a开始依次减一，如果aa+1对m取余为0，（nm=aa+1所以n和m都是aa+1的因数）则说明m找到了，至于n=aa+1/m即可；
AC代码：

#include<stdio.h>
int main()
{
	long long  a,i;
	while(scanf("%lld", &a)!=EOF)
	{
		for(i=a;;i--)
		{
			if((a*a+1)%i==0)
			{
			printf("%lld\n", 2*a+i+(a*a+1)/i);
			break;	
			}
		}
	}
	return 0;
}