The skewness of this distribution.
Computed from closed-form raw moments EX^r = a * p * S_r / (1-p)^r where S_1=1, S_2=1, S_3=(1+p), S_4=(1+4p+p²).