经过一系列的复杂计算及推演最终得出一个结论:任何x,y都会满足一个式子:
对任意K>=2,a[i] = pow(k,2) * a[i-1] - a[i-2],其中a[0]=0,a[1]=k。
那么我们只要枚举其中较大的那个数,打表后二分即可,由于a[2]=k^3,因此我们只需要枚举到1e6,时间复杂度1e6+Tlogn
#include<bits/stdc++.h>
using namespace std;
#define int long long
#define endl '\n'
vector<int> g;
const int maxn = 1e6;
int a[maxn];
signed main(){
int t;
cin >> t;
for (int k = 2; k <= 1e6; ++k) {
a[0] = 0;
a[1] = k;
for (int i = 2; i <= 60 ; ++i) {
a[i] = pow(k,2) * a[i-1] - a[i-2];
if(a[i]>1e18||a[i]<0) break;
g.push_back(a[i]);
}
}
sort(g.begin(),g.end());
g.erase(unique(g.begin(),g.end()),g.end());
while(t--){
int n;
cin >> n;
cout << upper_bound(g.begin(),g.end(),n)-g.begin()+1<<endl;
}
}