A 
geometric_distribution random number distribution
produces integer values 
i≥0
distributed according to
the discrete probability function
P(i|p)=p⋅(1−p)i .
template<class IntType = int>
  class geometric_distribution {
  public:
        using result_type = IntType;
    using param_type  = unspecified;
        geometric_distribution() : geometric_distribution(0.5) {}
    explicit geometric_distribution(double p);
    explicit geometric_distribution(const param_type& parm);
    void reset();
        template<class URBG>
      result_type operator()(URBG& g);
    template<class URBG>
      result_type operator()(URBG& g, const param_type& parm);
        double p() const;
    param_type param() const;
    void param(const param_type& parm);
    result_type min() const;
    result_type max() const;
  };