I was solving this question Pythagorean Theorem II
and saw some optimized solution based on this relation seems
Saw some solved optimized codes using this relations seems gcd(a,b,c)=1
int n,m=0;
cin>>n;
for(int i=1;i<=n;i++)
{
for(int j=1;(i*i+j*j<=n&&i>j);j++)
{
if((i-j)%2==1)
if(gcd(i,j)==1)
m+=n/(i*i+j*j);
}
}
cout<<m;
I know a,b,c are relatively prime but is it enough to know that ?
Can somebody tell me proof of this approach ?