/root/doris/be/src/util/random.h
| Line | Count | Source | 
| 1 |  | // Licensed to the Apache Software Foundation (ASF) under one | 
| 2 |  | // or more contributor license agreements.  See the NOTICE file | 
| 3 |  | // distributed with this work for additional information | 
| 4 |  | // regarding copyright ownership.  The ASF licenses this file | 
| 5 |  | // to you under the Apache License, Version 2.0 (the | 
| 6 |  | // "License"); you may not use this file except in compliance | 
| 7 |  | // with the License.  You may obtain a copy of the License at | 
| 8 |  | // | 
| 9 |  | //   http://www.apache.org/licenses/LICENSE-2.0 | 
| 10 |  | // | 
| 11 |  | // Unless required by applicable law or agreed to in writing, | 
| 12 |  | // software distributed under the License is distributed on an | 
| 13 |  | // "AS IS" BASIS, WITHOUT WARRANTIES OR CONDITIONS OF ANY | 
| 14 |  | // KIND, either express or implied.  See the License for the | 
| 15 |  | // specific language governing permissions and limitations | 
| 16 |  | // under the License. | 
| 17 |  |  | 
| 18 |  | #pragma once | 
| 19 |  |  | 
| 20 |  | #include <stdint.h> | 
| 21 |  |  | 
| 22 |  | namespace doris { | 
| 23 |  |  | 
| 24 |  | // A very simple random number generator.  Not especially good at | 
| 25 |  | // generating truly random bits, but good enough for our needs in this | 
| 26 |  | // package. | 
| 27 |  | class Random { | 
| 28 |  | private: | 
| 29 |  |     uint32_t seed_; | 
| 30 |  |  | 
| 31 |  | public: | 
| 32 | 905 |     explicit Random(uint32_t s) : seed_(s & 0x7fffffffu) { | 
| 33 |  |         // Avoid bad seeds. | 
| 34 | 905 |         if (seed_ == 0 || seed_ == 2147483647L) { | 
| 35 | 0 |             seed_ = 1; | 
| 36 | 0 |         } | 
| 37 | 905 |     } | 
| 38 | 22.4k |     uint32_t Next() { | 
| 39 | 22.4k |         static const uint32_t M = 2147483647L; // 2^31-1 | 
| 40 | 22.4k |         static const uint64_t A = 16807;       // bits 14, 8, 7, 5, 2, 1, 0 | 
| 41 |  |         // We are computing | 
| 42 |  |         //       seed_ = (seed_ * A) % M,    where M = 2^31-1 | 
| 43 |  |         // | 
| 44 |  |         // seed_ must not be zero or M, or else all subsequent computed values | 
| 45 |  |         // will be zero or M respectively.  For all other values, seed_ will end | 
| 46 |  |         // up cycling through every number in [1,M-1] | 
| 47 | 22.4k |         uint64_t product = seed_ * A; | 
| 48 |  |  | 
| 49 |  |         // Compute (product % M) using the fact that ((x << 31) % M) == x. | 
| 50 | 22.4k |         seed_ = static_cast<uint32_t>((product >> 31) + (product & M)); | 
| 51 |  |         // The first reduction may overflow by 1 bit, so we may need to | 
| 52 |  |         // repeat.  mod == M is not possible; using > allows the faster | 
| 53 |  |         // sign-bit-based test. | 
| 54 | 22.4k |         if (seed_ > M) { | 
| 55 | 0 |             seed_ -= M; | 
| 56 | 0 |         } | 
| 57 | 22.4k |         return seed_; | 
| 58 | 22.4k |     } | 
| 59 |  |     // Returns a uniformly distributed value in the range [0..n-1] | 
| 60 |  |     // REQUIRES: n > 0 | 
| 61 | 11.5k |     uint32_t Uniform(int n) { return Next() % n; } | 
| 62 |  |  | 
| 63 |  |     // Randomly returns true ~"1/n" of the time, and false otherwise. | 
| 64 |  |     // REQUIRES: n > 0 | 
| 65 | 0 |     bool OneIn(int n) { return (Next() % n) == 0; } | 
| 66 |  |  | 
| 67 |  |     // Skewed: pick "base" uniformly from range [0,max_log] and then | 
| 68 |  |     // return "base" random bits.  The effect is to pick a number in the | 
| 69 |  |     // range [0,2^max_log-1] with exponential bias towards smaller numbers. | 
| 70 | 0 |     uint32_t Skewed(int max_log) { return Uniform(1 << Uniform(max_log + 1)); } | 
| 71 |  | }; | 
| 72 |  |  | 
| 73 |  | } // namespace doris |