Exercise 18
A baker blends 600 raisins and 400 chocolate chips into a dough mix
and, from this, makes 500 cookies. (a) Find the probability that a
randomly picked cookie will have no raisins. (b) Find the probability
that a randomly picked cookie will have exactly two chocolate chips. (c)
Find the probability that a randomly chosen cookie will have at least
two bits (raisins or chips) in it.
#Lets use a function to calculate combinations using gmp
choose_gmp <- function(n, k) {
as.bigq(factorialZ(n) / (factorialZ(k) * factorialZ(n - k)))
}
# (a) Probability of no raisins
prob_no_raisins <- as.numeric(choose_gmp(400, 0) * choose_gmp(600, 500) / choose_gmp(1000, 500))
# (b) Probability of exactly 2 chocolate chips
prob_2_chips <- as.numeric(choose_gmp(400, 2) * choose_gmp(600, 498) / choose_gmp(1000, 500))
# (c) Probability of at least two bits
prob_at_least_two_bits <- 1 - prob_no_raisins - as.numeric(choose_gmp(400, 1) * choose_gmp(600, 499) / choose_gmp(1000, 500))
#Displaying the results
prob_no_raisins
## [1] 4.111948e-184
## [1] 7.946932e-178
## [1] 1
LS0tDQp0aXRsZTogIkRpc2N1c3Npb24gNSINCmF1dGhvcjogIkxhdXJhIEIiDQpkYXRlOiAiYHIgU3lzLkRhdGUoKWAiDQpvdXRwdXQ6IG9wZW5pbnRybzo6bGFiX3JlcG9ydA0KLS0tDQoNCmBgYHtyIGxvYWQtcGFja2FnZXMsIG1lc3NhZ2U9RkFMU0V9DQoNCmxpYnJhcnkoZ21wKQ0KDQpgYGANCg0KIyMjIEV4ZXJjaXNlIDE4DQpBIGJha2VyIGJsZW5kcyA2MDAgcmFpc2lucyBhbmQgNDAwIGNob2NvbGF0ZSBjaGlwcyBpbnRvIGEgZG91Z2ggbWl4IGFuZCwNCmZyb20gdGhpcywgbWFrZXMgNTAwIGNvb2tpZXMuDQooYSkgRmluZCB0aGUgcHJvYmFiaWxpdHkgdGhhdCBhIHJhbmRvbWx5IHBpY2tlZCBjb29raWUgd2lsbCBoYXZlIG5vIHJhaXNpbnMuDQooYikgRmluZCB0aGUgcHJvYmFiaWxpdHkgdGhhdCBhIHJhbmRvbWx5IHBpY2tlZCBjb29raWUgd2lsbCBoYXZlIGV4YWN0bHkgdHdvDQpjaG9jb2xhdGUgY2hpcHMuDQooYykgRmluZCB0aGUgcHJvYmFiaWxpdHkgdGhhdCBhIHJhbmRvbWx5IGNob3NlbiBjb29raWUgd2lsbCBoYXZlIGF0IGxlYXN0DQp0d28gYml0cyAocmFpc2lucyBvciBjaGlwcykgaW4gaXQuDQoNCg0KDQpgYGB7ciBjb2RlLWNodW5rLWxhYmVsfQ0KDQojTGV0cyB1c2UgYSBmdW5jdGlvbiB0byBjYWxjdWxhdGUgY29tYmluYXRpb25zIHVzaW5nIGdtcA0KY2hvb3NlX2dtcCA8LSBmdW5jdGlvbihuLCBrKSB7DQogIGFzLmJpZ3EoZmFjdG9yaWFsWihuKSAvIChmYWN0b3JpYWxaKGspICogZmFjdG9yaWFsWihuIC0gaykpKQ0KfQ0KDQojIChhKSBQcm9iYWJpbGl0eSBvZiBubyByYWlzaW5zDQpwcm9iX25vX3JhaXNpbnMgPC0gYXMubnVtZXJpYyhjaG9vc2VfZ21wKDQwMCwgMCkgKiBjaG9vc2VfZ21wKDYwMCwgNTAwKSAvIGNob29zZV9nbXAoMTAwMCwgNTAwKSkNCg0KIyAoYikgUHJvYmFiaWxpdHkgb2YgZXhhY3RseSAyIGNob2NvbGF0ZSBjaGlwcw0KcHJvYl8yX2NoaXBzIDwtIGFzLm51bWVyaWMoY2hvb3NlX2dtcCg0MDAsIDIpICogY2hvb3NlX2dtcCg2MDAsIDQ5OCkgLyBjaG9vc2VfZ21wKDEwMDAsIDUwMCkpDQoNCiMgKGMpIFByb2JhYmlsaXR5IG9mIGF0IGxlYXN0IHR3byBiaXRzDQpwcm9iX2F0X2xlYXN0X3R3b19iaXRzIDwtIDEgLSBwcm9iX25vX3JhaXNpbnMgLSBhcy5udW1lcmljKGNob29zZV9nbXAoNDAwLCAxKSAqIGNob29zZV9nbXAoNjAwLCA0OTkpIC8gY2hvb3NlX2dtcCgxMDAwLCA1MDApKQ0KDQojRGlzcGxheWluZyB0aGUgcmVzdWx0cw0KcHJvYl9ub19yYWlzaW5zDQpwcm9iXzJfY2hpcHMNCnByb2JfYXRfbGVhc3RfdHdvX2JpdHMNCg0KYGBgDQoNCg0KDQo=