#include #include #include #include #include #include #include // 内部定数 #define D_VAL_MAX 2000005 // 最大値 // 内部変数 static FILE *szpFpI; // 入力 static int siN; // N static double sdRate; // 残る確率 static double sd1Rate[D_VAL_MAX]; // 残る確率[値] // 内部変数 - テスト用 #ifdef D_TEST static int siRes; static FILE *szpFpA; #endif // 実行メイン int fMain( int piTNo // テスト番号 1~ ) { int i, j; 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 // N・残る確率 - 取得 fgets(lc1Buf, sizeof(lc1Buf), szpFpI); sscanf(lc1Buf, "%d%lf", &siN, &sdRate); sdRate = 1.0 - sdRate; // 残る確率[値] - 初期化 for (i = 2; i <= siN; i++) { sd1Rate[i] = 1.0; } // 残る確率[値] - 更新 for (i = 2; i <= siN; i++) { for (j = i * 2; j <= siN; j += i) { sd1Rate[j] *= sdRate; } } // 個数 - 取得 double ldCnt = 0.0; for (i = 2; i <= siN; i++) { ldCnt += sd1Rate[i]; } // 結果 - セット sprintf(lc1Out, "%.7lf\n", ldCnt); // 結果 - 表示 #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; }