#include #include #include #include #include #include #include // 内部定数 #define D_PRM_MAX 100000 // 最大素数数 #define D_DIV_MAX 100000 // 最大約数数 // 内部構造体 - 素数情報 typedef struct Prm { int miVal; // 値 int miCnt; // 個数 } Prm; // 内部変数 static FILE *szpFpI; // 入力 static Prm sz1Prm[D_PRM_MAX]; // 素数 static int siPCnt; // 素数数 static int si1Div[D_DIV_MAX]; // 約数 static int siDCnt; // 約数数 // 内部変数 - テスト用 #ifdef D_TEST static int siRes; static FILE *szpFpA; static int siTNo; #endif // 出力 int fOut( char *pcpLine // 1行 ) { char lc1Buf[1024]; #ifdef D_TEST fgets(lc1Buf, sizeof(lc1Buf), szpFpA); if (strcmp(lc1Buf, pcpLine)) { siRes = -1; } #else printf("%s", pcpLine); #endif return 0; } // ソート関数 - int昇順 int fSortFnc( const void *pzpVal1 // 値1 , const void *pzpVal2 // 値2 ) { int *lipVal1 = (int *)pzpVal1; int *lipVal2 = (int *)pzpVal2; // int昇順 if (*lipVal1 > *lipVal2) { return 1; } else if (*lipVal1 < *lipVal2) { return -1; } return 0; } // 素因数分解 int fDivPrm( int piVal // 対象値 ) { int i; // 初期化 memset(sz1Prm, 0, sizeof(sz1Prm)); siPCnt = 0; // 上限 int liLimit = (int)sqrt((double)piVal); // 素因数分解 for (i = 2; i <= liLimit; i++) { // 分解終了 if (piVal == 1) { break; } // 分解 while (piVal % i == 0) { piVal /= i; sz1Prm[siPCnt].miCnt++; } // 分解あり if (sz1Prm[siPCnt].miCnt > 0) { sz1Prm[siPCnt].miVal = i; siPCnt++; } } // 上限を超える素数 if (piVal != 1) { sz1Prm[siPCnt].miCnt = 1; sz1Prm[siPCnt].miVal = piVal; siPCnt++; } return 0; } // 約数リスト - 作成 int fMakeDivK( int piPNo // 素数 0~ , int piVal // 値 ) { int i; // 選択中 if (piPNo < siPCnt) { // 個数 for (i = 0; i <= sz1Prm[piPNo].miCnt; i++) { fMakeDivK(piPNo + 1, piVal); piVal *= sz1Prm[piPNo].miVal; } return 0; } // 約数 - 追加 si1Div[siDCnt] = piVal; siDCnt++; return 0; } int fMakeDiv( int piVal // 対象値 ) { // 初期化 siDCnt = 0; // 素因数分解 fDivPrm(piVal); // 約数リスト - 作成 fMakeDivK(0, 1); return 0; } // 実行メイン int fMain( ) { int i; char lc1Buf[1024]; // A, B - 取得 int liA, liB; fgets(lc1Buf, sizeof(lc1Buf), szpFpI); sscanf(lc1Buf, "%d%d", &liA, &liB); // 約数リスト - 作成 fMakeDiv(liA + liB); // 約数リスト - ソート qsort(si1Div, siDCnt, sizeof(int), fSortFnc); // 約数 - チェック for (i = 0; i < siDCnt; i++) { if (si1Div[i] == liA) { continue; } if (si1Div[i] == liB) { continue; } if ((si1Div[i] + liA) % liB != 0) { continue; } if ((si1Div[i] + liB) % liA != 0) { continue; } return si1Div[i]; } return -1; } // 1回実行 int fOne( ) { int liRet; char lc1Buf[1024]; // 入力 - セット #ifdef D_TEST sprintf(lc1Buf, ".\\Test\\T%d.txt", siTNo); szpFpI = fopen(lc1Buf, "r"); sprintf(lc1Buf, ".\\Test\\A%d.txt", siTNo); szpFpA = fopen(lc1Buf, "r"); siRes = 0; #else szpFpI = stdin; #endif // 実行メイン liRet = fMain(); // 出力 sprintf(lc1Buf, "%d\n", liRet); fOut(lc1Buf); // 残データ有無 #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", siTNo); } else { printf("NG %d\n", siTNo); } #endif return 0; } // プログラム開始 int main() { #ifdef D_TEST int i; for (i = D_TEST_SNO; i <= D_TEST_ENO; i++) { siTNo = i; fOne(); } #else fOne(); #endif return 0; }