A,B,CA,B,CA,B,C都是整数,给出一个公式arctanA+arctanB=arctan(A+B1−A×B)arctanA+arctanB=arctan(\frac{A+B}{1-A\times B})arctanA+arctanB=arctan(1−A×BA+B)
现在给你arctan1A=arctan1B+arctan1Carctan\frac1 A=arctan\frac1 B+arctan\frac1CarctanA1=arctanB1+arctanC1输入AAA,问B+CB+CB+C的最小解
- 根据公式很容易能够得到arctan1A=arctanB+C1−B×Carctan\frac1A=arctan\frac{B+C}{1-B\times C}arctanA1=arctan1−B×CB+C根据函数单调性容易得到AB+AC+BC=1AB+AC+BC=1AB+AC+BC=1现在已知AAA,我们要求B+CB+CB+C的最小值
- 如果想提取出AAA通过不等式得到最小值是错误的,因为在那样求出的最小值不是整数意义下的最小值,取等只有在B=CB=CB=C时候才成立,这也是高中数学的一个容易陷入的误区
- 那么我们怎么找最小值呢?我们可以做这样一个变形(B−A)×(C−A)=1+A2(B-A)\times(C-A)=1+A^2(B−A)×(C−A)=1+A2那么如果设m=B−A,n=C−Am=B-A,n=C-Am=B−A,n=C−A,那么有mn=1+A2mn=1+A^2mn=1+A2,因为都是整数,所以我们可以枚举所有1+A21+A^21+A2的因数寻找B+CB+CB+C的最小值
#include <iostream>
#include <cstring>
#include <algorithm>
#include <queue>
#include <stack>
#include <vector>
#include <cmath>
#include <cstdio>
#include <map>
#include <iomanip>
using namespace std;
typedef long long ll;
const int MAXN = 2e5 + 100;
const int INF = 0x3f3f3f3f;
int a[MAXN];
typedef long long ll;
int main(){
#ifdef LOCAL
freopen("input.txt", "r", stdin);
freopen("output.txt", "w", stdout);
#endif
ios::sync_with_stdio(false);
cin.tie(0);
cout.tie(0);
ll a;
cin >> a;
ll ans = 0x7fffffff;
ll tmp = 1 + a * a;
for(ll m=1;m<=a;m++){
if(tmp % m == 0){
ll n = tmp / m;
ans = min(ans, m + a + n + a);
}
}
cout << ans;
return 0;
}