#include #include #include #include #include #include #include // 内部定数 #define D_ON 1 // 汎用フラグ - ON #define D_OFF 0 // 汎用フラグ - OFF #define D_VAL_MAX 20005 // 最大値 #define D_PRM_MAX D_VAL_MAX // 最大素数リスト数 // 内部変数 static FILE *szpFpI; // 入力 static int siSum; // 和 static int si2Cnt[D_PRM_MAX / 5][D_VAL_MAX]; // 回数 static int si1Prm[D_PRM_MAX]; // 素数リスト static char sc1PWork[D_PRM_MAX + 5]; // 素数リスト作成用ワーク static int siPCnt; // 素数リスト数 // 内部変数 - テスト用 #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; } // 素数リスト - 追加 si1Prm[siPCnt] = liNow; siPCnt++; // ワークに素数以外をセット int liVal = liNow; while (liVal <= piMax) { sc1PWork[liVal] = D_ON; liVal += liNow; } } // 残った素数を追加 while (liNow <= piMax) { if (sc1PWork[liNow] == D_OFF) { si1Prm[siPCnt] = liNow; siPCnt++; } liNow++; } return 0; } // 素数 - 選択 int fSelPrm( int piPNo // 素数 0~ , int piSum // 合計 , int piCnt // 回数 ) { // 選択終了 if (piPNo >= siPCnt) { return 0; } // 合計 - チェック if (piSum > siSum) { return 0; } // 回数 - チェック・更新 if (si2Cnt[piPNo][piSum] >= piCnt) { return 0; } si2Cnt[piPNo][piSum] = piCnt; // 選択 fSelPrm(piPNo + 1, piSum + si1Prm[piPNo], piCnt + 1); // 選択なし fSelPrm(piPNo + 1, piSum, piCnt); return 0; } // 実行メイン int fMain( ) { char lc1Buf[1024]; // 和 - 取得 fgets(lc1Buf, sizeof(lc1Buf), szpFpI); sscanf(lc1Buf, "%d", &siSum); // 素数リスト - 作成 fMakePrm(siSum); // 素数 - 選択 fSelPrm(0, 0, 1); return 0; } // 1回実行 int fOne( ) { int liRet; char lc1Buf[1024]; // データ - 初期化 memset(si2Cnt, 0, sizeof(si2Cnt)); // 回数 // 入力 - セット #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", si2Cnt[siPCnt - 1][siSum] - 1); // 結果 - 出力 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; }