t=Table[RandomVariate[HypergeometricDistribution[129,167,299]],
1000000];//AbsoluteTiming
p1=Count[t,u_/;u>=86]/1000000.
{136.942,Null}
0.001054

sf=Join[s=Table[1,167],f=Table[0,132]];
t=Table[RandomSample[sf,129]//Total,10000000];//AbsoluteTiming
p2=Count[t,u_/;u>=86]/10000000.
{18.5029,Null}
0.000744

p3=1-CDF[HypergeometricDistribution[129,167,299],86-1]//N
0.00074213

p1/p3
p2/p3
1.42024
1.00252