source: mainline/uspace/lib/cpp/include/impl/algorithm.hpp@ 09416c12

lfn serial ticket/834-toolchain-update topic/msim-upgrade topic/simplify-dev-export
Last change on this file since 09416c12 was 09416c12, checked in by Dzejrou <dzejrou@…>, 7 years ago

cpp: added basic algorithm tests and fixed bugs they found

  • Property mode set to 100644
File size: 28.3 KB
RevLine 
[3457e11]1/*
2 * Copyright (c) 2017 Jaroslav Jindrak
3 * All rights reserved.
4 *
5 * Redistribution and use in source and binary forms, with or without
6 * modification, are permitted provided that the following conditions
7 * are met:
8 *
9 * - Redistributions of source code must retain the above copyright
10 * notice, this list of conditions and the following disclaimer.
11 * - Redistributions in binary form must reproduce the above copyright
12 * notice, this list of conditions and the following disclaimer in the
13 * documentation and/or other materials provided with the distribution.
14 * - The name of the author may not be used to endorse or promote products
15 * derived from this software without specific prior written permission.
16 *
17 * THIS SOFTWARE IS PROVIDED BY THE AUTHOR ``AS IS'' AND ANY EXPRESS OR
18 * IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES
19 * OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED.
20 * IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR ANY DIRECT, INDIRECT,
21 * INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT
22 * NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
23 * DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
24 * THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
25 * (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF
26 * THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
27 */
28
29#ifndef LIBCPP_ALGORITHM
30#define LIBCPP_ALGORITHM
31
[f041811]32#include <iterator>
[b6d68a3]33#include <utility>
34
[3457e11]35namespace std
36{
[f9823e2]37 template<class T>
38 struct less;
39
[b6d68a3]40 /**
41 * 25.2, non-modyfing sequence operations:
42 */
43
44 /**
45 * 25.2.1, all_of:
46 */
47
48 template<class InputIterator, class Predicate>
49 bool all_of(InputIterator first, InputIterator last, Predicate pred)
50 {
51 while (first != last)
52 {
53 if (!pred(*first++))
54 return false;
55 }
56
57 return true;
58 }
59
60 /**
61 * 25.2.2, any_of:
62 */
63
64 template<class InputIterator, class Predicate>
65 bool any_of(InputIterator first, InputIterator last, Predicate pred)
66 {
67 while (first != last)
68 {
69 if (pred(*first++))
70 return true;
71 }
72
73 return false;
74 }
75
76 /**
77 * 25.2.3, none_of:
78 */
79
80 template<class InputIterator, class Predicate>
81 bool none_of(InputIterator first, InputIterator last, Predicate pred)
82 {
83 return !any_of(first, last, pred);
84 }
85
86 /**
87 * 25.2.4, for_each:
88 */
89
90 template<class InputIterator, class Function>
91 Function for_each(InputIterator first, InputIterator last, Function f)
92 {
93 while (first != last)
94 f(*first++);
95
96 return move(f);
97 }
98
99 /**
100 * 25.2.5, find:
101 */
102
103 template<class InputIterator, class T>
104 InputIterator find(InputIterator first, InputIterator last, const T& value)
105 {
106 while (first != last)
107 {
108 if (*first == value)
109 return first;
110 ++first;
111 }
112
113 return last;
114 }
115
116 template<class InputIterator, class Predicate>
117 InputIterator find_if(InputIterator first, InputIterator last, Predicate pred)
118 {
119 while (first != last)
120 {
121 if (pred(*first))
122 return first;
123 ++first;
124 }
125
126 return last;
127 }
128
129 template<class InputIterator, class Predicate>
130 InputIterator find_if_not(InputIterator first, InputIterator last, Predicate pred)
131 {
132 while (first != last)
133 {
134 if (!pred(*first))
135 return first;
136 ++first;
137 }
138
139 return last;
140 }
141
142 /**
143 * 25.2.6, find_end:
144 */
145
146 // TODO: implement
147
148 /**
149 * 25.2.7, find_first:
150 */
151
152 // TODO: implement
153
154 /**
155 * 25.2.8, adjacent_find:
156 */
157
158 template<class ForwardIterator>
159 ForwardIterator adjacent_find(ForwardIterator first, ForwardIterator last)
160 {
161 while (first != last)
162 {
163 if (*first == *(first + 1))
164 return first;
165 ++first;
166 }
167
168 return last;
169 }
170
171 template<class ForwardIterator, class Predicate>
172 ForwardIterator adjacent_find(ForwardIterator first, ForwardIterator last, Predicate pred)
173 {
174 while (first != last)
175 {
176 if (pred(*first, *(first + 1)))
177 return first;
178 ++first;
179 }
180
181 return last;
182 }
183
184 /**
185 * 25.2.9, count:
186 */
187
188 template<class InputIterator, class T>
189 typename iterator_traits<InputIterator>::difference_type
190 count(InputIterator first, InputIterator last, const T& value)
191 {
192 typename iterator_traits<InputIterator>::difference_type cnt{};
193
194 while (first != last)
195 {
196 if (*first++ == value)
197 ++cnt;
198 }
199
200 return cnt;
201 }
202
203 template<class InputIterator, class Predicate>
204 typename iterator_traits<InputIterator>::difference_type
[09416c12]205 count_if(InputIterator first, InputIterator last, Predicate pred)
[b6d68a3]206 {
207 typename iterator_traits<InputIterator>::difference_type cnt{};
208
209 while (first != last)
210 {
211 if (pred(*first++))
212 ++cnt;
213 }
214
215 return cnt;
216 }
217
218 /**
219 * 25.2.10, mismatch:
220 */
221
222 template<class InputIterator1, class InputIterator2>
223 pair<InputIterator1, InputIterator2> mismatch(InputIterator1 first1, InputIterator1 last1,
224 InputIterator2 first2)
225 {
226 while (first1 != last1 && *first1++ == *first2++)
227 { /* DUMMY BODY */ }
228
229 return make_pair(first1, first2);
230 }
231
232 template<class InputIterator1, class InputIterator2, class BinaryPredicate>
233 pair<InputIterator1, InputIterator2> mismatch(InputIterator1 first1, InputIterator1 last1,
234 InputIterator2 first2, BinaryPredicate pred)
235 {
236 while (first1 != last1 && pred(*first1++, *first2++))
237 { /* DUMMY BODY */ }
238
239 return make_pair(first1, first2);
240 }
241
242 template<class InputIterator1, class InputIterator2>
243 pair<InputIterator1, InputIterator2> mismatch(InputIterator1 first1, InputIterator1 last1,
244 InputIterator2 first2, InputIterator2 last2)
245 {
246 while (first1 != last1 && first2 != last2 && *first1++ == *first2++)
247 { /* DUMMY BODY */ }
248
249 return make_pair(first1, first2);
250 }
251
252 template<class InputIterator1, class InputIterator2, class BinaryPredicate>
253 pair<InputIterator1, InputIterator2> mismatch(InputIterator1 first1, InputIterator1 last1,
254 InputIterator2 first2, InputIterator2 last2,
255 BinaryPredicate pred)
256 {
257 while (first1 != last1 && first2 != last2 && pred(*first1++, *first2++))
258 { /* DUMMY BODY */ }
259
260 return make_pair(first1, first2);
261 }
262
263 /**
264 * 25.2.11, equal:
265 */
266
267 template<class InputIterator1, class InputIterator2>
268 bool equal(InputIterator1 first1, InputIterator1 last1, InputIterator2 first2)
269 {
270 auto last2 = first2 + (last1 - first1);
271
272 return equal(first1, last1, first2, last2);
273 }
274
275 template<class InputIterator1, class InputIterator2>
276 bool equal(InputIterator1 first1, InputIterator1 last1,
277 InputIterator2 first2, InputIterator2 last2)
278 {
279 if ((last1 - first1) != (last2 - first2))
280 return false;
281
282 while (first1 != last1)
283 {
284 if (*first1++ != *first2++)
285 return false;
286 }
287
288 return true;
289 }
290
291 template<class InputIterator1, class InputIterator2, class BinaryPredicate>
292 bool equal(InputIterator1 first1, InputIterator1 last1,
293 InputIterator2 first2, BinaryPredicate pred)
294 {
295 auto last2 = first2 + (last1 - first1);
296
297 return equal(first1, last1, first2, last2, pred);
298 }
299
300 template<class InputIterator1, class InputIterator2, class BinaryPredicate>
301 bool equal(InputIterator1 first1, InputIterator1 last1,
302 InputIterator2 first2, InputIterator2 last2,
303 BinaryPredicate pred)
304 {
305 if ((last1 - first1) != (last2 - first2))
306 return false;
307
308 while (first1 != last1)
309 {
310 if (!pred(*first1++, *first2++))
311 return false;
312 }
313
314 return true;
315 }
316
317 /**
318 * 25.2.12, is_permutation:
319 */
320
321 // TODO: implement
322
323 /**
324 * 25.2.13, search:
325 */
326
327 // TODO: implement
328
329 /**
330 * 25.3, mutating sequence operations:
331 */
332
[3457e11]333 /**
334 * 25.3.1, copy:
335 */
336
337 template<class InputIterator, class OutputIterator>
338 OutputIterator copy(InputIterator first, InputIterator last, OutputIterator result)
339 {
[b6d68a3]340 while (first != last)
[4f202b20]341 *result++ = *first++;
[3457e11]342
343 return result;
344 }
345
[b6d68a3]346 template<class InputIterator, class Size, class OutputIterator>
347 OutputIterator copy_n(InputIterator first, Size count, OutputIterator result)
348 {
349 for (Size i = 0; i < count; ++i, ++first, ++result)
350 *result = *first;
351
352 return result;
353 }
354
355 template<class InputIterator, class OutputIterator, class Predicate>
356 OutputIterator copy_if(InputIterator first, InputIterator last,
357 OutputIterator result, Predicate pred)
358 {
359 while (first != last)
360 {
361 if (pred(*first))
362 *result++ = *first;
363 ++first;
364 }
365
366 return result;
367 }
368
369 template<class BidirectionalIterator1, class BidirectionalIterator2>
370 BidirectionalIterator2 copy_backward(BidirectionalIterator1 first, BidirectionalIterator1 last,
371 BidirectionalIterator2 result)
372 {
373 // Note: We're copying [first, last) so we need to skip the initial value of last.
374 while (last-- != first)
375 *result-- = *last;
[4f202b20]376
377 return result;
[b6d68a3]378 }
379
[3457e11]380 /**
[b6d68a3]381 * 25.3.2, move:
[3457e11]382 */
383
[b6d68a3]384 template<class InputIterator, class OutputIterator>
385 OutputIterator move(InputIterator first, InputIterator last, OutputIterator result)
[3457e11]386 {
[b6d68a3]387 while (first != last)
388 *result++ = move(first++);
389
390 return result;
391 }
392
393 template<class BidirectionalIterator1, class BidirectionalIterator2>
394 BidirectionalIterator2 move_backward(BidirectionalIterator1 first, BidirectionalIterator1 last,
395 BidirectionalIterator2 result)
396 {
397 // Note: We're copying [first, last) so we need to skip the initial value of last.
398 while (last-- != first)
399 *result-- = move(*last);
400 }
401
402 /**
403 * 25.3.3, swap:
404 */
405
406 template<class ForwardIterator1, class ForwardIterator2>
407 ForwardIterator2 swap_ranges(ForwardIterator1 first1, ForwardIterator1 last1,
408 ForwardIterator2 first2)
409 {
410 while (first1 != last1)
411 swap(*first1++, *first2++);
412
413 return first2;
414 }
415
416 template<class ForwardIterator1, class ForwardIterator2>
417 void iter_swap(ForwardIterator1 iter1, ForwardIterator2 iter2)
418 {
419 swap(*iter1, *iter2);
420 }
421
422 /**
423 * 25.3.4, transform:
424 */
425
426 template<class InputIterator, class OutputIterator, class UnaryOperation>
427 OutputIterator transform(InputIterator first, InputIterator last,
428 OutputIterator result, UnaryOperation op)
429 {
430 while (first != last)
431 *result++ = op(*first++);
432 }
433
434 template<class InputIterator1, class InputIterator2,
435 class OutputIterator, class BinaryOperation>
436 OutputIterator transform(InputIterator1 first1, InputIterator1 last1,
437 InputIterator2 first2, OutputIterator result,
438 BinaryOperation op)
439 {
440 while (first1 != last1)
441 *result++ = op(*first1++, *first2++);
442 }
443
444 /**
445 * 25.3.5, replace:
446 */
447
448 template<class ForwardIterator, class T>
449 void replace(ForwardIterator first, ForwardIterator last,
450 const T& old_value, const T& new_value)
451 {
452 while (first != last)
453 {
454 if (*first == old_value)
455 *first = new_value;
456 ++first;
457 }
458 }
459
460 template<class ForwardIterator, class Predicate, class T>
461 void replace_if(ForwardIterator first, ForwardIterator last,
462 Predicate pred, const T& new_value)
463 {
464 while (first != last)
465 {
466 if (pred(*first))
467 *first = new_value;
468 ++first;
469 }
470 }
471
472 template<class InputIterator, class OutputIterator, class T>
473 OutputIterator replace_copy(InputIterator first, InputIterator last,
474 OutputIterator result, const T& old_value,
475 const T& new_value)
476 {
477 while (first != last)
478 {
479 if (*first == old_value)
480 *result = new_value;
481 else
482 *result = *first;
483
484 ++first;
485 ++result;
486 }
487 }
488
489 template<class InputIterator, class OutputIterator, class Predicate, class T>
490 OutputIterator replace_copy_if(InputIterator first, InputIterator last,
491 OutputIterator result, Predicate pred,
492 const T& new_value)
493 {
494 while (first != last)
495 {
496 if (pred(*first))
497 *result = new_value;
498 else
499 *result = *first;
500
501 ++first;
502 ++result;
503 }
[3457e11]504 }
505
[b6d68a3]506 /**
507 * 25.3.6, fill:
508 */
509
510 template<class ForwardIterator, class T>
511 void fill(ForwardIterator first, ForwardIterator last, const T& value)
512 {
513 while (first != last)
514 *first++ = value;
515 }
516
517 template<class InputIterator, class Size, class T>
518 void fill_n(InputIterator first, Size count, const T& value)
519 {
520 for (Size i = 0; i < count; ++i)
521 *first++ = value;
522 }
523
524 /**
525 * 25.3.7, generate:
526 */
527
528 template<class ForwardIterator, class Generator>
529 void generate(ForwardIterator first, ForwardIterator last,
530 Generator gen)
531 {
532 while (first != last)
533 *first++ = gen();
534 }
535
536 template<class OutputIterator, class Size, class Generator>
537 void generate(OutputIterator first, Size count, Generator gen)
538 {
539 for (Size i = 0; i < count; ++i)
540 *first++ = gen();
541 }
542
543 /**
544 * 25.3.8, remove:
545 */
546
547 template<class ForwardIterator, class T>
548 ForwardIterator remove(ForwardIterator first, ForwardIterator last,
549 const T& value)
550 {
551 auto it = first;
552 while (it != last)
553 {
554 if (*it != value)
555 *first++ = move(*it);
556 }
557
558 return first;
559 }
560
561 template<class ForwardIterator, class Predicate>
562 ForwardIterator remove_if(ForwardIterator first, ForwardIterator last,
563 Predicate pred)
564 {
565 auto it = first;
566 while (it != last)
567 {
568 if (!pred(*it))
569 *first++ = move(*it);
570 }
571
572 return first;
573 }
574
575 template<class InputIterator, class OutputIterator, class T>
576 OutputIterator remove_copy(InputIterator first, InputIterator last,
577 OutputIterator result, const T& value)
578 {
579 while (first != last)
580 {
581 if (*first != value)
582 *result++ = *first;
583 ++first;
584 }
585
586 return result;
587 }
588
589 template<class InputIterator, class OutputIterator, class Predicate>
590 OutputIterator remove_copy_if(InputIterator first, InputIterator last,
591 OutputIterator result, Predicate pred)
592 {
593 while (first != last)
594 {
595 if (!pred(*first))
596 *result++ = *first;
597 ++first;
598 }
599
600 return result;
601 }
602
603 /**
604 * 25.3.9, unique:
605 */
606
607 // TODO: implement
608
609 /**
610 * 25.3.10, reverse:
611 */
612
613 template<class BidirectionalIterator>
614 void reverse(BidirectionalIterator first, BidirectionalIterator last)
615 {
616 if (first == last)
617 return;
618 auto mid_count = (last - first) / 2;
619
620 --last;
621 for (decltype(mid_count) i = 0; i < mid_count; ++i)
622 iter_swap(first++, last--);
623 }
624
625 template<class BidirectionalIterator, class OutputIterator>
626 OutputIterator reverse_copy(BidirectionalIterator first,
627 BidirectionalIterator last,
628 OutputIterator result)
629 {
630 while (--last != first)
631 *result++ = *last;
632 }
633
634 /**
635 * 25.3.11, rotate:
636 */
637
638 // TODO: implement
639
640 /**
641 * 25.3.12, shuffle:
642 */
643
644 // TODO: implement
645
646 /**
647 * 25.3.13, partitions:
648 */
649
650 // TODO: implement
651
652 /**
653 * 25.4, sorting and related operations:
654 */
655
656 /**
657 * 25.4.1, sorting:
658 */
659
660 /**
661 * 25.4.1.1, sort:
662 */
663
[7ea90cf]664 template<class RandomAccessIterator, class Compare>
665 void make_heap(RandomAccessIterator, RandomAccessIterator,
666 Compare);
667
668 template<class RandomAccessIterator, class Compare>
669 void sort_heap(RandomAccessIterator, RandomAccessIterator,
670 Compare);
671
672 template<class RandomAccessIterator>
673 void sort(RandomAccessIterator first, RandomAccessIterator last)
674 {
675 using value_type = typename iterator_traits<RandomAccessIterator>::value_type;
676
677 sort(first, last, less<value_type>{});
678 }
679
680 template<class RandomAccessIterator, class Compare>
681 void sort(RandomAccessIterator first, RandomAccessIterator last,
682 Compare comp)
683 {
684 /**
685 * Note: This isn't the most effective approach,
686 * but since we already have these two functions
687 * and they satisfy asymptotic limitations
688 * imposed by the standard, we're using them at
689 * the moment. Might be good to change it to qsort
690 * or merge sort later.
691 */
692
693 make_heap(first, last, comp);
694 sort_heap(first, last, comp);
695 }
[b6d68a3]696
697 /**
698 * 25.4.1.2, stable_sort:
699 */
700
701 // TODO: implement
702
703 /**
704 * 25.4.1.3, partial_sort:
705 */
706
707 // TODO: implement
708
709 /**
710 * 25.4.1.4, partial_sort_copy:
711 */
712
713 // TODO: implement
714
715 /**
716 * 25.4.1.5, is_sorted:
717 */
718
719 template<class ForwardIterator>
720 bool is_sorted(ForwardIterator first, ForwardIterator last)
721 {
722 return is_sorted_until(first, last) == last;
723 }
724
725 template<class ForwardIterator, class Comp>
726 bool is_sorted(ForwardIterator first, ForwardIterator last,
727 Comp comp)
728 {
729 return is_sorted_until(first, last, comp) == last;
730 }
731
732 template<class ForwardIterator>
733 ForwardIterator is_sorted_until(ForwardIterator first, ForwardIterator last)
734 {
735 if (distance(first, last) < 2)
736 return last;
737
738 while (first != last)
739 {
740 if (*first > *(++first))
741 return first;
742 }
743
744 return last;
745 }
746
747 template<class ForwardIterator, class Comp>
748 ForwardIterator is_sorted_until(ForwardIterator first, ForwardIterator last,
749 Comp comp)
750 {
751 if (distance(first, last) < 2)
752 return last;
753
754 while (first != last)
755 {
756 if (!comp(*first, *(++first)))
757 return first;
758 }
759
760 return last;
761 }
762
763 /**
764 * 25.4.2, nth_element:
765 */
766
767 // TODO: implement
768
769 /**
770 * 25.4.3, binary search:
771 */
772
773 /**
774 * 25.4.3.1, lower_bound
775 */
776
777 // TODO: implement
778
779 /**
780 * 25.4.3.2, upper_bound
781 */
782
783 // TODO: implement
784
785 /**
786 * 25.4.3.3, equal_range:
787 */
788
789 // TODO: implement
790
791 /**
792 * 25.4.3.4, binary_search:
793 */
794
795 // TODO: implement
796
797 /**
798 * 25.4.4, merge:
799 */
800
801 // TODO: implement
802
803 /**
804 * 25.4.5, set operations on sorted structures:
805 */
806
807 /**
808 * 25.4.5.1, includes:
809 */
810
811 // TODO: implement
812
813 /**
814 * 25.4.5.2, set_union:
815 */
816
817 // TODO: implement
818
819 /**
820 * 25.4.5.3, set_intersection:
821 */
822
823 // TODO: implement
824
825 /**
826 * 25.4.5.4, set_difference:
827 */
828
829 // TODO: implement
830
831 /**
832 * 25.4.5.5, set_symmetric_difference:
833 */
834
835 // TODO: implement
836
837 /**
838 * 25.4.6, heap operations:
839 */
840
[f9823e2]841 namespace aux
842 {
843 template<class T>
844 T heap_parent(T idx)
845 {
846 return (idx - 1) / 2;
847 }
848
849 template<class T>
850 T heap_left_child(T idx)
851 {
852 return 2 * idx + 1;
853 }
854
855 template<class T>
856 T heap_right_child(T idx)
857 {
858 return 2 * idx + 2;
859 }
860
861 template<class RandomAccessIterator, class Size, class Compare>
862 void correct_children(RandomAccessIterator first,
863 Size idx, Size count, Compare comp)
864 {
865 using aux::heap_left_child;
866 using aux::heap_right_child;
867
868 auto left = heap_left_child(idx);
869 auto right = heap_right_child(idx);
870
871 bool left_incorrect{comp(first[idx], first[left])};
872 bool right_incorrect{comp(first[idx], first[right])};
873 while ((left < count && left_incorrect) ||
874 (right < count && right_incorrect))
875 {
876 if (right >= count || (left_incorrect && comp(first[right], first[left])))
877 {
878 swap(first[idx], first[left]);
879
880 idx = left;
881 }
882 else if (right < count && right_incorrect)
883 {
884 swap(first[idx], first[right]);
885
886 idx = right;
887 } // Else should not happen because of the while condition.
888
889 left = heap_left_child(idx);
890 right = heap_right_child(idx);
891
892 left_incorrect = comp(first[idx], first[left]);
893 right_incorrect = comp(first[idx], first[right]);
894 }
895 }
896 }
897
[b6d68a3]898 /**
899 * 25.4.6.1, push_heap:
900 */
901
[f9823e2]902 template<class RandomAccessIterator>
903 void push_heap(RandomAccessIterator first,
904 RandomAccessIterator last)
905 {
906 using value_type = typename iterator_traits<RandomAccessIterator>::value_type;
907
908 push_heap(first, last, less<value_type>{});
909 }
910
911 template<class RandomAccessIterator, class Compare>
912 void push_heap(RandomAccessIterator first,
913 RandomAccessIterator last,
914 Compare comp)
915 {
916 using aux::heap_parent;
917
918 auto count = distance(first, last);
919 if (count <= 1)
920 return;
921
922 auto idx = count - 1;
923 auto parent = heap_parent(idx);
924 while (idx > 0 && comp(first[parent], first[idx]))
925 {
926 swap(first[idx], first[parent]);
927
928 idx = parent;
929 parent = heap_parent(idx);
930 }
931 }
[b6d68a3]932
933 /**
934 * 25.4.6.2, pop_heap:
935 */
936
[f9823e2]937 template<class RandomAccessIterator>
938 void pop_heap(RandomAccessIterator first,
939 RandomAccessIterator last)
940 {
941 using value_type = typename iterator_traits<RandomAccessIterator>::value_type;
942
943 pop_heap(first, last, less<value_type>{});
944 }
945
946 template<class RandomAccessIterator, class Compare>
947 void pop_heap(RandomAccessIterator first,
948 RandomAccessIterator last,
949 Compare comp)
950 {
951 auto count = distance(first, last);
952 if (count <= 1)
953 return;
954
955 swap(first[0], first[count - 1]);
956 aux::correct_children(first, decltype(count){}, count - 2, comp);
957 }
[b6d68a3]958
959 /**
960 * 25.4.6.3, make_heap:
961 */
962
[f9823e2]963 template<class RandomAccessIterator>
964 void make_heap(RandomAccessIterator first,
965 RandomAccessIterator last)
966 {
967 using value_type = typename iterator_traits<RandomAccessIterator>::value_type;
968
969 make_heap(first, last, less<value_type>{});
970 }
971
972 template<class RandomAccessIterator, class Compare>
973 void make_heap(RandomAccessIterator first,
974 RandomAccessIterator last,
975 Compare comp)
976 {
977 auto count = distance(first, last);
978 if (count <= 1)
979 return;
980
981 for (auto i = count; i > 0; --i)
982 {
983 auto idx = i - 1;
984
985 aux::correct_children(first, idx, count, comp);
986 }
987 }
[b6d68a3]988
989 /**
990 * 25.4.6.4, sort_heap:
991 */
992
[f9823e2]993 template<class RandomAccessIterator>
994 void sort_heap(RandomAccessIterator first,
995 RandomAccessIterator last)
996 {
997 using value_type = typename iterator_traits<RandomAccessIterator>::value_type;
998
999 sort_heap(first, last, less<value_type>{});
1000 }
1001
1002 template<class RandomAccessIterator, class Compare>
1003 void sort_heap(RandomAccessIterator first,
1004 RandomAccessIterator last,
1005 Compare comp)
1006 {
1007 while (first != last)
1008 pop_heap(first, last--, comp);
1009 }
[b6d68a3]1010
1011 /**
1012 * 25.4.6.5, is_heap:
1013 */
1014
[f9823e2]1015 template<class RandomAccessIterator>
1016 auto is_heap_until(RandomAccessIterator first, RandomAccessIterator last)
1017 {
1018 using value_type = typename iterator_traits<RandomAccessIterator>::value_type;
1019
1020 return is_heap_until(first, last, less<value_type>{});
1021 }
1022
1023 template<class RandomAccessIterator, class Compare>
1024 auto is_heap_until(RandomAccessIterator first, RandomAccessIterator last,
1025 Compare comp)
1026 {
1027 using aux::heap_left_child;
1028 using aux::heap_right_child;
1029
1030 auto count = distance(first, last);
1031 if (count < 2)
1032 return last;
1033
1034 auto res = first;
1035 for (decltype(count) idx = 0; idx < count; ++idx)
1036 {
1037 auto left = heap_left_child(idx);
1038 auto right = heap_right_child(idx);
1039
1040 if (left < count && comp(first[idx], first[left]))
1041 return res;
1042 if (right < count && comp(first[idx], first[right]))
1043 return res;
1044
1045 ++res;
1046 }
1047
1048 return res;
1049 }
1050
1051 template<class RandomAccessIterator>
1052 bool is_heap(RandomAccessIterator first, RandomAccessIterator last)
1053 {
1054 return is_heap_until(first, last) == last;
1055 }
1056
1057 template<class RandomAccessIterator, class Compare>
1058 bool is_heap(RandomAccessIterator first, RandomAccessIterator last,
1059 Compare comp)
1060 {
1061 return is_heap_until(first, last, comp) == last;
1062 }
[b6d68a3]1063
1064 /**
1065 * 25.4.7, minimum and maximum:
1066 * // TODO: implement container versions when we have
1067 * numeric limits and min/max element
1068 * // TODO: versions with comparators
1069 * // TODO: minmax
1070 */
1071
[3457e11]1072 template<class T>
1073 constexpr const T& min(const T& lhs, const T& rhs)
1074 {
1075 return (lhs < rhs) ? lhs : rhs;
1076 }
[b6d68a3]1077
1078 template<class T>
1079 constexpr const T& max(const T& lhs, const T& rhs)
1080 {
1081 return (lhs > rhs) ? lhs : rhs;
1082 }
1083
1084 /**
1085 * 25.4.8, lexicographical comparison:
1086 */
1087
[26d2990]1088 template<class InputIterator1, class InputIterator2>
1089 bool lexicographical_compare(InputIterator1 first1,
1090 InputIterator1 last1,
1091 InputIterator2 first2,
1092 InputIterator2 last2)
1093 {
1094 /**
1095 * *first1 and *first2 can have different types
1096 * so we use a transparent comparator.
1097 */
1098 return lexicographical_compare(
1099 first1, last1, first2, last2,
1100 less<void>{}
1101 );
1102 }
1103
1104 template<class InputIterator1, class InputIterator2, class Compare>
1105 bool lexicographical_compare(InputIterator1 first1,
1106 InputIterator1 last1,
1107 InputIterator2 first2,
1108 InputIterator2 last2,
1109 Compare comp)
1110 {
1111 while ((first1 != last1) && (first2 != last2))
1112 {
1113 if (comp(*first1, *first2))
1114 return true;
1115 if (comp(*first2, *first1))
1116 return false;
1117
1118 ++first1;
1119 ++first2;
1120 }
1121
1122 /**
1123 * Up until now they are same, so we have to check
1124 * if we reached the end on one.
1125 */
1126 return (first1 == last1) && (first2 != last2);
1127 }
[b6d68a3]1128
1129 /**
1130 * 25.4.9, permutation generators:
1131 */
1132
1133 // TODO: implement
[3457e11]1134}
1135
1136#endif
Note: See TracBrowser for help on using the repository browser.