#include<iostream>
#include<algorithm>
#include<cmath>
using namespace std;
int n,m;
const int N=2500010;
struct Sum{
int s,c,d;
bool operator <(const Sum&t)const{
if(s!=t.s)return s<t.s;
if(c!=t.c)return c<t.c;
return d<t.d;
}
}sum[N];
int main(){
cin>>n;
for(int c=0;c*c<=n;c++){
for(int d=c;d*d+c*c<=n;d++){
sum[m++]={c*c+d*d,c,d};
}
}
sort(sum,sum+m);
for(int a=0;a*a<n;a++){
for(int b=a;b*b+a*a<n;++b){
int t=n-a*a-b*b;
int l=0,r=m-1;
while(l<r){
int mid=(l+r)>>1;
if(sum[mid].s>=t)r=mid;
else l=mid+1;
}
if(sum[l].s==t){
cout<<a<<" "<<b<<" "<<sum[l].c<<" "<<sum[l].d;
return 0;
}
}
}
return 0;
}
小单刷题笔记之四平方和二分解法(空间换时间)
最新推荐文章于 2022-10-27 20:40:53 发布
这段代码的主要目的是找出所有使得a^2 + b^2 = c^2的整数三元组(a, b, c),其中a, b, c的平方和小于给定的整数n。它通过遍历所有可能的c和d值来生成所有可能的平方和,然后使用二分查找在排序后的平方和数组中找到匹配项。
摘要由CSDN通过智能技术生成