#include #include #include #include #include #include #include // 内部定数 #define D_ON 1 // 汎用フラグ - ON #define D_OFF 0 // 汎用フラグ - OFF #define D_PRM_MAX 10000000 // 最大素数リスト数 // 内部変数 static FILE *szpFpI; // 入力 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; } // 実行メイン long long fMain( ) { int i; char lc1Buf[1024]; // 鴨数・距離 - 取得 int liCnt, liLen; fgets(lc1Buf, sizeof(lc1Buf), szpFpI); sscanf(lc1Buf, "%d%d", &liCnt, &liLen); // 素数リスト - 作成 fMakePrm(liLen / 2); // パターン数 - 取得 long long llPCnt = 0; for (i = 0; i < siPCnt; i++) { // 最低距離 int liMin = (liCnt - 1) * si1Prm[i]; if (liMin > liLen) { break; } // パターン数 - 加算 llPCnt += liLen - liMin + 1; } return llPCnt; } // 1回実行 int fOne( ) { long long llRet; 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 // 実行メイン llRet = fMain(); // 結果 - セット sprintf(lc1Buf, "%lld\n", llRet); // 結果 - 出力 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; }