Các lượt nộp
    Danh sách bài
    Trang chủ
    Báo lỗi

    solution

    Đề bài: [Xác suất - Thống kê] Suy diễn hậu nghiệm trên mạng Bayes chuỗi

    Cho một mạng Bayes dạng chuỗi A→B→CA \to B \to CA→B→C, trong đó AAA có nAn_AnA​ trạng thái, BBB có nBn_BnB​ trạng thái, CCC có nCn_CnC​ trạng thái (đánh số từ 1). Cho:

    • Phân phối tiên nghiệm P(A)P(A)P(A): vector nAn_AnA​ số.
    • Bảng xác suất có điều kiện P(B∣A)P(B \mid A)P(B∣A): ma trận nA×nBn_A \times n_BnA​×nB​, hàng iii là phân phối của BBB khi A=iA=iA=i.
    • Bảng xác suất có điều kiện P(C∣B)P(C \mid B)P(C∣B): ma trận nB×nCn_B \times n_CnB​×nC​, hàng jjj là phân phối của CCC khi B=jB=jB=j.

    Biến CCC được quan sát ở trạng thái ccc. Vì CCC chỉ phụ thuộc BBB và BBB chỉ phụ thuộc AAA (tính chất Markov của mạng Bayes), hậu nghiệm của AAA được tính bằng cách liệt kê (enumerate) qua toàn bộ giá trị của BBB:

    P(A=a∣C=c)∝P(A=a)∑b=1nBP(B=b∣A=a) P(C=c∣B=b)P(A=a \mid C=c) \propto P(A=a) \sum_{b=1}^{n_B} P(B=b \mid A=a)\, P(C=c \mid B=b)P(A=a∣C=c)∝P(A=a)∑b=1nB​​P(B=b∣A=a)P(C=c∣B=b)

    Hãy tính phân phối hậu nghiệm đầy đủ P(A∣C=c)P(A \mid C=c)P(A∣C=c) cho mọi a=1,…,nAa = 1, \dots, n_Aa=1,…,nA​.

    • Định dạng đầu vào:

      Dòng 1: bốn số nguyên nA nB nC cn_A\ n_B\ n_C\ cnA​ nB​ nC​ c (1≤c≤nC1 \le c \le n_C1≤c≤nC​). Dòng 2: nAn_AnA​ số thực — P(A)P(A)P(A). nAn_AnA​ dòng tiếp theo: mỗi dòng nBn_BnB​ số thực — hàng của bảng P(B∣A)P(B\mid A)P(B∣A). nBn_BnB​ dòng tiếp theo: mỗi dòng nCn_CnC​ số thực — hàng của bảng P(C∣B)P(C\mid B)P(C∣B).

    • Định dạng đầu ra:

      In ra nAn_AnA​ số thực trên một dòng, cách nhau bởi khoảng trắng, là P(A=1∣C=c),…,P(A=nA∣C=c)P(A=1\mid C=c), \dots, P(A=n_A\mid C=c)P(A=1∣C=c),…,P(A=nA​∣C=c), mỗi số làm tròn 6 chữ số thập phân. Nếu P(C=c)=0P(C=c)=0P(C=c)=0 (tổng chuẩn hoá bằng 0, sự kiện không thể xảy ra), in ra UNDEFINED.

    Ví dụ:

    Đầu vào:

    2 2 2 1
    0.5 0.5
    0.9 0.1
    0.2 0.8
    0.8 0.2
    0.3 0.7
    

    Đầu ra:

    0.652174 0.347826
    

    Đầu vào:

    1 2 2 2
    1.0
    0.4 0.6
    0.9 0.1
    0.2 0.8
    

    Đầu ra:

    1.000000
    

    Đang tải editor...