Create a gist now

Instantly share code, notes, and snippets.

Embed
Bitcoin AttackerSuccessProbability
#include <math.h>
double AttackerSuccessProbability(double q, int z){
double p = 1.0 - q;
double lambda = z * (q / p);
double sum = 1.0;
int i, k;
for (k = 0; k <= z; k++){
double poisson = exp(-lambda);
for (i = 1; i <= k; i++)poisson *= lambda / i;
sum -= poisson * (1 - pow(q / p, z - k));
}
return sum;
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment