POJ1183 反正切函数的应用 数学 C语言
题目：http://poj.org/problem?id=1183
思路：知道要推导公式，根据公式(4)的两种形式推来推去，似乎没有太大结果。看了Discuss里大牛的推导，看懂之后就很简单了，关键就是b = a+m, c = a+n .
借用大牛的推导……
1/a = (1/b + 1/c)/ (1 - 1/(b*c)) => bc-1 = a(b+c) assume b=a+m and c=a+n (b and c is always bigger than a) (a+m)(a+n)-1=a(a+m+a+n) => a*a+a*n+a*m+m*n-1=2*a*a+m*a+n*a => m*n=a*a+1 and then for(m=a;m>=1;m--) if((a*a+1)%m==0) break; n=(a*a+1)/m
提交情况： WA 1次， AC 1次
总结：看a的范围觉得int足够了，结果WA。为了保险起见以后还都是用I64吧
AC code：
View Code
1 #include <stdio.h>
2 #include <stdlib.h>
3 #include <string.h>
4
5 #define I64 __int64
6
7 int main() {
8 I64 a;
9 scanf("%I64d", &a);
10 I64 m, n;
11 for(I64 i = a; ; --i) {
12 if(!((a * a + 1) % i)) {
13 m = i;
14 n = (a * a + 1) / m;
15 break;
16 }
17 }
18 printf("%I64d\n", 2 * a + m + n);
19 return 0;
20 }
