#include #include #include #include #include #include #include // 内部変数 static FILE *szpFpI; // 入力 // 内部変数 - テスト用 #ifdef D_TEST static int siRes; static FILE *szpFpA; #endif // 最大公約数 - 取得 int fGetGcd( int piVal1 // 値1 , int piVal2 // 値2 ) { int liWork; // ユークリッドの互除法 if (piVal1 < piVal2) { liWork = piVal1; piVal1 = piVal2; piVal2 = liWork; } while (piVal2 > 0) { liWork = piVal1 % piVal2; piVal1 = piVal2; piVal2 = liWork; } return piVal1; } // 実行メイン int fMain( int piTNo // テスト番号 1~ ) { char lc1Buf[1024], lc1Out[1024]; // 入力 - セット #ifdef D_TEST sprintf(lc1Buf, ".\\Test\\T%d.txt", piTNo); szpFpI = fopen(lc1Buf, "r"); sprintf(lc1Buf, ".\\Test\\A%d.txt", piTNo); szpFpA = fopen(lc1Buf, "r"); siRes = 0; #else szpFpI = stdin; #endif // 長さ - 取得 long long llLen; fgets(lc1Buf, sizeof(lc1Buf), szpFpI); sscanf(lc1Buf, "%lld", &llLen); // 1つ分の直角三角形の長さ llLen /= 4; // 上限値 int liLimit = (int)sqrt((double)llLen); // 原始ピタゴラス数の列挙 int liCnt = 0; for (int m = 2; m <= liLimit; m++) { for (int n = m - 1; n >= 1; n -= 2) { // 最大公約数 - 取得 int liGcd = fGetGcd(m, n); if (liGcd != 1) { // 互いに素ではない continue; } // 直角三角形の3辺 int liA = m * m - n * n; int liB = 2 * m * n; int liC = m * m + n * n; // 長さ合計 long long llSum = (long long)liA + (long long)liB + (long long)liC; if (llSum > llLen) { continue; } // 対象の数 liCnt++; } } // 結果 - セット sprintf(lc1Out, "%d\n", liCnt); // 結果 - 表示 #ifdef D_TEST fgets(lc1Buf, sizeof(lc1Buf), szpFpA); if (strcmp(lc1Buf, lc1Out)) { siRes = -1; } #else printf("%s", lc1Out); #endif // 残データ有無 #ifdef D_TEST lc1Buf[0] = '\0'; fgets(lc1Buf, sizeof(lc1Buf), szpFpA); if (strcmp(lc1Buf, "")) { siRes = -1; } #endif // テストファイルクローズ #ifdef D_TEST fclose(szpFpI); fclose(szpFpA); #endif // テスト結果 #ifdef D_TEST if (siRes == 0) { printf("OK %d\n", piTNo); } else { printf("NG %d\n", piTNo); } #endif return 0; } int main() { #ifdef D_TEST int i; for (i = D_TEST_SNO; i <= D_TEST_ENO; i++) { fMain(i); } #else fMain(0); #endif return 0; }