2 * Copyright (c) 2017 Mohamed Aslan <maslan@sce.carleton.ca>
4 * Permission to use, copy, modify, and distribute this software for any
5 * purpose with or without fee is hereby granted, provided that the above
6 * copyright notice and this permission notice appear in all copies.
8 * THE SOFTWARE IS PROVIDED "AS IS" AND THE AUTHOR DISCLAIMS ALL WARRANTIES
9 * WITH REGARD TO THIS SOFTWARE INCLUDING ALL IMPLIED WARRANTIES OF
10 * MERCHANTABILITY AND FITNESS. IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR
11 * ANY SPECIAL, DIRECT, INDIRECT, OR CONSEQUENTIAL DAMAGES OR ANY DAMAGES
12 * WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR PROFITS, WHETHER IN AN
13 * ACTION OF CONTRACT, NEGLIGENCE OR OTHER TORTIOUS ACTION, ARISING OUT OF
14 * OR IN CONNECTION WITH THE USE OR PERFORMANCE OF THIS SOFTWARE.
25 nCr_r(uint32_t n
, uint32_t r
)
31 return nCr_r(n
- 1, r
) + nCr_r(n
- 1, r
- 1);
35 nCr(uint32_t n
, uint32_t r
)
40 for (i
= 0 ; i
< 101 ; i
++)
41 C
[i
] = (uint64_t *)xcalloc(101, sizeof(uint64_t));
42 for (i
= 1 ; i
<= n
; i
++) {
43 for (j
= 1 ; j
<= r
; j
++) {
49 C
[i
][j
] = C
[i
- 1][j
] + C
[i
- 1][j
- 1];
52 for (i
= 0 ; i
< 101 ; i
++)
58 rw2phi(int r
, int w
, int n
)
61 return 1.0 - ((double)nCr(n
- w
, r
) / (double)nCr(n
, r
));
67 phi2rw(double phi
, int n
, int *r
, int *w
)
69 int i
, i_min
= 1, j
, j_min
= 1;
70 double d
, d_min
= 1.1;
72 for (i
= 1 ; i
< n
; i
++) {
73 for (j
= i
; j
<= n
- i
; j
++) {
74 d
= fabs(rw2phi(i
, j
, n
) - phi
);