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ê] Xấp xỉ nhị thức bằng CLT với hiệu chỉnh liên tục — truy vấn hàng loạt

    Thực hiện nnn phép thử Bernoulli độc lập, mỗi phép thử có xác suất thành công p=a/bp = a/bp=a/b (với a,ba, ba,b là số nguyên cho trước, 0≤a≤b0 \le a \le b0≤a≤b). Gọi XXX là số lần thành công trong nnn phép thử. Theo CLT, XXX xấp xỉ phân phối chuẩn N(np, np(1−p))N\big(np,\ np(1-p)\big)N(np, np(1−p)).

    Với mỗi truy vấn gồm cặp (k1,k2)(k_1, k_2)(k1​,k2​) (0≤k1≤k2≤n0 \le k_1 \le k_2 \le n0≤k1​≤k2​≤n, nguyên), hãy tính xấp xỉ P(k1≤X≤k2)P(k_1 \le X \le k_2)P(k1​≤X≤k2​) dùng hiệu chỉnh liên tục:

    P(k1≤X≤k2)≈Φ ⁣(k2+0.5−npnp(1−p))−Φ ⁣(k1−0.5−npnp(1−p))P(k_1\le X\le k_2) \approx \Phi\!\left(\frac{k_2+0.5-np}{\sqrt{np(1-p)}}\right) - \Phi\!\left(\frac{k_1-0.5-np}{\sqrt{np(1-p)}}\right)P(k1​≤X≤k2​)≈Φ(np(1−p)​k2​+0.5−np​)−Φ(np(1−p)​k1​−0.5−np​)

    Quy ước trường hợp suy biến: nếu p=0p=0p=0 hoặc p=1p=1p=1 (khi đó np(1−p)=0np(1-p)=0np(1−p)=0), XXX tất định bằng npnpnp: kết quả là 111 nếu k1≤np≤k2k_1 \le np \le k_2k1​≤np≤k2​, ngược lại là 000.

    Ví dụ: n=10,p=1/2n=10, p=1/2n=10,p=1/2. Với truy vấn (3,7)(3,7)(3,7): np=5np=5np=5, np(1−p)=2.5≈1.5811\sqrt{np(1-p)}=\sqrt{2.5}\approx1.5811np(1−p)​=2.5​≈1.5811; kết quả xấp xỉ ≈0.886154\approx 0.886154≈0.886154.

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

      Dòng đầu tiên chứa ba số nguyên n a bn\ a\ bn a b (1≤n≤1071 \le n \le 10^71≤n≤107, 0≤a≤b0 \le a \le b0≤a≤b, 1≤b≤1061 \le b \le 10^61≤b≤106) — số phép thử và phân số biểu diễn p=a/bp=a/bp=a/b. Dòng tiếp theo chứa số nguyên QQQ (1≤Q≤10001 \le Q \le 10001≤Q≤1000) — số truy vấn. QQQ dòng tiếp theo, mỗi dòng chứa hai số nguyên k1 k2k_1\ k_2k1​ k2​ (0≤k1≤k2≤n0 \le k_1 \le k_2 \le n0≤k1​≤k2​≤n).

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

      Với mỗi truy vấn, in ra một dòng là xác suất xấp xỉ, làm tròn đến đúng 6 chữ số thập phân.

    Ví dụ:

    Đầu vào:

    100 3 10
    3
    25 35
    0 100
    50 50
    

    Đầu ra:

    0.769938
    1.000000
    0.000007
    

    Đầu vào:

    10 1 2
    2
    3 7
    0 10
    

    Đầu ra:

    0.886154
    0.999496
    

    Đang tải editor...