Use the hypergeometric distribution

This will give us a factor about 2.5

The r-code is improved :
- use exact values for afc until 1754 (max that can be
  computed with "long double")
- use per-thread static variables to allow parallel run
  in different threads