結果

問題 No.1688 Veterinarian
ユーザー ID 21712
提出日時 2025-07-19 00:20:32
言語 Go
(1.23.4)
結果
AC  
実行時間 125 ms / 3,000 ms
コード長 1,361 bytes
コンパイル時間 13,242 ms
コンパイル使用メモリ 236,660 KB
実行使用メモリ 8,096 KB
最終ジャッジ日時 2025-07-19 00:20:47
合計ジャッジ時間 13,197 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 14
権限があれば一括ダウンロードができます

ソースコード

diff # 

package main

import . "fmt"

func main() {
	var a,b,c,n int
	Scan(&a,&b,&c,&n)
	dp := make([][][]float64, a+1)
	for i := range dp {
		dp[i] = make([][]float64, b+1)
		for j := range dp[i] {
			dp[i][j] = make([]float64, c+1)
		}
	}
	dp[0][0][0] = 1
	for i := 0; i < n; i++ {
		tmp := dp
		dp = make([][][]float64, a+1)
		for i := range dp {
			dp[i] = make([][]float64, b+1)
			for j := range dp[i] {
				dp[i][j] = make([]float64, c+1)
			}
		}
		for aa := 0; aa <= a; aa++ {
			for bb := 0; bb <= b; bb++ {
				for cc := 0; cc <= c; cc++ {
					cnt := a+b+c-(aa+bb+cc)
					var ap,bp,cp float64
					if a-aa>=2 {
						ap = float64((a-aa)*(a-aa-1))/float64(cnt*(cnt-1))
						dp[aa+1][bb][cc] += ap*tmp[aa][bb][cc]
					}
					if b-bb>=2 {
						bp = float64((b-bb)*(b-bb-1))/float64(cnt*(cnt-1))
						dp[aa][bb+1][cc] += bp*tmp[aa][bb][cc]
					}
					if c-cc>=2 {
						cp = float64((c-cc)*(c-cc-1))/float64(cnt*(cnt-1))
						dp[aa][bb][cc+1] += cp*tmp[aa][bb][cc]
					}
					dp[aa][bb][cc] += (1-(ap+bp+cp))*tmp[aa][bb][cc]
				}
			}
		}
	}
	var ans [3]float64
	for aa := 0; aa <= a; aa++ {
		for bb := 0; bb <= b; bb++ {
			for cc := 0; cc <= c; cc++ {
				ans[0] += float64(aa) * dp[aa][bb][cc]
				ans[1] += float64(bb) * dp[aa][bb][cc]
				ans[2] += float64(cc) * dp[aa][bb][cc]
			}
		}
	}
	Printf("%.10f %.10f %.10f\n", ans[0], ans[1], ans[2])
}
0