Math
题意：
问你有多少对(x,y),1<=x<=y<=n,满足(x2 + y2)%(xy+1) == 0
题解：
这种题。。。直接打表芜湖~
通过打表发现：满足情况的为(i,i * i * i),但是也有不和谐的声音出现：当x=8时，会出现两个，一个是(8,30),另一个是(8,512)，后者依然满足规律，所以前者有问题，完美继续找发现27也是，不满足的是(27,240),再往下发现有(30,112),再往下看会发现，不满足规律的情况其实是很多条链：
(8,30)(30,112)（112，418）…
(27,240)（240，）…
(64,1020)(1020, )…
我们发现这些数很多是相连的，而不想连的数彼此之间开头都是i * i * i，再通过枚举几个总结规律，每个链的开始都是a=i * i * i,后面紧跟着是(a * i * i )-pre,pre是上一组
找到规律，我们预处理出所有情况，然后排序，找到每个情况的上限值，直接二分就行
代码：
#include<bits/stdc++.h>
#define debug(a,b) printf("%s = %d\n",a,b);
typedef long double ll;
using namespace std;
//Fe~Jozky
const ll INF=0x3f3f3f3f;
inline int read(){
int s=0,w=1;
char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')w=-1;ch=getchar();}
while(ch>='0'&&ch<='9') s=s*10+ch-'0',ch=getchar();//s=(s<<3)+(s<<1)+(ch^48);
return s*w;
}
ll inf=1e18;
const int N=1e6;
long double a[7000007];
int tot=0;
void init(){
for(int i=2;i<=N;i++){
ll x=1ll*i*i*i;
a[++tot]=x;
ll t=x*i*i-i;
if(t>inf)continue;
ll prea=x;
while(t<inf){
a[++tot]=t;
t=1ll*t*i*i-prea;
prea=a[tot];
}
}
sort(a+1,a+1+tot);
}
int main()
{
init();
int t=read();
while(t--){
long double x;
cin>>x;
int id=lower_bound(a+1,a+1+tot,x)-a;
if(a[id]>x)id--;
cout<<id+1<<endl;
}
return 0;
}