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ê] Entropy và thông tin tương hỗ từ phân phối đồng thời

    Cho phân phối đồng thời của hai biến ngẫu nhiên rời rạc XXX (nhận rrr giá trị) và YYY (nhận ccc giá trị), cho bởi bảng số đếm nguyên không âm nijn_{ij}nij​ (i=1..r, j=1..ci=1..r,\ j=1..ci=1..r, j=1..c), trong đó P(X=xi,Y=yj)=nij/NP(X=x_i, Y=y_j) = n_{ij}/NP(X=xi​,Y=yj​)=nij​/N với N=∑i,jnij>0N = \sum_{i,j} n_{ij} > 0N=∑i,j​nij​>0. Hãy tính entropy Shannon (đơn vị bit, logarit cơ số 2, quy ước 0log⁡20=00\log_2 0 = 00log2​0=0):

    H(X)=−∑iP(xi)log⁡2P(xi),H(Y)=−∑jP(yj)log⁡2P(yj)H(X) = -\sum_i P(x_i)\log_2 P(x_i), \qquad H(Y) = -\sum_j P(y_j)\log_2 P(y_j)H(X)=−∑i​P(xi​)log2​P(xi​),H(Y)=−∑j​P(yj​)log2​P(yj​) H(X,Y)=−∑i,jP(xi,yj)log⁡2P(xi,yj)H(X,Y) = -\sum_{i,j} P(x_i,y_j)\log_2 P(x_i,y_j)H(X,Y)=−∑i,j​P(xi​,yj​)log2​P(xi​,yj​) I(X;Y)=H(X)+H(Y)−H(X,Y)(thoˆng tin tương hoˆ˜ — mutual information)I(X;Y) = H(X) + H(Y) - H(X,Y) \quad \text{(thông tin tương hỗ — mutual information)}I(X;Y)=H(X)+H(Y)−H(X,Y)(thoˆng tin tương hoˆ˜ — mutual information)

    Ví dụ: với bảng 2×22\times 22×2 toàn số 111 (tức X,YX,YX,Y độc lập, mỗi ô xác suất 1/41/41/4): H(X)=H(Y)=1H(X)=H(Y)=1H(X)=H(Y)=1, H(X,Y)=2H(X,Y)=2H(X,Y)=2, I(X;Y)=0I(X;Y)=0I(X;Y)=0 (đúng như lý thuyết: hai biến độc lập có thông tin tương hỗ bằng 0).

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

      Dòng 1: hai số nguyên r cr\ cr c (1≤r,c≤201 \le r, c \le 201≤r,c≤20). rrr dòng tiếp theo, mỗi dòng gồm ccc số nguyên không âm ni1 … nicn_{i1}\ \dots\ n_{ic}ni1​ … nic​ (0≤nij≤1060 \le n_{ij} \le 10^60≤nij​≤106). Đảm bảo N=∑nij>0N=\sum n_{ij} > 0N=∑nij​>0.

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

      In 4 số thực trên 4 dòng riêng biệt theo đúng thứ tự H(X)H(X)H(X), H(Y)H(Y)H(Y), H(X,Y)H(X,Y)H(X,Y), I(X;Y)I(X;Y)I(X;Y), mỗi số làm tròn và hiển thị đúng 6 chữ số thập phân.

    Ví dụ:

    Đầu vào:

    2 2
    1 1
    1 1
    

    Đầu ra:

    1.000000
    1.000000
    2.000000
    0.000000
    

    Đầu vào:

    2 2
    1 0
    0 1
    

    Đầu ra:

    1.000000
    1.000000
    1.000000
    1.000000
    

    Đang tải editor...