hash_table8.hpp 56 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499500501502503504505506507508509510511512513514515516517518519520521522523524525526527528529530531532533534535536537538539540541542543544545546547548549550551552553554555556557558559560561562563564565566567568569570571572573574575576577578579580581582583584585586587588589590591592593594595596597598599600601602603604605606607608609610611612613614615616617618619620621622623624625626627628629630631632633634635636637638639640641642643644645646647648649650651652653654655656657658659660661662663664665666667668669670671672673674675676677678679680681682683684685686687688689690691692693694695696697698699700701702703704705706707708709710711712713714715716717718719720721722723724725726727728729730731732733734735736737738739740741742743744745746747748749750751752753754755756757758759760761762763764765766767768769770771772773774775776777778779780781782783784785786787788789790791792793794795796797798799800801802803804805806807808809810811812813814815816817818819820821822823824825826827828829830831832833834835836837838839840841842843844845846847848849850851852853854855856857858859860861862863864865866867868869870871872873874875876877878879880881882883884885886887888889890891892893894895896897898899900901902903904905906907908909910911912913914915916917918919920921922923924925926927928929930931932933934935936937938939940941942943944945946947948949950951952953954955956957958959960961962963964965966967968969970971972973974975976977978979980981982983984985986987988989990991992993994995996997998999100010011002100310041005100610071008100910101011101210131014101510161017101810191020102110221023102410251026102710281029103010311032103310341035103610371038103910401041104210431044104510461047104810491050105110521053105410551056105710581059106010611062106310641065106610671068106910701071107210731074107510761077107810791080108110821083108410851086108710881089109010911092109310941095109610971098109911001101110211031104110511061107110811091110111111121113111411151116111711181119112011211122112311241125112611271128112911301131113211331134113511361137113811391140114111421143114411451146114711481149115011511152115311541155115611571158115911601161116211631164116511661167116811691170117111721173117411751176117711781179118011811182118311841185118611871188118911901191119211931194119511961197119811991200120112021203120412051206120712081209121012111212121312141215121612171218121912201221122212231224122512261227122812291230123112321233123412351236123712381239124012411242124312441245124612471248124912501251125212531254125512561257125812591260126112621263126412651266126712681269127012711272127312741275127612771278127912801281128212831284128512861287128812891290129112921293129412951296129712981299130013011302130313041305130613071308130913101311131213131314131513161317131813191320132113221323132413251326132713281329133013311332133313341335133613371338133913401341134213431344134513461347134813491350135113521353135413551356135713581359136013611362136313641365136613671368136913701371137213731374137513761377137813791380138113821383138413851386138713881389139013911392139313941395139613971398139914001401140214031404140514061407140814091410141114121413141414151416141714181419142014211422142314241425142614271428142914301431143214331434143514361437143814391440144114421443144414451446144714481449145014511452145314541455145614571458145914601461146214631464146514661467146814691470147114721473147414751476147714781479148014811482148314841485148614871488148914901491149214931494149514961497149814991500150115021503150415051506150715081509151015111512151315141515151615171518151915201521152215231524152515261527152815291530153115321533153415351536153715381539154015411542154315441545154615471548154915501551155215531554155515561557155815591560156115621563156415651566156715681569157015711572157315741575157615771578157915801581158215831584158515861587158815891590159115921593159415951596159715981599160016011602160316041605160616071608160916101611161216131614161516161617161816191620162116221623162416251626162716281629163016311632163316341635163616371638163916401641164216431644164516461647164816491650165116521653165416551656165716581659166016611662166316641665166616671668166916701671167216731674167516761677167816791680168116821683168416851686168716881689169016911692169316941695169616971698169917001701170217031704170517061707170817091710171117121713171417151716171717181719172017211722172317241725172617271728172917301731173217331734173517361737173817391740174117421743174417451746174717481749175017511752175317541755175617571758175917601761176217631764176517661767176817691770177117721773177417751776177717781779178017811782178317841785178617871788
  1. // emhash8::HashMap for C++11/14/17
  2. // version 1.6.3
  3. //
  4. // Licensed under the MIT License <http://opensource.org/licenses/MIT>.
  5. // SPDX-License-Identifier: MIT
  6. // Copyright (c) 2019-2022 Huang Yuanbing & bailuzhou AT 163.com
  7. //
  8. // Permission is hereby granted, free of charge, to any person obtaining a copy
  9. // of this software and associated documentation files (the "Software"), to deal
  10. // in the Software without restriction, including without limitation the rights
  11. // to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
  12. // copies of the Software, and to permit persons to whom the Software is
  13. // furnished to do so, subject to the following conditions:
  14. //
  15. // The above copyright notice and this permission notice shall be included in all
  16. // copies or substantial portions of the Software.
  17. //
  18. // THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
  19. // IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
  20. // FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
  21. // AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
  22. // LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
  23. // OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
  24. // SOFTWARE
  25. #pragma once
  26. #include <cstring>
  27. #include <string>
  28. #include <cstdlib>
  29. #include <type_traits>
  30. #include <cassert>
  31. #include <utility>
  32. #include <cstdint>
  33. #include <functional>
  34. #include <iterator>
  35. #include <algorithm>
  36. #ifdef EMH_KEY
  37. #undef EMH_KEY
  38. #undef EMH_VAL
  39. #undef EMH_KV
  40. #undef EMH_BUCKET
  41. #undef EMH_NEW
  42. #undef EMH_EMPTY
  43. #undef EMH_PREVET
  44. #undef EMH_LIKELY
  45. #undef EMH_UNLIKELY
  46. #endif
  47. // likely/unlikely
  48. #if defined(__GNUC__) || defined(__INTEL_COMPILER) || defined(__clang__)
  49. # define EMH_LIKELY(condition) __builtin_expect(condition, 1)
  50. # define EMH_UNLIKELY(condition) __builtin_expect(condition, 0)
  51. #else
  52. # define EMH_LIKELY(condition) (condition)
  53. # define EMH_UNLIKELY(condition) (condition)
  54. #endif
  55. #define EMH_KEY(p, n) p[n].first
  56. #define EMH_VAL(p, n) p[n].second
  57. #define EMH_KV(p, n) p[n]
  58. #define EMH_INDEX(i, n) i[n]
  59. #define EMH_BUCKET(i, n) i[n].bucket
  60. #define EMH_HSLOT(i, n) i[n].slot
  61. #define EMH_SLOT(i, n) (i[n].slot & _mask)
  62. #define EMH_PREVET(i, n) i[n].slot
  63. #define EMH_KEYMASK(key, mask) ((size_type)(key) & ~mask)
  64. #define EMH_EQHASH(n, key_hash) (EMH_KEYMASK(key_hash, _mask) == (_index[n].slot & ~_mask))
  65. #define EMH_NEW(key, val, bucket, key_hash) \
  66. new(_pairs + _num_filled) value_type(key, val); \
  67. _etail = bucket; \
  68. _index[bucket] = {bucket, _num_filled++ | EMH_KEYMASK(key_hash, _mask)}
  69. #define EMH_EMPTY(i, n) (0 > (int)i[n].bucket)
  70. namespace emhash8 {
  71. #ifndef EMH_DEFAULT_LOAD_FACTOR
  72. constexpr static float EMH_DEFAULT_LOAD_FACTOR = 0.80f;
  73. constexpr static float EMH_MIN_LOAD_FACTOR = 0.25f; //< 0.5
  74. #endif
  75. #if EMH_CACHE_LINE_SIZE < 32
  76. constexpr static uint32_t EMH_CACHE_LINE_SIZE = 64;
  77. #endif
  78. template <typename KeyT, typename ValueT, typename HashT = std::hash<KeyT>, typename EqT = std::equal_to<KeyT>>
  79. class HashMap
  80. {
  81. public:
  82. using htype = HashMap<KeyT, ValueT, HashT, EqT>;
  83. using value_type = std::pair<KeyT, ValueT>;
  84. using key_type = KeyT;
  85. using mapped_type = ValueT;
  86. #ifdef EMH_SMALL_TYPE
  87. using size_type = uint16_t;
  88. #elif EMH_SIZE_TYPE == 0
  89. using size_type = uint32_t;
  90. #else
  91. using size_type = size_t;
  92. #endif
  93. using hasher = HashT;
  94. using key_equal = EqT;
  95. constexpr static size_type INACTIVE = 0-1u;
  96. //constexpr uint32_t END = 0-0x1u;
  97. constexpr static size_type EAD = 2;
  98. struct Index
  99. {
  100. size_type bucket;
  101. size_type slot;
  102. };
  103. class const_iterator;
  104. class iterator
  105. {
  106. public:
  107. using iterator_category = std::bidirectional_iterator_tag;
  108. using difference_type = std::ptrdiff_t;
  109. using value_type = typename htype::value_type;
  110. using pointer = value_type*;
  111. using const_pointer = const value_type* ;
  112. using reference = value_type&;
  113. using const_reference = const value_type&;
  114. iterator() : kv_(nullptr) {}
  115. iterator(const_iterator& cit) {
  116. kv_ = cit.kv_;
  117. }
  118. iterator(const htype* hash_map, size_type bucket) {
  119. kv_ = hash_map->_pairs + (int)bucket;
  120. }
  121. iterator& operator++()
  122. {
  123. kv_ ++;
  124. return *this;
  125. }
  126. iterator operator++(int)
  127. {
  128. auto cur = *this; kv_ ++;
  129. return cur;
  130. }
  131. iterator& operator--()
  132. {
  133. kv_ --;
  134. return *this;
  135. }
  136. iterator operator--(int)
  137. {
  138. auto cur = *this; kv_ --;
  139. return cur;
  140. }
  141. reference operator*() const { return *kv_; }
  142. pointer operator->() const { return kv_; }
  143. bool operator == (const iterator& rhs) const { return kv_ == rhs.kv_; }
  144. bool operator != (const iterator& rhs) const { return kv_ != rhs.kv_; }
  145. bool operator == (const const_iterator& rhs) const { return kv_ == rhs.kv_; }
  146. bool operator != (const const_iterator& rhs) const { return kv_ != rhs.kv_; }
  147. public:
  148. value_type* kv_;
  149. };
  150. class const_iterator
  151. {
  152. public:
  153. using iterator_category = std::bidirectional_iterator_tag;
  154. using value_type = typename htype::value_type;
  155. using difference_type = std::ptrdiff_t;
  156. using pointer = value_type*;
  157. using const_pointer = const value_type*;
  158. using reference = value_type&;
  159. using const_reference = const value_type&;
  160. const_iterator(const iterator& it) {
  161. kv_ = it.kv_;
  162. }
  163. const_iterator (const htype* hash_map, size_type bucket) {
  164. kv_ = hash_map->_pairs + (int)bucket;
  165. }
  166. const_iterator& operator++()
  167. {
  168. kv_ ++;
  169. return *this;
  170. }
  171. const_iterator operator++(int)
  172. {
  173. auto cur = *this; kv_ ++;
  174. return cur;
  175. }
  176. const_iterator& operator--()
  177. {
  178. kv_ --;
  179. return *this;
  180. }
  181. const_iterator operator--(int)
  182. {
  183. auto cur = *this; kv_ --;
  184. return cur;
  185. }
  186. const_reference operator*() const { return *kv_; }
  187. const_pointer operator->() const { return kv_; }
  188. bool operator == (const iterator& rhs) const { return kv_ == rhs.kv_; }
  189. bool operator != (const iterator& rhs) const { return kv_ != rhs.kv_; }
  190. bool operator == (const const_iterator& rhs) const { return kv_ == rhs.kv_; }
  191. bool operator != (const const_iterator& rhs) const { return kv_ != rhs.kv_; }
  192. public:
  193. const value_type* kv_;
  194. };
  195. void init(size_type bucket, float mlf = EMH_DEFAULT_LOAD_FACTOR)
  196. {
  197. _pairs = nullptr;
  198. _index = nullptr;
  199. _mask = _num_buckets = 0;
  200. _num_filled = 0;
  201. max_load_factor(mlf);
  202. rehash(bucket);
  203. }
  204. HashMap(size_type bucket = 2, float mlf = EMH_DEFAULT_LOAD_FACTOR)
  205. {
  206. init(bucket, mlf);
  207. }
  208. HashMap(const HashMap& rhs)
  209. {
  210. if (rhs.load_factor() > EMH_MIN_LOAD_FACTOR) {
  211. _pairs = alloc_bucket((size_type)(rhs._num_buckets * rhs.max_load_factor()) + 4);
  212. _index = alloc_index(rhs._num_buckets);
  213. clone(rhs);
  214. } else {
  215. init(rhs._num_filled + 2, EMH_DEFAULT_LOAD_FACTOR);
  216. for (auto it = rhs.begin(); it != rhs.end(); ++it)
  217. insert_unique(it->first, it->second);
  218. }
  219. }
  220. HashMap(HashMap&& rhs) noexcept
  221. {
  222. init(0);
  223. *this = std::move(rhs);
  224. }
  225. HashMap(std::initializer_list<value_type> ilist)
  226. {
  227. init((size_type)ilist.size());
  228. for (auto it = ilist.begin(); it != ilist.end(); ++it)
  229. do_insert(*it);
  230. }
  231. template<class InputIt>
  232. HashMap(InputIt first, InputIt last, size_type bucket_count=4)
  233. {
  234. init(std::distance(first, last) + bucket_count);
  235. for (; first != last; ++first)
  236. emplace(*first);
  237. }
  238. HashMap& operator=(const HashMap& rhs)
  239. {
  240. if (this == &rhs)
  241. return *this;
  242. if (rhs.load_factor() < EMH_MIN_LOAD_FACTOR) {
  243. clear(); free(_pairs); _pairs = nullptr;
  244. rehash(rhs._num_filled + 2);
  245. for (auto it = rhs.begin(); it != rhs.end(); ++it)
  246. insert_unique(it->first, it->second);
  247. return *this;
  248. }
  249. clearkv();
  250. if (_num_buckets != rhs._num_buckets) {
  251. free(_pairs); free(_index);
  252. _index = alloc_index(rhs._num_buckets);
  253. _pairs = alloc_bucket((size_type)(rhs._num_buckets * rhs.max_load_factor()) + 4);
  254. }
  255. clone(rhs);
  256. return *this;
  257. }
  258. HashMap& operator=(HashMap&& rhs) noexcept
  259. {
  260. if (this != &rhs) {
  261. swap(rhs);
  262. rhs.clear();
  263. }
  264. return *this;
  265. }
  266. template<typename Con>
  267. bool operator == (const Con& rhs) const
  268. {
  269. if (size() != rhs.size())
  270. return false;
  271. for (auto it = begin(), last = end(); it != last; ++it) {
  272. auto oi = rhs.find(it->first);
  273. if (oi == rhs.end() || it->second != oi->second)
  274. return false;
  275. }
  276. return true;
  277. }
  278. template<typename Con>
  279. bool operator != (const Con& rhs) const { return !(*this == rhs); }
  280. ~HashMap() noexcept
  281. {
  282. clearkv();
  283. free(_pairs);
  284. free(_index);
  285. }
  286. void clone(const HashMap& rhs)
  287. {
  288. _hasher = rhs._hasher;
  289. // _eq = rhs._eq;
  290. _num_buckets = rhs._num_buckets;
  291. _num_filled = rhs._num_filled;
  292. _mlf = rhs._mlf;
  293. _last = rhs._last;
  294. _mask = rhs._mask;
  295. #if EMH_HIGH_LOAD
  296. _ehead = rhs._ehead;
  297. #endif
  298. _etail = rhs._etail;
  299. auto opairs = rhs._pairs;
  300. memcpy((char*)_index, (char*)rhs._index, (_num_buckets + EAD) * sizeof(Index));
  301. if (is_copy_trivially()) {
  302. if (opairs)
  303. memcpy((char*)_pairs, (char*)opairs, _num_filled * sizeof(value_type));
  304. } else {
  305. for (size_type slot = 0; slot < _num_filled; slot++)
  306. new(_pairs + slot) value_type(opairs[slot]);
  307. }
  308. }
  309. void swap(HashMap& rhs)
  310. {
  311. // std::swap(_eq, rhs._eq);
  312. std::swap(_hasher, rhs._hasher);
  313. std::swap(_pairs, rhs._pairs);
  314. std::swap(_index, rhs._index);
  315. std::swap(_num_buckets, rhs._num_buckets);
  316. std::swap(_num_filled, rhs._num_filled);
  317. std::swap(_mask, rhs._mask);
  318. std::swap(_mlf, rhs._mlf);
  319. std::swap(_last, rhs._last);
  320. #if EMH_HIGH_LOAD
  321. std::swap(_ehead, rhs._ehead);
  322. #endif
  323. std::swap(_etail, rhs._etail);
  324. }
  325. // -------------------------------------------------------------
  326. inline iterator first() const { return {this, 0}; }
  327. inline iterator last() const { return {this, _num_filled - 1}; }
  328. inline iterator begin() { return first(); }
  329. inline const_iterator cbegin() const { return first(); }
  330. inline const_iterator begin() const { return first(); }
  331. inline iterator end() { return {this, _num_filled}; }
  332. inline const_iterator cend() const { return {this, _num_filled}; }
  333. inline const_iterator end() const { return cend(); }
  334. inline const value_type* values() const { return _pairs; }
  335. inline const Index* index() const { return _index; }
  336. inline size_type size() const { return _num_filled; }
  337. inline bool empty() const { return _num_filled == 0; }
  338. inline size_type bucket_count() const { return _num_buckets; }
  339. /// Returns average number of elements per bucket.
  340. inline float load_factor() const { return static_cast<float>(_num_filled) / (_mask + 1); }
  341. inline HashT& hash_function() const { return _hasher; }
  342. inline EqT& key_eq() const { return _eq; }
  343. void max_load_factor(float mlf)
  344. {
  345. if (mlf < 0.991 && mlf > EMH_MIN_LOAD_FACTOR) {
  346. _mlf = (uint32_t)((1 << 27) / mlf);
  347. if (_num_buckets > 0) rehash(_num_buckets);
  348. }
  349. }
  350. inline constexpr float max_load_factor() const { return (1 << 27) / (float)_mlf; }
  351. inline constexpr size_type max_size() const { return (1ull << (sizeof(size_type) * 8 - 1)); }
  352. inline constexpr size_type max_bucket_count() const { return max_size(); }
  353. #if EMH_STATIS
  354. //Returns the bucket number where the element with key k is located.
  355. size_type bucket(const KeyT& key) const
  356. {
  357. const auto bucket = hash_bucket(key);
  358. const auto next_bucket = EMH_BUCKET(_index, bucket);
  359. if ((int)next_bucket < 0)
  360. return 0;
  361. else if (bucket == next_bucket)
  362. return bucket + 1;
  363. return hash_main(bucket) + 1;
  364. }
  365. //Returns the number of elements in bucket n.
  366. size_type bucket_size(const size_type bucket) const
  367. {
  368. auto next_bucket = EMH_BUCKET(_index, bucket);
  369. if ((int)next_bucket < 0)
  370. return 0;
  371. next_bucket = hash_main(bucket);
  372. size_type ibucket_size = 1;
  373. //iterator each item in current main bucket
  374. while (true) {
  375. const auto nbucket = EMH_BUCKET(_index, next_bucket);
  376. if (nbucket == next_bucket) {
  377. break;
  378. }
  379. ibucket_size ++;
  380. next_bucket = nbucket;
  381. }
  382. return ibucket_size;
  383. }
  384. size_type get_main_bucket(const size_type bucket) const
  385. {
  386. auto next_bucket = EMH_BUCKET(_index, bucket);
  387. if ((int)next_bucket < 0)
  388. return INACTIVE;
  389. return hash_main(bucket);
  390. }
  391. size_type get_diss(size_type bucket, size_type next_bucket, const size_type slots) const
  392. {
  393. auto pbucket = reinterpret_cast<uint64_t>(&_pairs[bucket]);
  394. auto pnext = reinterpret_cast<uint64_t>(&_pairs[next_bucket]);
  395. if (pbucket / EMH_CACHE_LINE_SIZE == pnext / EMH_CACHE_LINE_SIZE)
  396. return 0;
  397. size_type diff = pbucket > pnext ? (pbucket - pnext) : (pnext - pbucket);
  398. if (diff / EMH_CACHE_LINE_SIZE < slots - 1)
  399. return diff / EMH_CACHE_LINE_SIZE + 1;
  400. return slots - 1;
  401. }
  402. int get_bucket_info(const size_type bucket, size_type steps[], const size_type slots) const
  403. {
  404. auto next_bucket = EMH_BUCKET(_index, bucket);
  405. if ((int)next_bucket < 0)
  406. return -1;
  407. const auto main_bucket = hash_main(bucket);
  408. if (next_bucket == main_bucket)
  409. return 1;
  410. else if (main_bucket != bucket)
  411. return 0;
  412. steps[get_diss(bucket, next_bucket, slots)] ++;
  413. size_type ibucket_size = 2;
  414. //find a empty and linked it to tail
  415. while (true) {
  416. const auto nbucket = EMH_BUCKET(_index, next_bucket);
  417. if (nbucket == next_bucket)
  418. break;
  419. steps[get_diss(nbucket, next_bucket, slots)] ++;
  420. ibucket_size ++;
  421. next_bucket = nbucket;
  422. }
  423. return (int)ibucket_size;
  424. }
  425. void dump_statics() const
  426. {
  427. const size_type slots = 128;
  428. size_type buckets[slots + 1] = {0};
  429. size_type steps[slots + 1] = {0};
  430. for (size_type bucket = 0; bucket < _num_buckets; ++bucket) {
  431. auto bsize = get_bucket_info(bucket, steps, slots);
  432. if (bsize > 0)
  433. buckets[bsize] ++;
  434. }
  435. size_type sumb = 0, collision = 0, sumc = 0, finds = 0, sumn = 0;
  436. puts("============== buckets size ration =========");
  437. for (size_type i = 0; i < sizeof(buckets) / sizeof(buckets[0]); i++) {
  438. const auto bucketsi = buckets[i];
  439. if (bucketsi == 0)
  440. continue;
  441. sumb += bucketsi;
  442. sumn += bucketsi * i;
  443. collision += bucketsi * (i - 1);
  444. finds += bucketsi * i * (i + 1) / 2;
  445. printf(" %2u %8u %2.2lf| %.2lf\n", i, bucketsi, bucketsi * 100.0 * i / _num_filled, sumn * 100.0 / _num_filled);
  446. }
  447. puts("========== collision miss ration ===========");
  448. for (size_type i = 0; i < sizeof(steps) / sizeof(steps[0]); i++) {
  449. sumc += steps[i];
  450. if (steps[i] <= 2)
  451. continue;
  452. printf(" %2u %8u %.2lf %.2lf\n", i, steps[i], steps[i] * 100.0 / collision, sumc * 100.0 / collision);
  453. }
  454. if (sumb == 0) return;
  455. printf(" _num_filled/bucket_size/packed collision/cache_miss/hit_find = %u/%.2lf/%zd/ %.2lf%%/%.2lf%%/%.2lf\n",
  456. _num_filled, _num_filled * 1.0 / sumb, sizeof(value_type), (collision * 100.0 / _num_filled), (collision - steps[0]) * 100.0 / _num_filled, finds * 1.0 / _num_filled);
  457. assert(sumn == _num_filled);
  458. assert(sumc == collision);
  459. puts("============== buckets size end =============");
  460. }
  461. #endif
  462. // ------------------------------------------------------------
  463. template<typename K=KeyT>
  464. inline iterator find(const K& key) noexcept
  465. {
  466. return {this, find_filled_slot(key)};
  467. }
  468. template<typename K=KeyT>
  469. inline const_iterator find(const K& key) const noexcept
  470. {
  471. return {this, find_filled_slot(key)};
  472. }
  473. template<typename K=KeyT>
  474. ValueT& at(const K& key)
  475. {
  476. const auto slot = find_filled_slot(key);
  477. //throw
  478. return EMH_VAL(_pairs, slot);
  479. }
  480. template<typename K=KeyT>
  481. const ValueT& at(const K& key) const
  482. {
  483. const auto slot = find_filled_slot(key);
  484. //throw
  485. return EMH_VAL(_pairs, slot);
  486. }
  487. template<typename K=KeyT>
  488. inline bool contains(const K& key) const noexcept
  489. {
  490. return find_filled_slot(key) != _num_filled;
  491. }
  492. template<typename K=KeyT>
  493. inline size_type count(const K& key) const noexcept
  494. {
  495. return find_filled_slot(key) == _num_filled ? 0 : 1;
  496. //return find_sorted_bucket(key) == END ? 0 : 1;
  497. //return find_hash_bucket(key) == END ? 0 : 1;
  498. }
  499. template<typename K=KeyT>
  500. std::pair<iterator, iterator> equal_range(const K& key)
  501. {
  502. const auto found = find(key);
  503. if (found.second == _num_filled)
  504. return { found, found };
  505. else
  506. return { found, std::next(found) };
  507. }
  508. void merge(HashMap& rhs)
  509. {
  510. if (empty()) {
  511. *this = std::move(rhs);
  512. return;
  513. }
  514. for (auto rit = rhs.begin(); rit != rhs.end(); ) {
  515. auto fit = find(rit->first);
  516. if (fit == end()) {
  517. insert_unique(rit->first, std::move(rit->second));
  518. rit = rhs.erase(rit);
  519. } else {
  520. ++rit;
  521. }
  522. }
  523. }
  524. /// Returns the matching ValueT or nullptr if k isn't found.
  525. bool try_get(const KeyT& key, ValueT& val) const noexcept
  526. {
  527. const auto slot = find_filled_slot(key);
  528. const auto found = slot != _num_filled;
  529. if (found) {
  530. val = EMH_VAL(_pairs, slot);
  531. }
  532. return found;
  533. }
  534. /// Returns the matching ValueT or nullptr if k isn't found.
  535. ValueT* try_get(const KeyT& key) noexcept
  536. {
  537. const auto slot = find_filled_slot(key);
  538. return slot != _num_filled ? &EMH_VAL(_pairs, slot) : nullptr;
  539. }
  540. /// Const version of the above
  541. ValueT* try_get(const KeyT& key) const noexcept
  542. {
  543. const auto slot = find_filled_slot(key);
  544. return slot != _num_filled ? &EMH_VAL(_pairs, slot) : nullptr;
  545. }
  546. /// set value if key exist
  547. bool try_set(const KeyT& key, const ValueT& val) noexcept
  548. {
  549. const auto slot = find_filled_slot(key);
  550. if (slot == _num_filled)
  551. return false;
  552. EMH_VAL(_pairs, slot) = val;
  553. return true;
  554. }
  555. /// set value if key exist
  556. bool try_set(const KeyT& key, ValueT&& val) noexcept
  557. {
  558. const auto slot = find_filled_slot(key);
  559. if (slot == _num_filled)
  560. return false;
  561. EMH_VAL(_pairs, slot) = std::move(val);
  562. return true;
  563. }
  564. /// Convenience function.
  565. ValueT get_or_return_default(const KeyT& key) const noexcept
  566. {
  567. const auto slot = find_filled_slot(key);
  568. return slot == _num_filled ? ValueT() : EMH_VAL(_pairs, slot);
  569. }
  570. // -----------------------------------------------------
  571. std::pair<iterator, bool> do_insert(const value_type& value) noexcept
  572. {
  573. const auto key_hash = hash_key(value.first);
  574. const auto bucket = find_or_allocate(value.first, key_hash);
  575. const auto bempty = EMH_EMPTY(_index, bucket);
  576. if (bempty) {
  577. EMH_NEW(value.first, value.second, bucket, key_hash);
  578. }
  579. const auto slot = EMH_SLOT(_index, bucket);
  580. return { {this, slot}, bempty };
  581. }
  582. std::pair<iterator, bool> do_insert(value_type&& value) noexcept
  583. {
  584. const auto key_hash = hash_key(value.first);
  585. const auto bucket = find_or_allocate(value.first, key_hash);
  586. const auto bempty = EMH_EMPTY(_index, bucket);
  587. if (bempty) {
  588. EMH_NEW(std::move(value.first), std::move(value.second), bucket, key_hash);
  589. }
  590. const auto slot = EMH_SLOT(_index, bucket);
  591. return { {this, slot}, bempty };
  592. }
  593. template<typename K, typename V>
  594. std::pair<iterator, bool> do_insert(K&& key, V&& val) noexcept
  595. {
  596. const auto key_hash = hash_key(key);
  597. const auto bucket = find_or_allocate(key, key_hash);
  598. const auto bempty = EMH_EMPTY(_index, bucket);
  599. if (bempty) {
  600. EMH_NEW(std::forward<K>(key), std::forward<V>(val), bucket, key_hash);
  601. }
  602. const auto slot = EMH_SLOT(_index, bucket);
  603. return { {this, slot}, bempty };
  604. }
  605. template<typename K, typename V>
  606. std::pair<iterator, bool> do_assign(K&& key, V&& val) noexcept
  607. {
  608. check_expand_need();
  609. const auto key_hash = hash_key(key);
  610. const auto bucket = find_or_allocate(key, key_hash);
  611. const auto bempty = EMH_EMPTY(_index, bucket);
  612. if (bempty) {
  613. EMH_NEW(std::forward<K>(key), std::forward<V>(val), bucket, key_hash);
  614. } else {
  615. EMH_VAL(_pairs, EMH_SLOT(_index, bucket)) = std::move(val);
  616. }
  617. const auto slot = EMH_SLOT(_index, bucket);
  618. return { {this, slot}, bempty };
  619. }
  620. std::pair<iterator, bool> insert(const value_type& p)
  621. {
  622. check_expand_need();
  623. return do_insert(p);
  624. }
  625. std::pair<iterator, bool> insert(value_type && p)
  626. {
  627. check_expand_need();
  628. return do_insert(std::move(p));
  629. }
  630. void insert(std::initializer_list<value_type> ilist)
  631. {
  632. reserve(ilist.size() + _num_filled, false);
  633. for (auto it = ilist.begin(); it != ilist.end(); ++it)
  634. do_insert(*it);
  635. }
  636. template <typename Iter>
  637. void insert(Iter first, Iter last)
  638. {
  639. reserve(std::distance(first, last) + _num_filled, false);
  640. for (; first != last; ++first)
  641. do_insert(first->first, first->second);
  642. }
  643. #if 0
  644. template <typename Iter>
  645. void insert_unique(Iter begin, Iter end)
  646. {
  647. reserve(std::distance(begin, end) + _num_filled, false);
  648. for (; begin != end; ++begin) {
  649. insert_unique(*begin);
  650. }
  651. }
  652. #endif
  653. template<typename K, typename V>
  654. size_type insert_unique(K&& key, V&& val)
  655. {
  656. check_expand_need();
  657. const auto key_hash = hash_key(key);
  658. auto bucket = find_unique_bucket(key_hash);
  659. EMH_NEW(std::forward<K>(key), std::forward<V>(val), bucket, key_hash);
  660. return bucket;
  661. }
  662. size_type insert_unique(value_type&& value)
  663. {
  664. return insert_unique(std::move(value.first), std::move(value.second));
  665. }
  666. inline size_type insert_unique(const value_type& value)
  667. {
  668. return insert_unique(value.first, value.second);
  669. }
  670. template <class... Args>
  671. inline std::pair<iterator, bool> emplace(Args&&... args) noexcept
  672. {
  673. check_expand_need();
  674. return do_insert(std::forward<Args>(args)...);
  675. }
  676. //no any optimize for position
  677. template <class... Args>
  678. iterator emplace_hint(const_iterator hint, Args&&... args)
  679. {
  680. (void)hint;
  681. check_expand_need();
  682. return do_insert(std::forward<Args>(args)...).first;
  683. }
  684. template<class... Args>
  685. std::pair<iterator, bool> try_emplace(const KeyT& k, Args&&... args)
  686. {
  687. check_expand_need();
  688. return do_insert(k, std::forward<Args>(args)...);
  689. }
  690. template<class... Args>
  691. std::pair<iterator, bool> try_emplace(KeyT&& k, Args&&... args)
  692. {
  693. check_expand_need();
  694. return do_insert(std::move(k), std::forward<Args>(args)...);
  695. }
  696. template <class... Args>
  697. inline size_type emplace_unique(Args&&... args)
  698. {
  699. return insert_unique(std::forward<Args>(args)...);
  700. }
  701. std::pair<iterator, bool> insert_or_assign(const KeyT& key, ValueT&& val) { return do_assign(key, std::forward<ValueT>(val)); }
  702. std::pair<iterator, bool> insert_or_assign(KeyT&& key, ValueT&& val) { return do_assign(std::move(key), std::forward<ValueT>(val)); }
  703. /// Return the old value or ValueT() if it didn't exist.
  704. ValueT set_get(const KeyT& key, const ValueT& val)
  705. {
  706. check_expand_need();
  707. const auto key_hash = hash_key(key);
  708. const auto bucket = find_or_allocate(key, key_hash);
  709. if (EMH_EMPTY(_index, bucket)) {
  710. EMH_NEW(key, val, bucket, key_hash);
  711. return ValueT();
  712. } else {
  713. const auto slot = EMH_SLOT(_index, bucket);
  714. ValueT old_value(val);
  715. std::swap(EMH_VAL(_pairs, slot), old_value);
  716. return old_value;
  717. }
  718. }
  719. /// Like std::map<KeyT, ValueT>::operator[].
  720. ValueT& operator[](const KeyT& key) noexcept
  721. {
  722. check_expand_need();
  723. const auto key_hash = hash_key(key);
  724. const auto bucket = find_or_allocate(key, key_hash);
  725. if (EMH_EMPTY(_index, bucket)) {
  726. /* Check if inserting a value rather than overwriting an old entry */
  727. EMH_NEW(key, std::move(ValueT()), bucket, key_hash);
  728. }
  729. const auto slot = EMH_SLOT(_index, bucket);
  730. return EMH_VAL(_pairs, slot);
  731. }
  732. ValueT& operator[](KeyT&& key) noexcept
  733. {
  734. check_expand_need();
  735. const auto key_hash = hash_key(key);
  736. const auto bucket = find_or_allocate(key, key_hash);
  737. if (EMH_EMPTY(_index, bucket)) {
  738. EMH_NEW(std::move(key), std::move(ValueT()), bucket, key_hash);
  739. }
  740. const auto slot = EMH_SLOT(_index, bucket);
  741. return EMH_VAL(_pairs, slot);
  742. }
  743. /// Erase an element from the hash table.
  744. /// return 0 if element was not found
  745. size_type erase(const KeyT& key) noexcept
  746. {
  747. const auto key_hash = hash_key(key);
  748. const auto sbucket = find_filled_bucket(key, key_hash);
  749. if (sbucket == INACTIVE)
  750. return 0;
  751. const auto main_bucket = key_hash & _mask;
  752. erase_slot(sbucket, (size_type)main_bucket);
  753. return 1;
  754. }
  755. //iterator erase(const_iterator begin_it, const_iterator end_it)
  756. iterator erase(const const_iterator& cit) noexcept
  757. {
  758. const auto slot = (size_type)(cit.kv_ - _pairs);
  759. size_type main_bucket;
  760. const auto sbucket = find_slot_bucket(slot, main_bucket); //TODO
  761. erase_slot(sbucket, main_bucket);
  762. return {this, slot};
  763. }
  764. //only last >= first
  765. iterator erase(const_iterator first, const_iterator last) noexcept
  766. {
  767. auto esize = long(last.kv_ - first.kv_);
  768. auto tsize = long((_pairs + _num_filled) - last.kv_); //last to tail size
  769. auto next = first;
  770. while (tsize -- > 0) {
  771. if (esize-- <= 0)
  772. break;
  773. next = ++erase(next);
  774. }
  775. //fast erase from last
  776. next = this->last();
  777. while (esize -- > 0)
  778. next = --erase(next);
  779. return {this, size_type(next.kv_ - _pairs)};
  780. }
  781. template<typename Pred>
  782. size_type erase_if(Pred pred)
  783. {
  784. auto old_size = size();
  785. for (auto it = begin(); it != end();) {
  786. if (pred(*it))
  787. it = erase(it);
  788. else
  789. ++it;
  790. }
  791. return old_size - size();
  792. }
  793. static constexpr bool is_triviall_destructable()
  794. {
  795. #if __cplusplus >= 201402L || _MSC_VER > 1600
  796. return !(std::is_trivially_destructible<KeyT>::value && std::is_trivially_destructible<ValueT>::value);
  797. #else
  798. return !(std::is_pod<KeyT>::value && std::is_pod<ValueT>::value);
  799. #endif
  800. }
  801. static constexpr bool is_copy_trivially()
  802. {
  803. #if __cplusplus >= 201103L || _MSC_VER > 1600
  804. return (std::is_trivially_copyable<KeyT>::value && std::is_trivially_copyable<ValueT>::value);
  805. #else
  806. return (std::is_pod<KeyT>::value && std::is_pod<ValueT>::value);
  807. #endif
  808. }
  809. void clearkv()
  810. {
  811. if (is_triviall_destructable()) {
  812. while (_num_filled --)
  813. _pairs[_num_filled].~value_type();
  814. }
  815. }
  816. /// Remove all elements, keeping full capacity.
  817. void clear() noexcept
  818. {
  819. clearkv();
  820. if (_num_filled > 0)
  821. memset((char*)_index, INACTIVE, sizeof(_index[0]) * _num_buckets);
  822. _last = _num_filled = 0;
  823. _etail = INACTIVE;
  824. #if EMH_HIGH_LOAD
  825. _ehead = 0;
  826. #endif
  827. }
  828. void shrink_to_fit(const float min_factor = EMH_DEFAULT_LOAD_FACTOR / 4)
  829. {
  830. if (load_factor() < min_factor && bucket_count() > 10) //safe guard
  831. rehash(_num_filled + 1);
  832. }
  833. #if EMH_HIGH_LOAD
  834. void set_empty()
  835. {
  836. auto prev = 0;
  837. for (int32_t bucket = 1; bucket < _num_buckets; ++bucket) {
  838. if (EMH_EMPTY(_index, bucket)) {
  839. if (prev != 0) {
  840. EMH_PREVET(_index, bucket) = prev;
  841. EMH_BUCKET(_index, prev) = -bucket;
  842. }
  843. else
  844. _ehead = bucket;
  845. prev = bucket;
  846. }
  847. }
  848. EMH_PREVET(_index, _ehead) = prev;
  849. EMH_BUCKET(_index, prev) = 0-_ehead;
  850. _ehead = 0-EMH_BUCKET(_index, _ehead);
  851. }
  852. void clear_empty()
  853. {
  854. auto prev = EMH_PREVET(_index, _ehead);
  855. while (prev != _ehead) {
  856. EMH_BUCKET(_index, prev) = INACTIVE;
  857. prev = EMH_PREVET(_index, prev);
  858. }
  859. EMH_BUCKET(_index, _ehead) = INACTIVE;
  860. _ehead = 0;
  861. }
  862. //prev-ehead->next
  863. size_type pop_empty(const size_type bucket)
  864. {
  865. const auto prev_bucket = EMH_PREVET(_index, bucket);
  866. const int next_bucket = 0-EMH_BUCKET(_index, bucket);
  867. EMH_PREVET(_index, next_bucket) = prev_bucket;
  868. EMH_BUCKET(_index, prev_bucket) = -next_bucket;
  869. _ehead = next_bucket;
  870. return bucket;
  871. }
  872. //ehead->bucket->next
  873. void push_empty(const int32_t bucket)
  874. {
  875. const int next_bucket = 0-EMH_BUCKET(_index, _ehead);
  876. assert(next_bucket > 0);
  877. EMH_PREVET(_index, bucket) = _ehead;
  878. EMH_BUCKET(_index, bucket) = -next_bucket;
  879. EMH_PREVET(_index, next_bucket) = bucket;
  880. EMH_BUCKET(_index, _ehead) = -bucket;
  881. // _ehead = bucket;
  882. }
  883. #endif
  884. /// Make room for this many elements
  885. bool reserve(uint64_t num_elems, bool force)
  886. {
  887. (void)force;
  888. #if EMH_HIGH_LOAD == 0
  889. const auto required_buckets = num_elems * _mlf >> 27;
  890. if (EMH_LIKELY(required_buckets < _mask)) // && !force
  891. return false;
  892. #elif EMH_HIGH_LOAD
  893. const auto required_buckets = num_elems + num_elems * 1 / 9;
  894. if (EMH_LIKELY(required_buckets < _mask))
  895. return false;
  896. else if (_num_buckets < 16 && _num_filled < _num_buckets)
  897. return false;
  898. else if (_num_buckets > EMH_HIGH_LOAD) {
  899. if (_ehead == 0) {
  900. set_empty();
  901. return false;
  902. } else if (/*_num_filled + 100 < _num_buckets && */EMH_BUCKET(_index, _ehead) != 0-_ehead) {
  903. return false;
  904. }
  905. }
  906. #endif
  907. #if EMH_STATIS
  908. if (_num_filled > EMH_STATIS) dump_statics();
  909. #endif
  910. //assert(required_buckets < max_size());
  911. rehash(required_buckets + 2);
  912. return true;
  913. }
  914. static value_type* alloc_bucket(size_type num_buckets)
  915. {
  916. auto new_pairs = (char*)malloc((uint64_t)num_buckets * sizeof(value_type));
  917. return (value_type *)(new_pairs);
  918. }
  919. static Index* alloc_index(size_type num_buckets)
  920. {
  921. auto new_index = (char*)malloc((uint64_t)(EAD + num_buckets) * sizeof(Index));
  922. return (Index *)(new_index);
  923. }
  924. bool reserve(size_type required_buckets) noexcept
  925. {
  926. if (_num_filled != required_buckets)
  927. return reserve(required_buckets, true);
  928. _last = 0;
  929. #if EMH_HIGH_LOAD
  930. _ehead = 0;
  931. #endif
  932. #if EMH_SORT
  933. std::sort(_pairs, _pairs + _num_filled, [this](const value_type & l, const value_type & r) {
  934. const auto hashl = (size_type)hash_key(l.first) & _mask, hashr = (size_type)hash_key(r.first) & _mask;
  935. return hashl < hashr;
  936. //return l.first < r.first;
  937. });
  938. #endif
  939. memset((char*)_index, INACTIVE, sizeof(_index[0]) * _num_buckets);
  940. for (size_type slot = 0; slot < _num_filled; slot++) {
  941. const auto& key = EMH_KEY(_pairs, slot);
  942. const auto key_hash = hash_key(key);
  943. const auto bucket = size_type(key_hash & _mask);
  944. auto& next_bucket = EMH_BUCKET(_index, bucket);
  945. if ((int)next_bucket < 0)
  946. EMH_INDEX(_index, bucket) = {1, slot | EMH_KEYMASK(key_hash, _mask)};
  947. else {
  948. EMH_HSLOT(_index, bucket) |= EMH_KEYMASK(key_hash, _mask);
  949. next_bucket ++;
  950. }
  951. }
  952. return true;
  953. }
  954. void rebuild(size_type num_buckets) noexcept
  955. {
  956. free(_index);
  957. auto new_pairs = (value_type*)alloc_bucket((size_type)(num_buckets * max_load_factor()) + 4);
  958. if (is_copy_trivially()) {
  959. memcpy((char*)new_pairs, (char*)_pairs, _num_filled * sizeof(value_type));
  960. } else {
  961. for (size_type slot = 0; slot < _num_filled; slot++) {
  962. new(new_pairs + slot) value_type(std::move(_pairs[slot]));
  963. if (is_triviall_destructable())
  964. _pairs[slot].~value_type();
  965. }
  966. }
  967. free(_pairs);
  968. _pairs = new_pairs;
  969. _index = (Index*)alloc_index (num_buckets);
  970. memset((char*)_index, INACTIVE, sizeof(_index[0]) * num_buckets);
  971. memset((char*)(_index + num_buckets), 0, sizeof(_index[0]) * EAD);
  972. }
  973. void rehash(uint64_t required_buckets)
  974. {
  975. if (required_buckets < _num_filled)
  976. return;
  977. assert(required_buckets < max_size());
  978. auto num_buckets = _num_filled > (1u << 16) ? (1u << 16) : 4u;
  979. while (num_buckets < required_buckets) { num_buckets *= 2; }
  980. #if EMH_REHASH_LOG
  981. auto last = _last;
  982. size_type collision = 0;
  983. #endif
  984. #if EMH_HIGH_LOAD
  985. _ehead = 0;
  986. #endif
  987. _last = _mask / 4;
  988. _mask = num_buckets - 1;
  989. #if EMH_PACK_TAIL > 1
  990. _last = _mask;
  991. num_buckets += num_buckets * EMH_PACK_TAIL / 100; //add more 5-10%
  992. #endif
  993. _num_buckets = num_buckets;
  994. rebuild(num_buckets);
  995. #ifdef EMH_SORT
  996. std::sort(_pairs, _pairs + _num_filled, [this](const value_type & l, const value_type & r) {
  997. const auto hashl = hash_key(l.first), hashr = hash_key(r.first);
  998. auto diff = int64_t((hashl & _mask) - (hashr & _mask));
  999. if (diff != 0)
  1000. return diff < 0;
  1001. return hashl < hashr;
  1002. // return l.first < r.first;
  1003. });
  1004. #endif
  1005. _etail = INACTIVE;
  1006. for (size_type slot = 0; slot < _num_filled; ++slot) {
  1007. const auto& key = EMH_KEY(_pairs, slot);
  1008. const auto key_hash = hash_key(key);
  1009. const auto bucket = find_unique_bucket(key_hash);
  1010. EMH_INDEX(_index, bucket) = {bucket, slot | EMH_KEYMASK(key_hash, _mask)};
  1011. #if EMH_REHASH_LOG
  1012. if (bucket != hash_main(bucket))
  1013. collision ++;
  1014. #endif
  1015. }
  1016. #if EMH_REHASH_LOG
  1017. if (_num_filled > EMH_REHASH_LOG) {
  1018. auto mbucket = _num_filled - collision;
  1019. char buff[255] = {0};
  1020. sprintf(buff, " _num_filled/aver_size/K.V/pack/collision|last = %u/%.2lf/%s.%s/%zd|%.2lf%%,%.2lf%%",
  1021. _num_filled, double (_num_filled) / mbucket, typeid(KeyT).name(), typeid(ValueT).name(), sizeof(_pairs[0]), collision * 100.0 / _num_filled, last * 100.0 / _num_buckets);
  1022. #ifdef EMH_LOG
  1023. static uint32_t ihashs = 0; EMH_LOG() << "hash_nums = " << ihashs ++ << "|" <<__FUNCTION__ << "|" << buff << endl;
  1024. #else
  1025. puts(buff);
  1026. #endif
  1027. }
  1028. #endif
  1029. }
  1030. private:
  1031. // Can we fit another element?
  1032. inline bool check_expand_need()
  1033. {
  1034. return reserve(_num_filled, false);
  1035. }
  1036. size_type slot_to_bucket(const size_type slot) const noexcept
  1037. {
  1038. size_type main_bucket;
  1039. return find_slot_bucket(slot, main_bucket); //TODO
  1040. }
  1041. //very slow
  1042. void erase_slot(const size_type sbucket, const size_type main_bucket) noexcept
  1043. {
  1044. const auto slot = EMH_SLOT(_index, sbucket);
  1045. const auto ebucket = erase_bucket(sbucket, main_bucket);
  1046. const auto last_slot = --_num_filled;
  1047. if (EMH_LIKELY(slot != last_slot)) {
  1048. const auto last_bucket = (_etail == INACTIVE || ebucket == _etail)
  1049. ? slot_to_bucket(last_slot) : _etail;
  1050. EMH_KV(_pairs, slot) = std::move(EMH_KV(_pairs, last_slot));
  1051. EMH_HSLOT(_index, last_bucket) = slot | (EMH_HSLOT(_index, last_bucket) & ~_mask);
  1052. }
  1053. if (is_triviall_destructable())
  1054. _pairs[last_slot].~value_type();
  1055. _etail = INACTIVE;
  1056. EMH_INDEX(_index, ebucket) = {INACTIVE, 0};
  1057. #if EMH_HIGH_LOAD
  1058. if (_ehead) {
  1059. if (10 * _num_filled < 8 * _num_buckets)
  1060. clear_empty();
  1061. else if (ebucket)
  1062. push_empty(ebucket);
  1063. }
  1064. #endif
  1065. }
  1066. size_type erase_bucket(const size_type bucket, const size_type main_bucket) noexcept
  1067. {
  1068. const auto next_bucket = EMH_BUCKET(_index, bucket);
  1069. if (bucket == main_bucket) {
  1070. if (main_bucket != next_bucket) {
  1071. const auto nbucket = EMH_BUCKET(_index, next_bucket);
  1072. EMH_INDEX(_index, main_bucket) = {
  1073. (nbucket == next_bucket) ? main_bucket : nbucket,
  1074. EMH_HSLOT(_index, next_bucket)
  1075. };
  1076. }
  1077. return next_bucket;
  1078. }
  1079. const auto prev_bucket = find_prev_bucket(main_bucket, bucket);
  1080. EMH_BUCKET(_index, prev_bucket) = (bucket == next_bucket) ? prev_bucket : next_bucket;
  1081. return bucket;
  1082. }
  1083. // Find the slot with this key, or return bucket size
  1084. size_type find_slot_bucket(const size_type slot, size_type& main_bucket) const
  1085. {
  1086. const auto key_hash = hash_key(EMH_KEY(_pairs, slot));
  1087. const auto bucket = main_bucket = size_type(key_hash & _mask);
  1088. if (slot == EMH_SLOT(_index, bucket))
  1089. return bucket;
  1090. auto next_bucket = EMH_BUCKET(_index, bucket);
  1091. while (true) {
  1092. if (EMH_LIKELY(slot == EMH_SLOT(_index, next_bucket)))
  1093. return next_bucket;
  1094. next_bucket = EMH_BUCKET(_index, next_bucket);
  1095. }
  1096. return INACTIVE;
  1097. }
  1098. // Find the slot with this key, or return bucket size
  1099. size_type find_filled_bucket(const KeyT& key, uint64_t key_hash) const noexcept
  1100. {
  1101. const auto bucket = size_type(key_hash & _mask);
  1102. auto next_bucket = EMH_BUCKET(_index, bucket);
  1103. if (EMH_UNLIKELY((int)next_bucket < 0))
  1104. return INACTIVE;
  1105. if (EMH_EQHASH(bucket, key_hash)) {
  1106. const auto slot = EMH_SLOT(_index, bucket);
  1107. if (EMH_LIKELY(_eq(key, EMH_KEY(_pairs, slot))))
  1108. return bucket;
  1109. }
  1110. if (next_bucket == bucket)
  1111. return INACTIVE;
  1112. while (true) {
  1113. if (EMH_EQHASH(next_bucket, key_hash)) {
  1114. const auto slot = EMH_SLOT(_index, next_bucket);
  1115. if (EMH_LIKELY(_eq(key, EMH_KEY(_pairs, slot))))
  1116. return next_bucket;
  1117. }
  1118. const auto nbucket = EMH_BUCKET(_index, next_bucket);
  1119. if (nbucket == next_bucket)
  1120. return INACTIVE;
  1121. next_bucket = nbucket;
  1122. }
  1123. return INACTIVE;
  1124. }
  1125. // Find the slot with this key, or return bucket size
  1126. template<typename K=KeyT>
  1127. size_type find_filled_slot(const K& key) const noexcept
  1128. {
  1129. const auto key_hash = hash_key(key);
  1130. const auto bucket = size_type(key_hash & _mask);
  1131. auto next_bucket = EMH_BUCKET(_index, bucket);
  1132. if ((int)next_bucket < 0)
  1133. return _num_filled;
  1134. if (EMH_EQHASH(bucket, key_hash)) {
  1135. const auto slot = EMH_SLOT(_index, bucket);
  1136. if (EMH_LIKELY(_eq(key, EMH_KEY(_pairs, slot))))
  1137. return slot;
  1138. }
  1139. if (next_bucket == bucket)
  1140. return _num_filled;
  1141. while (true) {
  1142. if (EMH_EQHASH(next_bucket, key_hash)) {
  1143. const auto slot = EMH_SLOT(_index, next_bucket);
  1144. if (EMH_LIKELY(_eq(key, EMH_KEY(_pairs, slot))))
  1145. return slot;
  1146. }
  1147. const auto nbucket = EMH_BUCKET(_index, next_bucket);
  1148. if (nbucket == next_bucket)
  1149. return _num_filled;
  1150. next_bucket = nbucket;
  1151. }
  1152. return _num_filled;
  1153. }
  1154. #if EMH_SORT
  1155. size_type find_hash_bucket(const KeyT& key) const noexcept
  1156. {
  1157. const auto key_hash = hash_key(key);
  1158. const auto bucket = size_type(key_hash & _mask);
  1159. const auto next_bucket = EMH_BUCKET(_index, bucket);
  1160. if ((int)next_bucket < 0)
  1161. return END;
  1162. auto slot = EMH_SLOT(_index, bucket);
  1163. if (_eq(key, EMH_KEY(_pairs, slot++)))
  1164. return slot;
  1165. else if (next_bucket == bucket)
  1166. return END;
  1167. while (true) {
  1168. const auto& okey = EMH_KEY(_pairs, slot++);
  1169. if (_eq(key, okey))
  1170. return slot;
  1171. const auto hasho = hash_key(okey);
  1172. if ((hasho & _mask) != bucket)
  1173. break;
  1174. else if (hasho > key_hash)
  1175. break;
  1176. else if (EMH_UNLIKELY(slot >= _num_filled))
  1177. break;
  1178. }
  1179. return END;
  1180. }
  1181. //only for find/can not insert
  1182. size_type find_sorted_bucket(const KeyT& key) const noexcept
  1183. {
  1184. const auto key_hash = hash_key(key);
  1185. const auto bucket = size_type(key_hash & _mask);
  1186. const auto slots = (int)(EMH_BUCKET(_index, bucket)); //TODO
  1187. if (slots < 0 /**|| key < EMH_KEY(_pairs, slot)*/)
  1188. return END;
  1189. const auto slot = EMH_SLOT(_index, bucket);
  1190. auto ormask = _index[bucket].slot & ~_mask;
  1191. auto hmask = EMH_KEYMASK(key_hash, _mask);
  1192. if ((hmask | ormask) != ormask)
  1193. return END;
  1194. if (_eq(key, EMH_KEY(_pairs, slot)))
  1195. return slot;
  1196. else if (slots == 1 || key < EMH_KEY(_pairs, slot))
  1197. return END;
  1198. #if EMH_SORT
  1199. if (key < EMH_KEY(_pairs, slot) || key > EMH_KEY(_pairs, slots + slot - 1))
  1200. return END;
  1201. #endif
  1202. for (size_type i = 1; i < slots; ++i) {
  1203. const auto& okey = EMH_KEY(_pairs, slot + i);
  1204. if (_eq(key, okey))
  1205. return slot + i;
  1206. // else if (okey > key)
  1207. // return END;
  1208. }
  1209. return END;
  1210. }
  1211. #endif
  1212. //kick out bucket and find empty to occpuy
  1213. //it will break the orgin link and relnik again.
  1214. //before: main_bucket-->prev_bucket --> bucket --> next_bucket
  1215. //atfer : main_bucket-->prev_bucket --> (removed)--> new_bucket--> next_bucket
  1216. size_type kickout_bucket(const size_type kmain, const size_type bucket) noexcept
  1217. {
  1218. const auto next_bucket = EMH_BUCKET(_index, bucket);
  1219. const auto new_bucket = find_empty_bucket(next_bucket, 2);
  1220. const auto prev_bucket = find_prev_bucket(kmain, bucket);
  1221. const auto last = next_bucket == bucket ? new_bucket : next_bucket;
  1222. EMH_INDEX(_index, new_bucket) = {last, EMH_HSLOT(_index, bucket)};
  1223. EMH_BUCKET(_index, prev_bucket) = new_bucket;
  1224. EMH_BUCKET(_index, bucket) = INACTIVE;
  1225. return bucket;
  1226. }
  1227. /*
  1228. ** inserts a new key into a hash table; first, check whether key's main
  1229. ** bucket/position is free. If not, check whether colliding node/bucket is in its main
  1230. ** position or not: if it is not, move colliding bucket to an empty place and
  1231. ** put new key in its main position; otherwise (colliding bucket is in its main
  1232. ** position), new key goes to an empty position.
  1233. */
  1234. template<typename K=KeyT>
  1235. size_type find_or_allocate(const K& key, uint64_t key_hash) noexcept
  1236. {
  1237. const auto bucket = size_type(key_hash & _mask);
  1238. auto next_bucket = EMH_BUCKET(_index, bucket);
  1239. if ((int)next_bucket < 0) {
  1240. #if EMH_HIGH_LOAD
  1241. if (next_bucket != INACTIVE)
  1242. pop_empty(bucket);
  1243. #endif
  1244. return bucket;
  1245. }
  1246. const auto slot = EMH_SLOT(_index, bucket);
  1247. if (EMH_EQHASH(bucket, key_hash))
  1248. if (EMH_LIKELY(_eq(key, EMH_KEY(_pairs, slot))))
  1249. return bucket;
  1250. //check current bucket_key is in main bucket or not
  1251. const auto kmain = hash_bucket(EMH_KEY(_pairs, slot));
  1252. if (kmain != bucket)
  1253. return kickout_bucket(kmain, bucket);
  1254. else if (next_bucket == bucket)
  1255. return EMH_BUCKET(_index, next_bucket) = find_empty_bucket(next_bucket, 1);
  1256. uint32_t csize = 1;
  1257. //find next linked bucket and check key
  1258. while (true) {
  1259. const auto eslot = EMH_SLOT(_index, next_bucket);
  1260. if (EMH_EQHASH(next_bucket, key_hash)) {
  1261. if (EMH_LIKELY(_eq(key, EMH_KEY(_pairs, eslot))))
  1262. return next_bucket;
  1263. }
  1264. csize += 1;
  1265. const auto nbucket = EMH_BUCKET(_index, next_bucket);
  1266. if (nbucket == next_bucket)
  1267. break;
  1268. next_bucket = nbucket;
  1269. }
  1270. //find a empty and link it to tail
  1271. const auto new_bucket = find_empty_bucket(next_bucket, csize);
  1272. return EMH_BUCKET(_index, next_bucket) = new_bucket;
  1273. }
  1274. size_type find_unique_bucket(uint64_t key_hash) noexcept
  1275. {
  1276. const auto bucket = size_type(key_hash & _mask);
  1277. auto next_bucket = EMH_BUCKET(_index, bucket);
  1278. if ((int)next_bucket < 0) {
  1279. #if EMH_HIGH_LOAD
  1280. if (next_bucket != INACTIVE)
  1281. pop_empty(bucket);
  1282. #endif
  1283. return bucket;
  1284. }
  1285. //check current bucket_key is in main bucket or not
  1286. const auto kmain = hash_main(bucket);
  1287. if (EMH_UNLIKELY(kmain != bucket))
  1288. return kickout_bucket(kmain, bucket);
  1289. else if (EMH_UNLIKELY(next_bucket != bucket))
  1290. next_bucket = find_last_bucket(next_bucket);
  1291. return EMH_BUCKET(_index, next_bucket) = find_empty_bucket(next_bucket, 2);
  1292. }
  1293. /***
  1294. Different probing techniques usually provide a trade-off between memory locality and avoidance of clustering.
  1295. Since Robin Hood hashing is relatively resilient to clustering (both primary and secondary), linear probing is the most cache friendly alternativeis typically used.
  1296. It's the core algorithm of this hash map with highly optimization/benchmark.
  1297. normaly linear probing is inefficient with high load factor, it use a new 3-way linear
  1298. probing strategy to search empty slot. from benchmark even the load factor > 0.9, it's more 2-3 timer fast than
  1299. one-way search strategy.
  1300. 1. linear or quadratic probing a few cache line for less cache miss from input slot "bucket_from".
  1301. 2. the first search slot from member variant "_last", init with 0
  1302. 3. the second search slot from calculated pos "(_num_filled + _last) & _mask", it's like a rand value
  1303. */
  1304. // key is not in this mavalue. Find a place to put it.
  1305. size_type find_empty_bucket(const size_type bucket_from, uint32_t csize) noexcept
  1306. {
  1307. #if EMH_HIGH_LOAD
  1308. if (_ehead)
  1309. return pop_empty(_ehead);
  1310. #endif
  1311. auto bucket = bucket_from;
  1312. if (EMH_EMPTY(_index, ++bucket) || EMH_EMPTY(_index, ++bucket))
  1313. return bucket;
  1314. #ifdef EMH_QUADRATIC
  1315. constexpr size_type linear_probe_length = 2 * EMH_CACHE_LINE_SIZE / sizeof(Index);//16
  1316. for (size_type offset = csize + 2, step = 4; offset <= linear_probe_length; ) {
  1317. bucket = (bucket_from + offset) & _mask;
  1318. if (EMH_EMPTY(_index, bucket) || EMH_EMPTY(_index, ++bucket))
  1319. return bucket;
  1320. offset += step; //7/8. 12. 16
  1321. }
  1322. #else
  1323. constexpr size_type quadratic_probe_length = 6u;
  1324. for (size_type offset = 4u, step = 3u; step < quadratic_probe_length; ) {
  1325. bucket = (bucket_from + offset) & _mask;
  1326. if (EMH_EMPTY(_index, bucket) || EMH_EMPTY(_index, ++bucket))
  1327. return bucket;
  1328. offset += step++;//3.4.5
  1329. }
  1330. #endif
  1331. #if EMH_PREFETCH
  1332. __builtin_prefetch(static_cast<const void*>(_index + _last + 1), 0, EMH_PREFETCH);
  1333. #endif
  1334. for (;;) {
  1335. #if EMH_PACK_TAIL
  1336. //find empty bucket and skip next
  1337. if (EMH_EMPTY(_index, _last++))// || EMH_EMPTY(_index, _last++))
  1338. return _last++ - 1;
  1339. if (EMH_UNLIKELY(_last >= _num_buckets))
  1340. _last = 0;
  1341. auto medium = (_mask / 4 + _last++) & _mask;
  1342. if (EMH_EMPTY(_index, medium))
  1343. return medium;
  1344. #else
  1345. if (EMH_EMPTY(_index, ++_last))// || EMH_EMPTY(_index, ++_last))
  1346. return _last++;
  1347. _last &= _mask;
  1348. auto medium = (_num_buckets / 2 + _last) & _mask;
  1349. if (EMH_EMPTY(_index, medium))// || EMH_EMPTY(_index, ++medium))
  1350. return _last = medium;
  1351. #endif
  1352. }
  1353. return 0;
  1354. }
  1355. size_type find_last_bucket(size_type main_bucket) const
  1356. {
  1357. auto next_bucket = EMH_BUCKET(_index, main_bucket);
  1358. if (next_bucket == main_bucket)
  1359. return main_bucket;
  1360. while (true) {
  1361. const auto nbucket = EMH_BUCKET(_index, next_bucket);
  1362. if (nbucket == next_bucket)
  1363. return next_bucket;
  1364. next_bucket = nbucket;
  1365. }
  1366. }
  1367. size_type find_prev_bucket(const size_type main_bucket, const size_type bucket) const
  1368. {
  1369. auto next_bucket = EMH_BUCKET(_index, main_bucket);
  1370. if (next_bucket == bucket)
  1371. return main_bucket;
  1372. while (true) {
  1373. const auto nbucket = EMH_BUCKET(_index, next_bucket);
  1374. if (nbucket == bucket)
  1375. return next_bucket;
  1376. next_bucket = nbucket;
  1377. }
  1378. }
  1379. inline size_type hash_bucket(const KeyT& key) const noexcept
  1380. {
  1381. return (size_type)hash_key(key) & _mask;
  1382. }
  1383. inline size_type hash_main(const size_type bucket) const noexcept
  1384. {
  1385. const auto slot = EMH_SLOT(_index, bucket);
  1386. return (size_type)hash_key(EMH_KEY(_pairs, slot)) & _mask;
  1387. }
  1388. #if EMH_INT_HASH
  1389. static constexpr uint64_t KC = UINT64_C(11400714819323198485);
  1390. static uint64_t hash64(uint64_t key)
  1391. {
  1392. #if __SIZEOF_INT128__ && EMH_INT_HASH == 1
  1393. __uint128_t r = key; r *= KC;
  1394. return (uint64_t)(r >> 64) + (uint64_t)r;
  1395. #elif EMH_INT_HASH == 2
  1396. //MurmurHash3Mixer
  1397. uint64_t h = key;
  1398. h ^= h >> 33;
  1399. h *= 0xff51afd7ed558ccd;
  1400. h ^= h >> 33;
  1401. h *= 0xc4ceb9fe1a85ec53;
  1402. h ^= h >> 33;
  1403. return h;
  1404. #elif _WIN64 && EMH_INT_HASH == 1
  1405. uint64_t high;
  1406. return _umul128(key, KC, &high) + high;
  1407. #elif EMH_INT_HASH == 3
  1408. auto ror = (key >> 32) | (key << 32);
  1409. auto low = key * 0xA24BAED4963EE407ull;
  1410. auto high = ror * 0x9FB21C651E98DF25ull;
  1411. auto mix = low + high;
  1412. return mix;
  1413. #elif EMH_INT_HASH == 1
  1414. uint64_t r = key * UINT64_C(0xca4bcaa75ec3f625);
  1415. return (r >> 32) + r;
  1416. #elif EMH_WYHASH64
  1417. return wyhash64(key, KC);
  1418. #else
  1419. uint64_t x = key;
  1420. x = (x ^ (x >> 30)) * UINT64_C(0xbf58476d1ce4e5b9);
  1421. x = (x ^ (x >> 27)) * UINT64_C(0x94d049bb133111eb);
  1422. x = x ^ (x >> 31);
  1423. return x;
  1424. #endif
  1425. }
  1426. #endif
  1427. #if EMH_WYHASH_HASH
  1428. //#define WYHASH_CONDOM 1
  1429. inline uint64_t wymix(uint64_t A, uint64_t B)
  1430. {
  1431. #if defined(__SIZEOF_INT128__)
  1432. __uint128_t r = A; r *= B;
  1433. #if WYHASH_CONDOM2
  1434. A ^= (uint64_t)r; B ^= (uint64_t)(r >> 64);
  1435. #else
  1436. A = (uint64_t)r; B = (uint64_t)(r >> 64);
  1437. #endif
  1438. #elif defined(_MSC_VER) && defined(_M_X64)
  1439. #if WYHASH_CONDOM2
  1440. uint64_t a, b;
  1441. a = _umul128(A, B, &b);
  1442. A ^= a; B ^= b;
  1443. #else
  1444. A = _umul128(A, B, &B);
  1445. #endif
  1446. #else
  1447. uint64_t ha = A >> 32, hb = B >> 32, la = (uint32_t)A, lb = (uint32_t)B, hi, lo;
  1448. uint64_t rh = ha * hb, rm0 = ha * lb, rm1 = hb * la, rl = la * lb, t = rl + (rm0 << 32), c = t < rl;
  1449. lo = t + (rm1 << 32); c += lo < t; hi = rh + (rm0 >> 32) + (rm1 >> 32) + c;
  1450. #if WYHASH_CONDOM2
  1451. A ^= lo; B ^= hi;
  1452. #else
  1453. A = lo; B = hi;
  1454. #endif
  1455. #endif
  1456. return A ^ B;
  1457. }
  1458. //multiply and xor mix function, aka MUM
  1459. static inline uint64_t wyr8(const uint8_t *p) { uint64_t v; memcpy(&v, p, 8); return v; }
  1460. static inline uint64_t wyr4(const uint8_t *p) { uint32_t v; memcpy(&v, p, 4); return v; }
  1461. static inline uint64_t wyr3(const uint8_t *p, size_t k) {
  1462. return (((uint64_t)p[0]) << 16) | (((uint64_t)p[k >> 1]) << 8) | p[k - 1];
  1463. }
  1464. static constexpr uint64_t secret[4] = {
  1465. 0xa0761d6478bd642full, 0xe7037ed1a0b428dbull,
  1466. 0x8ebc6af09c88c6e3ull, 0x589965cc75374cc3ull};
  1467. public:
  1468. //wyhash main function https://github.com/wangyi-fudan/wyhash
  1469. static uint64_t wyhashstr(const char *key, const size_t len)
  1470. {
  1471. uint64_t a = 0, b = 0, seed = secret[0];
  1472. const uint8_t *p = (const uint8_t*)key;
  1473. if (EMH_LIKELY(len <= 16)) {
  1474. if (EMH_LIKELY(len >= 4)) {
  1475. const auto half = (len >> 3) << 2;
  1476. a = (wyr4(p) << 32U) | wyr4(p + half); p += len - 4;
  1477. b = (wyr4(p) << 32U) | wyr4(p - half);
  1478. } else if (len) {
  1479. a = wyr3(p, len);
  1480. }
  1481. } else {
  1482. size_t i = len;
  1483. if (EMH_UNLIKELY(i > 48)) {
  1484. uint64_t see1 = seed, see2 = seed;
  1485. do {
  1486. seed = wymix(wyr8(p + 0) ^ secret[1], wyr8(p + 8) ^ seed);
  1487. see1 = wymix(wyr8(p + 16) ^ secret[2], wyr8(p + 24) ^ see1);
  1488. see2 = wymix(wyr8(p + 32) ^ secret[3], wyr8(p + 40) ^ see2);
  1489. p += 48; i -= 48;
  1490. } while (EMH_LIKELY(i > 48));
  1491. seed ^= see1 ^ see2;
  1492. }
  1493. while (i > 16) {
  1494. seed = wymix(wyr8(p) ^ secret[1], wyr8(p + 8) ^ seed);
  1495. i -= 16; p += 16;
  1496. }
  1497. a = wyr8(p + i - 16);
  1498. b = wyr8(p + i - 8);
  1499. }
  1500. return wymix(secret[1] ^ len, wymix(a ^ secret[1], b ^ seed));
  1501. }
  1502. #endif
  1503. private:
  1504. template<typename UType, typename std::enable_if<std::is_integral<UType>::value, uint32_t>::type = 0>
  1505. inline uint64_t hash_key(const UType key) const
  1506. {
  1507. #if EMH_INT_HASH
  1508. return hash64(key);
  1509. #elif EMH_IDENTITY_HASH
  1510. return key + (key >> 24);
  1511. #else
  1512. return _hasher(key);
  1513. #endif
  1514. }
  1515. template<typename UType, typename std::enable_if<std::is_same<UType, std::string>::value, uint32_t>::type = 0>
  1516. inline uint64_t hash_key(const UType& key) const
  1517. {
  1518. #if EMH_WYHASH_HASH
  1519. return wyhashstr(key.data(), key.size());
  1520. #else
  1521. return _hasher(key);
  1522. #endif
  1523. }
  1524. template<typename UType, typename std::enable_if<!std::is_integral<UType>::value && !std::is_same<UType, std::string>::value, uint32_t>::type = 0>
  1525. inline uint64_t hash_key(const UType& key) const
  1526. {
  1527. return _hasher(key);
  1528. }
  1529. private:
  1530. Index* _index;
  1531. value_type*_pairs;
  1532. HashT _hasher;
  1533. EqT _eq;
  1534. uint32_t _mlf;
  1535. size_type _mask;
  1536. size_type _num_buckets;
  1537. size_type _num_filled;
  1538. size_type _last;
  1539. #if EMH_HIGH_LOAD
  1540. size_type _ehead;
  1541. #endif
  1542. size_type _etail;
  1543. };
  1544. } // namespace emhash