题意:找出满足x²+y²=k(xy+1)的(x,y)在1<=x<=y<=n范围内共有多少个
思路:看到范围,显然这应该是一道打表规律题
于是我们开始推公式:把x视为常量,(x,y1)为一个点,(x,y2)又是一个点,由韦达定理可以推出下面结论:
起始条件y1=0(这个点不算在内)则k=x²(这个是以后也是这样的,推出k恒为x²),y2=kx
然后把y2视为新的x,x视为y1进行下一轮推算,见下图:
因此:
之后就判断就行
代码用到了一个知识:
参考文章:关于lower_bound( )和upper_bound( )的常见用法_brandong-CSDN博客_lower_bound
代码:
#include <bits/stdc++.h>
using namespace std;
//(x,kx-y);k=x2,这是带入原方程得到的,所以一直是这个数,而其他的x和y是根据上一个点推出来的
//(x,x^3)(x^3,x)(x^3,x^5-x),原来是这样算的,k=x^2是定的
typedef long long ll;
const ll M = 1e18 + 1, N = 1e6 + 10;
ll n;
int t;
ll ans[N];
vector<ll> v;//放答案,我也不知道有多少个,正常这里可以二分,但是既然vector中有这个函数,为什么不用呢
int main()
{
cin >> t;
//因为第一个点是(x,x^3),所以是x*x*x,所以数组只需要开6次方
for (int x = 2; 1ll * x * x * x <= M; x++)x=1是无解的,但是从1开始会死循环
{//嘶,我就想不明白了,为啥加上=M就能过,没有+就不过?????、
ans[0] = 0, ans[1] = x;
for (int i = 2;; i++)
{
__int128 tem = 1ll * x * x;
tem = tem * ans[i - 1] - ans[i - 2];
if (tem > M)
{
break;
}
ans[i] = tem;
v.push_back(tem);
}
}
sort(v.begin(), v.end());
while (t--)
{
scanf("%lld", &n);
int cnt = upper_bound(v.begin(), v.end(), n) - v.begin() + 1;
printf("%d\n", cnt);
}
return 0;
}