src/weighted_reservoir_sampler.h

// Copyright 2016 Google Inc. All rights reserved.
//
// Licensed under the Apache License, Version 2.0 (the "License");
// you may not use this file except in compliance with the License.
// You may obtain a copy of the License at
//
//     http://www.apache.org/licenses/LICENSE-2.0
//
// Unless required by applicable law or agreed to in writing, software
// distributed under the License is distributed on an "AS IS" BASIS,
// WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
// See the License for the specific language governing permissions and
// limitations under the License.

#ifndef SRC_WEIGHTED_RESERVOIR_SAMPLER_H_
#define SRC_WEIGHTED_RESERVOIR_SAMPLER_H_

#include <cassert>
#include <random>

namespace protobuf_mutator {

// Algorithm pick one item from the sequence of weighted items.
// https://en.wikipedia.org/wiki/Reservoir_sampling#Algorithm_A-Chao
//
// Example:
//   WeightedReservoirSampler<int> sampler;
//   for(int i = 0; i < size; ++i)
//     sampler.Pick(weight[i], i);
//   return sampler.GetSelected();
template <class T, class RandomEngine = std::default_random_engine>
class WeightedReservoirSampler {
 public:
  explicit WeightedReservoirSampler(RandomEngine* random) : random_(random) {}

  void Try(uint64_t weight, const T& item) {
    if (Pick(weight)) selected_ = item;
  }

  const T& selected() const { return selected_; }

  bool IsEmpty() const { return total_weight_ == 0; }

 private:
  bool Pick(uint64_t weight) {
    if (weight == 0) return false;
    total_weight_ += weight;
    return weight == total_weight_ || std::uniform_int_distribution<uint64_t>(
                                          1, total_weight_)(*random_) <= weight;
  }

  T selected_ = {};
  uint64_t total_weight_ = 0;
  RandomEngine* random_;
};

}  // namespace protobuf_mutator

#endif  // SRC_WEIGHTED_RESERVOIR_SAMPLER_H_