原题: http://poj.org/problem?id=1183
题目:
反正切函数的应用
Time Limit: 1000MS Memory Limit: 10000K
Total Submissions: 15009 Accepted: 5397
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
思路:
由arctan(1/a)=arctan(1/b)+arctan(1/c) 和公式(4)可以得到:
1/a=((1/b)+(1/c)) / (1-(1/b)*(1/c))=(b+c) / (b*c-1)
令b+c=x,则有c=x-b,代入上式得到:x=(b*b+1)/(b-a)
再令b-a=y,则b=y+a,代入上式有:x=((y+a)*(y+a)+1)/y=y+(a*a+1)/y+2*a
由双勾函数的特点我们知道,在sprt((a*a+1)/1)处取得最小值。
所以我们只需要从最小值处向两边寻找满足条件的最小整数解即可。
这里由于最大的a是6w,3.6*10^9会int溢出导致结果错误,usigned int却不会,而long long int好像并不能用sqrt函数。
因为对于任何的双勾曲线,在a*a+1处一定有整数解,所以我们定义的INF同样需要usigned int类型。
所以这题我们全用用usigned int类型。
代码:
#include"cstdio"
#include"math.h"
#include"iostream"
const int INF = 3600000002;
using namespace std;
typedef unsigned int uint;
int main()
{
uint a;
while(scanf("%I64d",&a)!=EOF)
{
uint l=floor(sqrt(a*a+1));
uint r=ceil(sqrt(a*a+1));
uint ansl=INF;
uint ansr=INF;
while(l>0)
{
if((a*a+1)%l==0)
{
ansl=l+(a*a+1)/l+2*a;
break;
}
l--;
}
while(r<=(a*a+1))
{
if((a*a+1)%r==0)
{
ansr=r+(a*a+1)/r+2*a;
break;
}
r++;
}
uint ans=min(ansl,ansr);
printf("%d\n",ans);
}
return 0;
}