... | ... |
@@ -614,257 +614,257 @@ |
614 | 614 |
|
615 | 615 |
private: |
616 | 616 |
const UnionFindEnum* unionFind; |
617 | 617 |
int idx, fdx; |
618 | 618 |
}; |
619 | 619 |
|
620 | 620 |
}; |
621 | 621 |
|
622 | 622 |
/// \ingroup auxdat |
623 | 623 |
/// |
624 | 624 |
/// \brief A \e Extend-Find data structure implementation which |
625 | 625 |
/// is able to enumerate the components. |
626 | 626 |
/// |
627 | 627 |
/// The class implements an \e Extend-Find data structure which is |
628 | 628 |
/// able to enumerate the components and the items in a |
629 | 629 |
/// component. The data structure is a simplification of the |
630 | 630 |
/// Union-Find structure, and it does not allow to merge two components. |
631 | 631 |
/// |
632 | 632 |
/// \pre You need to add all the elements by the \ref insert() |
633 | 633 |
/// method. |
634 | 634 |
template <typename IM> |
635 | 635 |
class ExtendFindEnum { |
636 | 636 |
public: |
637 | 637 |
|
638 | 638 |
///\e |
639 | 639 |
typedef IM ItemIntMap; |
640 | 640 |
///\e |
641 | 641 |
typedef typename ItemIntMap::Key Item; |
642 | 642 |
|
643 | 643 |
private: |
644 | 644 |
|
645 | 645 |
ItemIntMap& index; |
646 | 646 |
|
647 | 647 |
struct ItemT { |
648 | 648 |
int cls; |
649 | 649 |
Item item; |
650 | 650 |
int next, prev; |
651 | 651 |
}; |
652 | 652 |
|
653 | 653 |
std::vector<ItemT> items; |
654 | 654 |
int firstFreeItem; |
655 | 655 |
|
656 | 656 |
struct ClassT { |
657 | 657 |
int firstItem; |
658 | 658 |
int next, prev; |
659 | 659 |
}; |
660 | 660 |
|
661 | 661 |
std::vector<ClassT> classes; |
662 | 662 |
|
663 | 663 |
int firstClass, firstFreeClass; |
664 | 664 |
|
665 | 665 |
int newClass() { |
666 | 666 |
if (firstFreeClass != -1) { |
667 | 667 |
int cdx = firstFreeClass; |
668 | 668 |
firstFreeClass = classes[cdx].next; |
669 | 669 |
return cdx; |
670 | 670 |
} else { |
671 | 671 |
classes.push_back(ClassT()); |
672 | 672 |
return classes.size() - 1; |
673 | 673 |
} |
674 | 674 |
} |
675 | 675 |
|
676 | 676 |
int newItem() { |
677 | 677 |
if (firstFreeItem != -1) { |
678 | 678 |
int idx = firstFreeItem; |
679 | 679 |
firstFreeItem = items[idx].next; |
680 | 680 |
return idx; |
681 | 681 |
} else { |
682 | 682 |
items.push_back(ItemT()); |
683 | 683 |
return items.size() - 1; |
684 | 684 |
} |
685 | 685 |
} |
686 | 686 |
|
687 | 687 |
public: |
688 | 688 |
|
689 | 689 |
/// \brief Constructor |
690 | 690 |
ExtendFindEnum(ItemIntMap& _index) |
691 | 691 |
: index(_index), items(), firstFreeItem(-1), |
692 | 692 |
classes(), firstClass(-1), firstFreeClass(-1) {} |
693 | 693 |
|
694 | 694 |
/// \brief Inserts the given element into a new component. |
695 | 695 |
/// |
696 | 696 |
/// This method creates a new component consisting only of the |
697 | 697 |
/// given element. |
698 | 698 |
int insert(const Item& item) { |
699 | 699 |
int cdx = newClass(); |
700 | 700 |
classes[cdx].prev = -1; |
701 | 701 |
classes[cdx].next = firstClass; |
702 | 702 |
if (firstClass != -1) { |
703 | 703 |
classes[firstClass].prev = cdx; |
704 | 704 |
} |
705 | 705 |
firstClass = cdx; |
706 | 706 |
|
707 | 707 |
int idx = newItem(); |
708 | 708 |
items[idx].item = item; |
709 | 709 |
items[idx].cls = cdx; |
710 | 710 |
items[idx].prev = idx; |
711 | 711 |
items[idx].next = idx; |
712 | 712 |
|
713 | 713 |
classes[cdx].firstItem = idx; |
714 | 714 |
|
715 | 715 |
index.set(item, idx); |
716 | 716 |
|
717 | 717 |
return cdx; |
718 | 718 |
} |
719 | 719 |
|
720 | 720 |
/// \brief Inserts the given element into the given component. |
721 | 721 |
/// |
722 | 722 |
/// This methods inserts the element \e item a into the \e cls class. |
723 | 723 |
void insert(const Item& item, int cls) { |
724 | 724 |
int idx = newItem(); |
725 | 725 |
int rdx = classes[cls].firstItem; |
726 | 726 |
items[idx].item = item; |
727 | 727 |
items[idx].cls = cls; |
728 | 728 |
|
729 | 729 |
items[idx].prev = rdx; |
730 | 730 |
items[idx].next = items[rdx].next; |
731 | 731 |
items[items[rdx].next].prev = idx; |
732 | 732 |
items[rdx].next = idx; |
733 | 733 |
|
734 | 734 |
index.set(item, idx); |
735 | 735 |
} |
736 | 736 |
|
737 | 737 |
/// \brief Clears the union-find data structure |
738 | 738 |
/// |
739 | 739 |
/// Erase each item from the data structure. |
740 | 740 |
void clear() { |
741 | 741 |
items.clear(); |
742 |
classes.clear; |
|
742 |
classes.clear(); |
|
743 | 743 |
firstClass = firstFreeClass = firstFreeItem = -1; |
744 | 744 |
} |
745 | 745 |
|
746 | 746 |
/// \brief Gives back the class of the \e item. |
747 | 747 |
/// |
748 | 748 |
/// Gives back the class of the \e item. |
749 | 749 |
int find(const Item &item) const { |
750 | 750 |
return items[index[item]].cls; |
751 | 751 |
} |
752 | 752 |
|
753 | 753 |
/// \brief Gives back a representant item of the component. |
754 | 754 |
/// |
755 | 755 |
/// Gives back a representant item of the component. |
756 | 756 |
Item item(int cls) const { |
757 | 757 |
return items[classes[cls].firstItem].item; |
758 | 758 |
} |
759 | 759 |
|
760 | 760 |
/// \brief Removes the given element from the structure. |
761 | 761 |
/// |
762 | 762 |
/// Removes the element from its component and if the component becomes |
763 | 763 |
/// empty then removes that component from the component list. |
764 | 764 |
/// |
765 | 765 |
/// \warning It is an error to remove an element which is not in |
766 | 766 |
/// the structure. |
767 | 767 |
void erase(const Item &item) { |
768 | 768 |
int idx = index[item]; |
769 | 769 |
int cdx = items[idx].cls; |
770 | 770 |
|
771 | 771 |
if (idx == items[idx].next) { |
772 | 772 |
if (classes[cdx].prev != -1) { |
773 | 773 |
classes[classes[cdx].prev].next = classes[cdx].next; |
774 | 774 |
} else { |
775 | 775 |
firstClass = classes[cdx].next; |
776 | 776 |
} |
777 | 777 |
if (classes[cdx].next != -1) { |
778 | 778 |
classes[classes[cdx].next].prev = classes[cdx].prev; |
779 | 779 |
} |
780 | 780 |
classes[cdx].next = firstFreeClass; |
781 | 781 |
firstFreeClass = cdx; |
782 | 782 |
} else { |
783 | 783 |
classes[cdx].firstItem = items[idx].next; |
784 | 784 |
items[items[idx].next].prev = items[idx].prev; |
785 | 785 |
items[items[idx].prev].next = items[idx].next; |
786 | 786 |
} |
787 | 787 |
items[idx].next = firstFreeItem; |
788 | 788 |
firstFreeItem = idx; |
789 | 789 |
|
790 | 790 |
} |
791 | 791 |
|
792 | 792 |
|
793 | 793 |
/// \brief Removes the component of the given element from the structure. |
794 | 794 |
/// |
795 | 795 |
/// Removes the component of the given element from the structure. |
796 | 796 |
/// |
797 | 797 |
/// \warning It is an error to give an element which is not in the |
798 | 798 |
/// structure. |
799 | 799 |
void eraseClass(int cdx) { |
800 | 800 |
int idx = classes[cdx].firstItem; |
801 | 801 |
items[items[idx].prev].next = firstFreeItem; |
802 | 802 |
firstFreeItem = idx; |
803 | 803 |
|
804 | 804 |
if (classes[cdx].prev != -1) { |
805 | 805 |
classes[classes[cdx].prev].next = classes[cdx].next; |
806 | 806 |
} else { |
807 | 807 |
firstClass = classes[cdx].next; |
808 | 808 |
} |
809 | 809 |
if (classes[cdx].next != -1) { |
810 | 810 |
classes[classes[cdx].next].prev = classes[cdx].prev; |
811 | 811 |
} |
812 | 812 |
classes[cdx].next = firstFreeClass; |
813 | 813 |
firstFreeClass = cdx; |
814 | 814 |
} |
815 | 815 |
|
816 | 816 |
/// \brief LEMON style iterator for the classes. |
817 | 817 |
/// |
818 | 818 |
/// ClassIt is a lemon style iterator for the components. It iterates |
819 | 819 |
/// on the ids of classes. |
820 | 820 |
class ClassIt { |
821 | 821 |
public: |
822 | 822 |
/// \brief Constructor of the iterator |
823 | 823 |
/// |
824 | 824 |
/// Constructor of the iterator |
825 | 825 |
ClassIt(const ExtendFindEnum& ufe) : extendFind(&ufe) { |
826 | 826 |
cdx = extendFind->firstClass; |
827 | 827 |
} |
828 | 828 |
|
829 | 829 |
/// \brief Constructor to get invalid iterator |
830 | 830 |
/// |
831 | 831 |
/// Constructor to get invalid iterator |
832 | 832 |
ClassIt(Invalid) : extendFind(0), cdx(-1) {} |
833 | 833 |
|
834 | 834 |
/// \brief Increment operator |
835 | 835 |
/// |
836 | 836 |
/// It steps to the next representant item. |
837 | 837 |
ClassIt& operator++() { |
838 | 838 |
cdx = extendFind->classes[cdx].next; |
839 | 839 |
return *this; |
840 | 840 |
} |
841 | 841 |
|
842 | 842 |
/// \brief Conversion operator |
843 | 843 |
/// |
844 | 844 |
/// It converts the iterator to the current class id. |
845 | 845 |
operator int() const { |
846 | 846 |
return cdx; |
847 | 847 |
} |
848 | 848 |
|
849 | 849 |
/// \brief Equality operator |
850 | 850 |
/// |
851 | 851 |
/// Equality operator |
852 | 852 |
bool operator==(const ClassIt& i) { |
853 | 853 |
return i.cdx == cdx; |
854 | 854 |
} |
855 | 855 |
|
856 | 856 |
/// \brief Inequality operator |
857 | 857 |
/// |
858 | 858 |
/// Inequality operator |
859 | 859 |
bool operator!=(const ClassIt& i) { |
860 | 860 |
return i.cdx != cdx; |
861 | 861 |
} |
862 | 862 |
|
863 | 863 |
private: |
864 | 864 |
const ExtendFindEnum* extendFind; |
865 | 865 |
int cdx; |
866 | 866 |
}; |
867 | 867 |
|
868 | 868 |
/// \brief LEMON style iterator for the items of a component. |
869 | 869 |
/// |
870 | 870 |
/// ClassIt is a lemon style iterator for the components. It iterates |
0 comments (0 inline)