sum k=0 to infinity of 1/(3k+1)*binomial(4k,k)*p^k*(1-p)^(3k+1)