Programa de probabilidad

[Intael]

Hola,

Tengo que diseñar un programa que dé la probabilidad de sacar cara C veces al tirar una moneda T veces, asumiendo que 0=<t=<9 y 0=<c=<t.

Sé que esto se consigue a través de la distribución binomial, pero no acabo de entender cómo diseñar el algoritmo... ¿Alguna idea? Os dejo el enunciado aquí (en catalán):

Calculeu la probabilitat que, tirant una moneda t vegades, surti cara c vegades. Suposeu
que la moneda no està trucada.
Entrada
L’entrada consisteix en dos naturals t i c. Suposeu 0 t 9 i 0 c t.
Sortida
Escriviu, amb quatre decimals, la probabilitat d’obtenir c cares amb t tirades.
Observaci ´o
Si programeu en C++, feu servir el tipus double per fer els c`alculs, i poseu aquestes dues
l´ınies al principi del vostre main():
cout. setf ( ios :: fixed );
cout. precision (4);

[sam90]

Para calcular eso, hay una formula ya conocida:

Mira wikipedia: http://es.wikipedia.org/wiki/Distribuci%C3%B3n_binomial

[Intael]

Iniciado por sam90

Para calcular eso, hay una formula ya conocida:

Mira wikipedia: [url]http://es.wikipedia.org/wiki/Distribuci%C3%B3n_binomial[/url]

Pero implica el uso de factoriales... Existe alguna función para calcular factoriales? Por que si tengo que hacer una iteración para cada factorial, lo llevo claro.

[sam90]

Puede meter esa iteracion en una funcion.

Código C:

int factorial(int n)
{
    int a = 1,i;
    for (i = 1 ;i <= n ; i ++) a = a * i;
    return a;
}
 

[Intael]

Hola!

Usando la función que mencionas, me da 0 todo el rato... Qué he hecho mal? Aquí está el código:

#include <iostream>
#include <cmath>

using namespace std;
int n, k;
int factorial(int);

int main ()
{
int n, k, c;
double p, b;
cin>>n>>k;
p=0,5;
c=n-k;
b=((factorial(n))/(factorial(k)*factorial(c)))*pow(p,k)*pow(1-p,n-k);
cout<<b;

}

int factorial(int n)
{
int a = 1,i;
for (i = 1 ;i <= n ; i ++) a = a * i;
return a;
}

[CalgaryCorpus]

p = 0,5; asigna 0 a p. Usa el . para doubles.