POJ 1305 Fermat vs. Pythagoras （构造 毕达哥拉斯三元组）

题目链接：POJ 1305 Fermat vs. Pythagoras

勾股数a,b,c。求在a,b,c小于n的情况下，a,b,c互质的三元组的个数，和n中没有出现在勾股数中的数的个数。

思路：构造 毕达哥拉斯三元组 （一个三元组（a，b，c),其中a，b，c无公因数，且满足a² +b² =c²。），a=x^2-y^2,b=2*x*y,c=x^2+y^2.在,a,b互质,a,b奇偶性不同条件下成立。
所有的勾股数都可以通过三元组乘上系数得到。

AC代码：

#include<stdio.h>
#include<string.h>
#include<math.h>
bool vis[1000010];
int gcd(int a,int b)
{
	return b? gcd(b,a%b):a;
}
int main()
{
	int x,y,z,k;
	int n,i,j,w;
	while(scanf("%d",&n)!=EOF)
	{
		int sum=0,count=0;
		memset(vis,false,sizeof vis);
		for(i=1;i<=n;i++)//枚举i构造毕达哥拉斯三元组
		{
			for(j=1;j<i && j*j+i*i<=n;j++)
			{
				if((i%2)^(j%2) && gcd(i,j)==1)
				{
					count++;
					for(k=1;k*(i*i+j*j)<=n;k++)
						vis[k*(i*i-j*j)]=vis[k*2*i*j]=vis[k*(j*j+i*i)]=true;
				}
			}
			if(!vis[i]) sum++;
		}
		printf("%d %d\n",count,sum);
	}
	return 0;
}