Thực hiện n phép thử Bernoulli độc lập, mỗi phép thử có xác suất thành công p=a/b (với a,b là số nguyên cho trước, 0≤a≤b). Gọi X là số lần thành công trong n phép thử. Theo CLT, X xấp xỉ phân phối chuẩn N(np, np(1−p)).
Với mỗi truy vấn gồm cặp (k1,k2) (0≤k1≤k2≤n, nguyên), hãy tính xấp xỉ P(k1≤X≤k2) dùng hiệu chỉnh liên tục:
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=0 hoặc p=1 (khi đó np(1−p)=0), X tất định bằng np: kết quả là 1 nếu k1≤np≤k2, ngược lại là 0.
Ví dụ: n=10,p=1/2. Với truy vấn (3,7): np=5, np(1−p)=2.5≈1.5811; kết quả xấp xỉ ≈0.886154.
Dòng đầu tiên chứa ba số nguyên n a b (1≤n≤107, 0≤a≤b, 1≤b≤106) — số phép thử và phân số biểu diễn p=a/b. Dòng tiếp theo chứa số nguyên Q (1≤Q≤1000) — số truy vấn. Q dòng tiếp theo, mỗi dòng chứa hai số nguyên k1 k2 (0≤k1≤k2≤n).
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...