#14797 (Go1.4) No.162 8020運動

• 問題: No.162 8020運動
• ユーザー: yuki2006
• 提出日時: 2015-03-05 22:15:16
• 言語: Go1.4 (1.4.2)
• 結果: AC
• 実行時間: 165 ms / 5,000 ms
• コード長: 2,277 bytes
• コンパイル時間: 417 ms
• コンパイル使用メモリ: 35,304 KB
• 実行使用メモリ: 6,820 KB
• 最終ジャッジ日時: 2024-11-25 00:10:44
• 合計ジャッジ時間: 3,425 ms
• ジャッジサーバーID （参考情報）: judge4 / judge2

ファイルパターン	結果
sample	AC * 3
other	AC * 26

ソースコード

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
package main
import (
    "bufio"
    "os"
    "strconv"
    "fmt"
    "strings"
)
var s = bufio.NewScanner(os.Stdin)
func next() string {
    s.Split(bufio.ScanWords)
    s.Scan()
    return s.Text()
}
func nextInt() int {
    i, e := strconv.Atoi(next())
    if e != nil {
        panic(e)
    }
    return i
}
func nextLong() int64 {
    i, e := strconv.ParseInt(next(), 10, 64)
    if e != nil {
        panic(e)
    }
    return i
}
func nextIntArray() []int {
    s.Split(bufio.ScanLines)
    s.Scan()
    S := s.Text()
    vs := strings.Split(S, " ")
    N := len(vs)
    data := make([]int, N)
    for i := 0; i < N; i++ {
        v, e := strconv.Atoi(vs[i])
        if e != nil {
            panic(e)
        }
        data[i] = v
    }
    return data
}
/*
 あと N年で歯が連続でK本残っている場合のK本残る期待値
 */
func recv(N uint, K uint, P []int) float64 {
    if N == 0 || K == 0 {
        return float64(K)
    }
    if memo[N][K] != -1 {
        return memo[N][K]
    }
    if K == 1 {
        memo[N][K] = ((100.0-float64(P[0]))/100.0)*recv(N-1, K, P)
    }else {
        memo[N][K] = 0
        for mask := uint(0); mask < (1<<K); mask++ {
            prob := 1.0
            for i := uint(0); i < K; i++ {
                //左端と右端の確率
                bad := 0.0
                if i == 0 || i == K-1 {
                    //連続という前提なので 1 や K-2のは歯もある
                    bad = float64(P[1])/100.0
                }else {
                    bad = float64(P[2])/100.0
                }
                //歯が残る
                if (mask&(1<<i)) > 0 {
                    prob*=1-bad
                }else {
                    prob*=bad
                }
            }
            left := uint(0)
            if prob==0.0{
                continue
            }
            for i := uint(0); i < K; i++ {
                if (mask&(1<<i)) == 0 {
                    //歯の連続が途切れた時の確率
                    memo[N][K]+=prob*recv(N-1, left, P)
                    left = 0
                }else {
                    //歯が残っている
                    left++
                }
            }
            memo[N][K]+=prob*recv(N-1, left, P)
        }
    }
    return memo[N][K]
}
func fill(arr []float64, v float64) {
    for i := 0; i < len(arr); i++ {
        arr[i] = v
    }
}
var memo [][]float64
func main() {
    A := nextInt()
    //P[0] 両方の歯がないときの虫歯になる確率
    //P[1] 片方の歯がないときの虫歯になる確率
    //P[0] 両方の歯があるときの虫歯になる確率
    memo = make([][]float64, 80*20)
    P := nextIntArray()
    for i := 0; i < 80; i++ {
        memo[i] = make([]float64, 20)
        fill(memo[i], -1)
    }
    //  fmt.Println(memo)
    fmt.Println(2 * recv(uint(80 - A), uint(14), P))
}
 
 
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX