#include #include #include #include #include #include #include // 内部定数 #define D_ON 1 // 汎用フラグ - ON #define D_OFF 0 // 汎用フラグ - OFF #define D_PRM_MAX 100000 // 最大素数リスト数 // 内部構造体 - 素数情報 typedef struct Prm { int miVal; // 素数 int miCnt; // 数 } Prm; // 内部変数 static FILE *szpFpI; // 入力 static int siN; // N static int siK; // K static Prm sz1Prm[D_PRM_MAX]; // 素数リスト static char sc1PWork[D_PRM_MAX + 5]; // 素数リスト作成用ワーク static int siPCnt; // 素数リスト数 static int siAns; // 解答 - 値 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 fMakePrm( int piMax // 最大値 ) { // 初期化 siPCnt = 0; memset(sc1PWork, D_OFF, sizeof(sc1PWork)); // 開始値 int liNow = 2; // 上限値 int liLimit = (int)sqrt((double)piMax); // 作成 - 開始 while (1) { // 素数の検索 while (sc1PWork[liNow] != D_OFF) { liNow++; } // 上限チェック if (liNow > liLimit) { break; } // 素数リスト - 追加 sz1Prm[siPCnt].miVal = liNow; sz1Prm[siPCnt].miCnt = 0; siPCnt++; // ワークに素数以外をセット int liVal = liNow; while (liVal <= piMax) { sc1PWork[liVal] = D_ON; liVal += liNow; } } return 0; } // 素因数分解 int fDivPrm( int piVal // 値 ) { int liNo = 0; while (piVal > 1) { while (piVal % sz1Prm[liNo].miVal == 0) { piVal /= sz1Prm[liNo].miVal; sz1Prm[liNo].miCnt++; } if (liNo < siPCnt - 1) { liNo++; } else { break; } } return 0; } // 解答 - セット int fSetAns( int piPNo // 素数 0~ , int piVal // 値 , int piCCnt // 共通素数数 , int piDCnt // 約数数 ) { int i; // 選択中 if (piPNo < siPCnt) { // 個数でループ for (i = 0; ; i++) { // 値 if (piVal >= siN) { break; } // 共通素数数 int liCCnt; if (i < sz1Prm[piPNo].miCnt) { liCCnt = i; } else { liCCnt = sz1Prm[piPNo].miCnt; } // 下位へ fSetAns(piPNo + 1, piVal, piCCnt + liCCnt, piDCnt * (i + 1)); // 次の値 piVal *= sz1Prm[piPNo].miVal; } return 0; } // 共通素数数 if (piCCnt < siK) { return 0; } // 約数数 if (siDCnt > piDCnt) { return 0; } else if (siDCnt == piDCnt) { if (siAns < piVal) { return 0; } } // 解答 - セット siAns = piVal; siDCnt = piDCnt; return 0; } // 実行メイン int fMain( ) { char lc1Buf[1024]; // N, K - 取得 fgets(lc1Buf, sizeof(lc1Buf), szpFpI); sscanf(lc1Buf, "%d%d", &siN, &siK); // 素数リスト - 作成 fMakePrm(siN); // 素因数分解 fDivPrm(siN); // 解答 - セット fSetAns(0, 1, 0, 1); return 0; } // 1回実行 int fOne( ) { int liRet; char lc1Buf[1024]; // データ - 初期化 siAns = 0; // 解答 - 値 siDCnt = 0; // 解答 - 約数数 // 入力 - セット #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", siAns); 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; }