Submission #1797152

---

Source Code Expand

#include <cstdio>
#include <algorithm>
using namespace std;
#define MAX_N 100
#define MAX_A 70
#define MAX_B 40
#define MAX_C 30
#define p double(1) / double(6)
// dp[n][a][b][c] : サイコロをn回振った時に素因数2,3,5がそれぞれa,b,c個以上含まれる確率
double dp[MAX_N][MAX_A][MAX_B][MAX_C];

int main(int argc, const char * argv[]) {
    int N;
    long long D;
    scanf("%d%lld", &N, &D);
    
    const long long prime[] = {2, 3, 5};
    int r[3] = {0, 0, 0};  // r[i] : prime[i]が含まれる個数
    for (int i = 0; i < 3; i++) {
        while (D % prime[i] == 0) {
            D = D / prime[i];
            r[i] += 1;
        }
    }
    // dp[n][a][b][c] : サイコロをn回振った時に素因数2,3,5がそれぞれa,b,c個以上含まれる確率
    //float dp[MAX_N][MAX_A][MAX_B][MAX_C];
    for (int n = 0; n < N+1; n++) {
        for (int a = 0; a < r[0]+1; a++) {
            for (int b = 0; b < r[1]+1; b++) {
                for (int c = 0; c < r[2]+1; c++) {
                    if (n == 0) {
                        if (a == 0 && b == 0 && c == 0) {
                            dp[n][a][b][c] = double(1);
                        } else {
                            dp[n][a][b][c] = double(0);
                        }
                    } else {
                        dp[n][a][b][c] =
                              p * dp[n-1][a][b][c]
                            + p * dp[n-1][max(a-1, 0)][b][c]
                            + p * dp[n-1][a][max(b-1, 0)][c]
                            + p * dp[n-1][max(a-2, 0)][b][c]
                            + p * dp[n-1][a][b][max(c-1, 0)]
                            + p * dp[n-1][max(a-1, 0)][max(b-1, 0)][c];
                    }
                }
            }
        }
    }
    
    printf("%f\n", dp[N][r[0]][r[1]][r[2]]);
    
    return 0;
}

Submission Info

Compile Error

./Main.cpp: In function ‘int main(int, const char**)’:
./Main.cpp:15:28: warning: ignoring return value of ‘int scanf(const char*, ...)’, declared with attribute warn_unused_result [-Wunused-result]
     scanf("%d%lld", &N, &D);
                            ^

Judge Result