ARCTAN - Use of Function Arctan

题目地址请点击——

---

ARCTAN

---

Description

相信大家都看得懂英文吧，我就不翻译了吧…………

It is easy to know that

arctan(1/2)+arctan(1/3)=arctan(1)

$$arctan(1/2)+arctan(1/3)=arctan(1)$$
Given number $A$, write a program to calculate the minimum sum $B+C$.
$A$，$B$ and $C$ are all positive integers and satisfy the equation below:
$$arctan(1/A)=arctan(1/B)+arctan(1/C)$$

---

Input

The first line contains a integer number $T$(about $1000$).
$T$ lines follow,each contains a single integer $A$，$1\leq A\leq 60000$.

---

Output

$T$ lines，each contains a single integer which denotes to the minimum sum $B+C$.

---

Sample Input

1

---

Sample Output

5

---

Solution

解法一：附上扫雷神犇Cai的坑爹解法。
解法二：
由

tan(α+β)=1A,tanβ=1B,tanα=1C

$$tan(α+β)=\frac{1}{A},tanβ=\frac{1}{B},tanα=\frac{1}{C}$$
得 $$A=\frac{BC-1}{B+C}$$
设 $B+C=k$，则 $$BC=Ak+1$$ ，
则 $B、C$ 为一元二次方程 $$x^2-kx+Ak+1=0$$ 的两个根。
那么 $$k^2-4Ak-4$$ 为完全平方数。
设 $k^2-4Ak-4=s^2$，
$$\begin{align*} k&=\frac{4A+\sqrt{16A^2+16+4s^2}}{2}\\ &=2A+\sqrt{4A^2+s^2+4} \end{align*}$$
则 $\sqrt{4A^2+s^2+4}$ 为完全平方数。
设 $4A^2+s^2+4=r^2$，则 $4A^2+4=(r+s)(r-s)$
枚举 $4A^2+4$ 的约数即可。

---

Code

因为有代码长度限制，所以读者在复制代码的时候请注意删除缩进以及空格…………

#include<iostream>
#include<cmath>

typedef unsigned long long L;
using namespace std;

L t,i,T,a,b,k;

main(){
    cin>>t;
    while(t--){
        cin>>a;
        T=a*a*4+4;
        b=sqrt(T);
        i=b+1;
        while(--i)
            if(!(T%i)&&!((i+T/i)&1)){
                k=(T/i-i)/2;
                break;
            }
        cout<<2*a+L(sqrt(T+k*k)+0.5)<<endl;
    }
    return 0;
}