Implementation of Simon-Park-Krug algorithm in C++ to simulate a Wright-Fisher model (beneficial mutation rate U and advantage s_0 are constant).