|
그리고 중요한것은 a의 약수이며, a/2와 나머지가 같은 k가 아니라
a²+k²=2ck
(a²+k²)/2k의 나머지는 0 이라는 것이 더 정확하군요.
1~10000까지의 수중에 저번것과 이번것의 오차는 약 450000개입니다.
a²/2k + k/2의 나머지는 0이여야 합니다.
#include <stdio.h>
#define MAX 10000
main()
{
long unsigned int a,b,c,k;
FILE *fp;
fp=fopen("pitamath.txt","w");
for(a=3;a<=MAX;a++)
{
for(k=1;k<a;k++)
{
if((a*a+k*k)%(2*k)==0)
{
c=(a*a+k*k)/(2*k);
b=c-k;
fprintf(fp,"%d, %d, %d\n",a,b,c);
}
}
}
fclose(fp);
}
무려 21줄이나 되는군요.
여기서 제일 중요한것은 a²+k²=2ck 라는 것입니다. 즉 (a²+k²)/2ck는 정수여야하며,
(a*a+k*k)%(2*k)==0를 사용해야 합니다.
제작자 : BlackFrozen |
