Google OR-Tools: ortools/constraint_solver/pack.cc Source File
24#include "absl/strings/str_format.h"
25#include "absl/strings/str_join.h"
26#include "absl/types/span.h"
43 virtual void InitialPropagate(int bin_index, const std::vector<int>& forced,
44 const std::vector<int>& undecided) = 0;
46 const std::vector<int>& assigned, const std::vector<int>& unassigned) = 0;
48 virtual void Propagate(int bin_index, const std::vector<int>& forced,
49 const std::vector<int>& removed) = 0;
51 const std::vector<int>& unassigned) = 0;
53 std::string DebugString() const override { return "Dimension"; }
58 bool IsUndecided(int var_index, int bin_index) const {
62 bool IsPossible(int var_index, int bin_index) const {
82 void SetAssigned(int var_index) { pack_->SetAssigned(var_index); }
84 void SetUnassigned(int var_index) { pack_->SetUnassigned(var_index); }
103 Solver* const solver_;
104 Pack* const pack_;
114 unprocessed_(new RevBitMatrix(bins_ + 1, vars_.size())),
121 for (int i = 0; i < vars_.size(); ++i) {
129 for (int i = 0; i < vars_.size(); ++i) {
130 IntVar* const var = vars_[i];
131 if (!var->Bound()) {
145 for (int bin_index = 0; bin_index <= bins_; ++bin_index) {
146 forced_[bin_index].clear();
156 for (int i = 0; i < to_set_.size(); ++i) {
157 vars_[to_set_[i].first]->SetValue(to_set_[i].second);
159 for (int i = 0; i < to_unset_.size(); ++i) {
160 vars_[to_unset_[i].first]->RemoveValue(to_unset_[i].second);
167class InitialPropagateData : public BaseObject {
169 explicit InitialPropagateData(size_t num_bins)
170 : BaseObject(), undecided_(num_bins) {}
171 void PushAssigned(int index) { assigned_.push_back(index); }
172 void PushUnassigned(int index) { unassigned_.push_back(index); }
173 void PushUndecided(int bin, int index) {
174 undecided_.at(bin).push_back(index);
176 const std::vector<int>& undecided(int bin) const {
177 return undecided_.at(bin);
179 const std::vector<int>& assigned() const { return assigned_; }
180 const std::vector<int>& unassigned() const { return unassigned_; }
182 std::string DebugString() const override { return "InitialPropagateData"; }
185 std::vector<std::vector<int>> undecided_;
186 std::vector<int> unassigned_;
187 std::vector<int> assigned_;
196 InitialPropagateData* data = s->RevAlloc(new InitialPropagateData(bins_));
197 for (int var_index = 0; var_index < vars_.size(); ++var_index) {
198 IntVar* const var = vars_[var_index];
199 var->SetRange(0, bins_);
200 if (var->Bound()) {
201 const int64_t value = var->Min();
203 forced_[value].push_back(var_index);
204 data->PushAssigned(var_index);
206 data->PushUnassigned(var_index);
209 DCHECK_GT(bins_, var->Min())
210 << "The variable maximum should be at most " << bins_
211 << ", and it should be unbound, so its min is expected to be less "
213 if (var->Max() < bins_) {
214 data->PushAssigned(var_index);
218 if (value >= 0 && value <= bins_) {
219 unprocessed_->SetToOne(s, value, var_index);
221 data->PushUndecided(value, var_index);
227 for (int bin_index = 0; bin_index < bins_; ++bin_index) {
230 absl::StrFormat("Pack(bin %d, forced = [%s], undecided = [%s])",
231 bin_index, absl::StrJoin(forced_[bin_index], ", "),
232 absl::StrJoin(data->undecided(bin_index), ", ")));
235 for (int dim_index = 0; dim_index < dims_.size(); ++dim_index) {
238 "InitialProgateDimension(%s)", dims_[dim_index]->DebugString()));
240 dims_[dim_index]->InitialPropagate(bin_index, forced_[bin_index],
241 data->undecided(bin_index));
252 absl::StrFormat("Pack(assigned = [%s], unassigned = [%s])",
253 absl::StrJoin(data->assigned(), ", "),
254 absl::StrJoin(data->unassigned(), ", ")));
256 for (int dim_index = 0; dim_index < dims_.size(); ++dim_index) {
259 "InitialProgateDimension(%s)", dims_[dim_index]->DebugString()));
261 dims_[dim_index]->InitialPropagateUnassigned(data->assigned(),
279 DCHECK_EQ(stamp_, solver()->fail_stamp());
280 for (int bin_index = 0; bin_index < bins_; ++bin_index) {
281 if (!removed_[bin_index].empty() || !forced_[bin_index].empty()) {
284 absl::StrFormat("Pack(bin %d, forced = [%s], removed = [%s])",
285 bin_index, absl::StrJoin(forced_[bin_index], ", "),
286 absl::StrJoin(removed_[bin_index], ", ")));
289 for (int dim_index = 0; dim_index < dims_.size(); ++dim_index) {
292 "ProgateDimension(%s)", dims_[dim_index]->DebugString()));
294 dims_[dim_index]->Propagate(bin_index, forced_[bin_index],
305 if (!removed_[bins_].empty() || !forced_[bins_].empty()) {
308 absl::StrFormat("Pack(removed = [%s], forced = [%s])",
309 absl::StrJoin(removed_[bins_], ", "),
310 absl::StrJoin(forced_[bins_], ", ")));
313 for (int dim_index = 0; dim_index < dims_.size(); ++dim_index) {
316 "ProgateDimension(%s)", dims_[dim_index]->DebugString()));
318 dims_[dim_index]->PropagateUnassigned(removed_[bins_], forced_[bins_]);
327 for (int dim_index = 0; dim_index < dims_.size(); ++dim_index) {
339 const uint64_t current_stamp = s->fail_stamp();
340 if (stamp_ < current_stamp) {
344 IntVar* const var = vars_[var_index];
345 bool bound = var->Bound();
346 const int64_t oldmin = var->OldMin();
347 const int64_t oldmax = var->OldMax();
348 const int64_t vmin = var->Min();
349 const int64_t vmax = var->Max();
350 for (int64_t value = std::max(oldmin, int64_t{0});
351 value < std::min(vmin, bins_ + int64_t{1}); ++value) {
352 if (unprocessed_->IsSet(value, var_index)) {
353 unprocessed_->SetToZero(s, value, var_index);
354 removed_[value].push_back(var_index);
358 for (const int64_t value : InitAndGetValues(holes_[var_index])) {
359 if (value >= std::max(int64_t{0}, vmin) &&
360 value <= std::min(static_cast<int64_t>(bins_), vmax)) {
361 DCHECK(unprocessed_->IsSet(value, var_index));
362 unprocessed_->SetToZero(s, value, var_index);
363 removed_[value].push_back(var_index);
367 for (int64_t value = std::max(vmax + 1, int64_t{0});
368 value <= std::min(oldmax, static_cast<int64_t>(bins_)); ++value) {
369 if (unprocessed_->IsSet(value, var_index)) {
370 unprocessed_->SetToZero(s, value, var_index);
371 removed_[value].push_back(var_index);
375 unprocessed_->SetToZero(s, var->Min(), var_index);
376 forced_[var->Min()].push_back(var_index);
382 std::string result = "Pack([";
383 for (int i = 0; i < vars_.size(); ++i) {
384 result += vars_[i]->DebugString() + " ";
386 result += "], dimensions = [";
387 for (int i = 0; i < dims_.size(); ++i) {
388 result += dims_[i]->DebugString() + " ";
411 to_unset_.push_back(std::make_pair(var_index, bin_index));
419 to_set_.push_back(std::make_pair(var_index, bin_index));
439 to_unset_.push_back(std::make_pair(var_index, bins_));
453bool Pack::IsInProcess() const {
458 int var_index = unprocessed_->GetFirstBit(bin_index, 0);
459 while (var_index != -1 && var_index < vars_.size()) {
461 var_index = var_index == vars_.size() - 1
468 int var_index = unprocessed_->GetFirstBit(bin_index, 0);
469 while (var_index != -1 && var_index < vars_.size()) {
470 Assign(var_index, bin_index);
471 var_index = var_index == vars_.size() - 1
478 int var_index = unprocessed_->GetFirstBit(bin_index, 0);
479 if (var_index != -1 && var_index < vars_.size()) {
480 Assign(var_index, bin_index);
485 int var_index = unprocessed_->GetFirstBit(bins_, 0);
486 while (var_index != -1 && var_index < vars_.size()) {
488 var_index = var_index == vars_.size() - 1
495 int var_index = unprocessed_->GetFirstBit(bins_, 0);
496 while (var_index != -1 && var_index < vars_.size()) {
498 var_index = var_index == vars_.size() - 1
512 WeightContainer(int i, int64_t w) : index(i), weight(w) {}
513 bool operator<(const WeightContainer& c) const { return (weight < c.weight); }
516void SortWeightVector(std::vector<int>* const indices,
517 std::vector<WeightContainer>* const to_sort) {
518 std::sort(to_sort->begin(), to_sort->end());
519 for (int index = 0; index < to_sort->size(); ++index) {
520 (*indices)[index] = (*to_sort)[index].index;
522 indices->resize(to_sort->size());
525void SortIndexByWeight(std::vector<int>* const indices,
526 absl::Span<const int64_t> weights) {
527 std::vector<WeightContainer> to_sort;
528 for (int index = 0; index < indices->size(); ++index) {
529 if (weights[index] != 0) {
530 to_sort.push_back(WeightContainer((*indices)[index], weights[index]));
533 SortWeightVector(indices, &to_sort);
536void SortIndexByWeight(std::vector<int>* const indices,
538 std::vector<WeightContainer> to_sort;
539 for (int index = 0; index < indices->size(); ++index) {
540 const int w = weights(index);
541 if (w != 0) {
542 to_sort.push_back(WeightContainer((*indices)[index], w));
545 SortWeightVector(indices, &to_sort);
548void SortIndexByWeight(std::vector<int>* const indices,
550 std::vector<WeightContainer> to_sort;
551 for (int index = 0; index < indices->size(); ++index) {
552 const int w = weights(index, bin_index);
553 if (w != 0) {
554 to_sort.push_back(WeightContainer((*indices)[index], w));
557 SortWeightVector(indices, &to_sort);
560class DimensionLessThanConstant : public Dimension {
562 DimensionLessThanConstant(Solver* const s, Pack* const p,
563 const std::vector<int64_t>& weights,
564 const std::vector<int64_t>& upper_bounds)
566 vars_count_(weights.size()),
568 bins_count_(upper_bounds.size()),
569 upper_bounds_(upper_bounds),
570 first_unbound_backward_vector_(bins_count_, 0),
571 sum_of_bound_variables_vector_(bins_count_, 0LL),
573 for (int i = 0; i < vars_count_; ++i) {
576 SortIndexByWeight(&ranked_, weights_);
579 ~DimensionLessThanConstant() override {}
583 void PushFromTop(int bin_index) {
585 upper_bounds_[bin_index] - sum_of_bound_variables_vector_[bin_index];
589 int last_unbound = first_unbound_backward_vector_[bin_index];
590 for (; last_unbound >= 0; --last_unbound) {
591 const int var_index = ranked_[last_unbound];
592 if (IsUndecided(var_index, bin_index)) {
593 if (weights_[var_index] > slack) {
594 SetImpossible(var_index, bin_index);
600 first_unbound_backward_vector_.SetValue(solver(), bin_index, last_unbound);
603 void InitialPropagate(int bin_index, const std::vector<int>& forced,
604 const std::vector<int>& undecided) override {
605 Solver* const s = solver();
607 for (const int value : forced) {
610 sum_of_bound_variables_vector_.SetValue(s, bin_index, sum);
611 first_unbound_backward_vector_.SetValue(s, bin_index, ranked_.size() - 1);
615 void EndInitialPropagate() override {}
617 void Propagate(int bin_index, const std::vector<int>& forced,
618 const std::vector<int>& removed) override {
620 Solver* const s = solver();
621 int64_t sum = sum_of_bound_variables_vector_[bin_index];
622 for (const int value : forced) {
625 sum_of_bound_variables_vector_.SetValue(s, bin_index, sum);
629 void InitialPropagateUnassigned(const std::vector<int>& assigned,
630 const std::vector<int>& unassigned) override {
632 void PropagateUnassigned(const std::vector<int>& assigned,
633 const std::vector<int>& unassigned) override {}
635 void EndPropagate() override {}
637 void Accept(ModelVisitor* const visitor) const override {
648 const std::vector<int64_t> weights_;
650 const std::vector<int64_t> upper_bounds_;
651 RevArray<int> first_unbound_backward_vector_;
652 RevArray<int64_t> sum_of_bound_variables_vector_;
656class DimensionSumCallbackLessThanConstant : public Dimension {
658 DimensionSumCallbackLessThanConstant(Solver* const s, Pack* const p,
661 const std::vector<int64_t>& upper_bounds)
665 bins_count_(upper_bounds.size()),
666 upper_bounds_(upper_bounds),
667 first_unbound_backward_vector_(bins_count_, 0),
668 sum_of_bound_variables_vector_(bins_count_, 0LL),
672 for (int i = 0; i < vars_count_; ++i) {
675 SortIndexByWeight(&ranked_, weights_);
678 ~DimensionSumCallbackLessThanConstant() override {}
682 void PushFromTop(int bin_index) {
684 upper_bounds_[bin_index] - sum_of_bound_variables_vector_[bin_index];
688 int last_unbound = first_unbound_backward_vector_[bin_index];
689 for (; last_unbound >= 0; --last_unbound) {
690 const int var_index = ranked_[last_unbound];
691 if (IsUndecided(var_index, bin_index)) {
692 if (weights_(var_index) > slack) {
693 SetImpossible(var_index, bin_index);
699 first_unbound_backward_vector_.SetValue(solver(), bin_index, last_unbound);
702 void InitialPropagate(int bin_index, const std::vector<int>& forced,
703 const std::vector<int>& undecided) override {
704 Solver* const s = solver();
706 for (const int value : forced) {
709 sum_of_bound_variables_vector_.SetValue(s, bin_index, sum);
710 first_unbound_backward_vector_.SetValue(s, bin_index, ranked_.size() - 1);
714 void EndInitialPropagate() override {}
716 void Propagate(int bin_index, const std::vector<int>& forced,
717 const std::vector<int>& removed) override {
719 Solver* const s = solver();
720 int64_t sum = sum_of_bound_variables_vector_[bin_index];
721 for (const int value : forced) {
724 sum_of_bound_variables_vector_.SetValue(s, bin_index, sum);
728 void InitialPropagateUnassigned(const std::vector<int>& assigned,
729 const std::vector<int>& unassigned) override {
731 void PropagateUnassigned(const std::vector<int>& assigned,
732 const std::vector<int>& unassigned) override {}
734 void EndPropagate() override {}
736 void Accept(ModelVisitor* const visitor) const override {
750 const std::vector<int64_t> upper_bounds_;
751 RevArray<int> first_unbound_backward_vector_;
752 RevArray<int64_t> sum_of_bound_variables_vector_;
756class DimensionLessThanConstantCallback2 : public Dimension {
758 DimensionLessThanConstantCallback2(Solver* const s, Pack* const p,
761 const std::vector<int64_t>& upper_bounds)
765 bins_count_(upper_bounds.size()),
766 upper_bounds_(upper_bounds),
767 first_unbound_backward_vector_(bins_count_, 0),
768 sum_of_bound_variables_vector_(bins_count_, 0LL),
772 for (int b = 0; b < bins_count_; ++b) {
773 ranked_[b].resize(vars_count);
774 for (int i = 0; i < vars_count_; ++i) {
777 SortIndexByWeight(&ranked_[b], weights_, b);
781 ~DimensionLessThanConstantCallback2() override {}
785 void PushFromTop(int bin_index) {
787 upper_bounds_[bin_index] - sum_of_bound_variables_vector_[bin_index];
791 int last_unbound = first_unbound_backward_vector_[bin_index];
792 for (; last_unbound >= 0; --last_unbound) {
793 const int var_index = ranked_[bin_index][last_unbound];
794 if (IsUndecided(var_index, bin_index)) {
795 if (weights_(var_index, bin_index) > slack) {
796 SetImpossible(var_index, bin_index);
802 first_unbound_backward_vector_.SetValue(solver(), bin_index, last_unbound);
805 void InitialPropagate(int bin_index, const std::vector<int>& forced,
806 const std::vector<int>& undecided) override {
807 Solver* const s = solver();
809 for (const int value : forced) {
810 sum += weights_(value, bin_index);
812 sum_of_bound_variables_vector_.SetValue(s, bin_index, sum);
813 first_unbound_backward_vector_.SetValue(s, bin_index,
814 ranked_[bin_index].size() - 1);
818 void EndInitialPropagate() override {}
820 void Propagate(int bin_index, const std::vector<int>& forced,
821 const std::vector<int>& removed) override {
823 Solver* const s = solver();
824 int64_t sum = sum_of_bound_variables_vector_[bin_index];
825 for (const int value : forced) {
826 sum += weights_(value, bin_index);
828 sum_of_bound_variables_vector_.SetValue(s, bin_index, sum);
832 void InitialPropagateUnassigned(const std::vector<int>& assigned,
833 const std::vector<int>& unassigned) override {
835 void PropagateUnassigned(const std::vector<int>& assigned,
836 const std::vector<int>& unassigned) override {}
838 void EndPropagate() override {}
840 void Accept(ModelVisitor* const visitor) const override {
854 const std::vector<int64_t> upper_bounds_;
855 RevArray<int> first_unbound_backward_vector_;
856 RevArray<int64_t> sum_of_bound_variables_vector_;
857 std::vector<std::vector<int>> ranked_;
860class DimensionWeightedSumEqVar : public Dimension {
862 class VarDemon : public Demon {
864 VarDemon(DimensionWeightedSumEqVar* const dim, int index)
865 : dim_(dim), index_(index) {}
868 void Run(Solver* const s) override { dim_->PushFromTop(index_); }
871 DimensionWeightedSumEqVar* const dim_;
875 DimensionWeightedSumEqVar(Solver* const s, Pack* const p,
876 const std::vector<int64_t>& weights,
877 const std::vector<IntVar*>& loads)
879 vars_count_(weights.size()),
881 bins_count_(loads.size()),
883 first_unbound_backward_vector_(bins_count_, 0),
884 sum_of_bound_variables_vector_(bins_count_, 0LL),
885 sum_of_all_variables_vector_(bins_count_, 0LL),
887 DCHECK_GT(vars_count_, 0);
888 DCHECK_GT(bins_count_, 0);
889 for (int i = 0; i < vars_count_; ++i) {
892 SortIndexByWeight(&ranked_, weights_);
895 ~DimensionWeightedSumEqVar() override {}
897 std::string DebugString() const override {
898 return "DimensionWeightedSumEqVar";
902 for (int i = 0; i < bins_count_; ++i) {
903 Demon* const d = solver()->RevAlloc(new VarDemon(this, i));
904 loads_[i]->WhenRange(d);
908 void PushFromTop(int bin_index) {
909 IntVar* const load = loads_[bin_index];
910 const int64_t sum_min = sum_of_bound_variables_vector_[bin_index];
911 const int64_t sum_max = sum_of_all_variables_vector_[bin_index];
912 load->SetRange(sum_min, sum_max);
913 const int64_t slack_up = load->Max() - sum_min;
914 const int64_t slack_down = sum_max - load->Min();
917 int last_unbound = first_unbound_backward_vector_[bin_index];
918 for (; last_unbound >= 0; --last_unbound) {
919 const int var_index = ranked_[last_unbound];
920 const int64_t weight = weights_[var_index];
921 if (IsUndecided(var_index, bin_index)) {
923 SetImpossible(var_index, bin_index);
924 } else if (weight > slack_down) {
925 Assign(var_index, bin_index);
931 first_unbound_backward_vector_.SetValue(solver(), bin_index, last_unbound);
934 void InitialPropagate(int bin_index, const std::vector<int>& forced,
935 const std::vector<int>& undecided) override {
936 Solver* const s = solver();
938 for (const int value : forced) {
941 sum_of_bound_variables_vector_.SetValue(s, bin_index, sum);
942 for (const int value : undecided) {
945 sum_of_all_variables_vector_.SetValue(s, bin_index, sum);
946 first_unbound_backward_vector_.SetValue(s, bin_index, ranked_.size() - 1);
950 void EndInitialPropagate() override {}
952 void Propagate(int bin_index, const std::vector<int>& forced,
953 const std::vector<int>& removed) override {
954 Solver* const s = solver();
955 int64_t down = sum_of_bound_variables_vector_[bin_index];
956 for (const int value : forced) {
959 sum_of_bound_variables_vector_.SetValue(s, bin_index, down);
960 int64_t up = sum_of_all_variables_vector_[bin_index];
961 for (const int value : removed) {
964 sum_of_all_variables_vector_.SetValue(s, bin_index, up);
967 void InitialPropagateUnassigned(const std::vector<int>& assigned,
968 const std::vector<int>& unassigned) override {
970 void PropagateUnassigned(const std::vector<int>& assigned,
971 const std::vector<int>& unassigned) override {}
973 void EndPropagate() override {}
975 void Accept(ModelVisitor* const visitor) const override {
986 const std::vector<int64_t> weights_;
988 const std::vector<IntVar*> loads_;
989 RevArray<int> first_unbound_backward_vector_;
990 RevArray<int64_t> sum_of_bound_variables_vector_;
991 RevArray<int64_t> sum_of_all_variables_vector_;
995class DimensionWeightedCallback2SumEqVar : public Dimension {
997 class VarDemon : public Demon {
999 VarDemon(DimensionWeightedCallback2SumEqVar* const dim, int index)
1000 : dim_(dim), index_(index) {}
1003 void Run(Solver* const s) override { dim_->PushFromTop(index_); }
1006 DimensionWeightedCallback2SumEqVar* const dim_;
1010 DimensionWeightedCallback2SumEqVar(Solver* const s, Pack* const p,
1013 const std::vector<IntVar*>& loads)
1017 bins_count_(loads.size()),
1019 first_unbound_backward_vector_(bins_count_, 0),
1020 sum_of_bound_variables_vector_(bins_count_, 0LL),
1021 sum_of_all_variables_vector_(bins_count_, 0LL),
1024 DCHECK_GT(vars_count_, 0);
1025 DCHECK_GT(bins_count_, 0);
1026 for (int b = 0; b < bins_count_; ++b) {
1027 ranked_[b].resize(vars_count);
1028 for (int i = 0; i < vars_count_; ++i) {
1031 SortIndexByWeight(&ranked_[b], weights_, b);
1035 ~DimensionWeightedCallback2SumEqVar() override {}
1037 std::string DebugString() const override {
1038 return "DimensionWeightedCallback2SumEqVar";
1042 for (int i = 0; i < bins_count_; ++i) {
1043 Demon* const d = solver()->RevAlloc(new VarDemon(this, i));
1044 loads_[i]->WhenRange(d);
1048 void PushFromTop(int bin_index) {
1049 IntVar* const load = loads_[bin_index];
1050 const int64_t sum_min = sum_of_bound_variables_vector_[bin_index];
1051 const int64_t sum_max = sum_of_all_variables_vector_[bin_index];
1052 load->SetRange(sum_min, sum_max);
1053 const int64_t slack_up = load->Max() - sum_min;
1054 const int64_t slack_down = sum_max - load->Min();
1055 DCHECK_GE(slack_down, 0);
1057 int last_unbound = first_unbound_backward_vector_[bin_index];
1058 for (; last_unbound >= 0; --last_unbound) {
1059 const int var_index = ranked_[bin_index][last_unbound];
1060 const int64_t weight = weights_(var_index, bin_index);
1061 if (IsUndecided(var_index, bin_index)) {
1063 SetImpossible(var_index, bin_index);
1064 } else if (weight > slack_down) {
1065 Assign(var_index, bin_index);
1071 first_unbound_backward_vector_.SetValue(solver(), bin_index, last_unbound);
1074 void InitialPropagate(int bin_index, const std::vector<int>& forced,
1075 const std::vector<int>& undecided) override {
1076 Solver* const s = solver();
1078 for (const int value : forced) {
1079 sum += weights_(value, bin_index);
1081 sum_of_bound_variables_vector_.SetValue(s, bin_index, sum);
1082 for (const int value : undecided) {
1083 sum += weights_(value, bin_index);
1085 sum_of_all_variables_vector_.SetValue(s, bin_index, sum);
1086 first_unbound_backward_vector_.SetValue(s, bin_index,
1087 ranked_[bin_index].size() - 1);
1091 void EndInitialPropagate() override {}
1093 void Propagate(int bin_index, const std::vector<int>& forced,
1094 const std::vector<int>& removed) override {
1095 Solver* const s = solver();
1096 int64_t down = sum_of_bound_variables_vector_[bin_index];
1097 for (const int value : forced) {
1098 down += weights_(value, bin_index);
1100 sum_of_bound_variables_vector_.SetValue(s, bin_index, down);
1101 int64_t up = sum_of_all_variables_vector_[bin_index];
1102 for (const int value : removed) {
1103 up -= weights_(value, bin_index);
1105 sum_of_all_variables_vector_.SetValue(s, bin_index, up);
1108 void InitialPropagateUnassigned(const std::vector<int>& assigned,
1109 const std::vector<int>& unassigned) override {
1111 void PropagateUnassigned(const std::vector<int>& assigned,
1112 const std::vector<int>& unassigned) override {}
1114 void EndPropagate() override {}
1116 void Accept(ModelVisitor* const visitor) const override {
1130 const std::vector<IntVar*> loads_;
1131 RevArray<int> first_unbound_backward_vector_;
1132 RevArray<int64_t> sum_of_bound_variables_vector_;
1133 RevArray<int64_t> sum_of_all_variables_vector_;
1134 std::vector<std::vector<int>> ranked_;
1137class AssignedWeightedSumDimension : public Dimension {
1139 class VarDemon : public Demon {
1141 explicit VarDemon(AssignedWeightedSumDimension* const dim) : dim_(dim) {}
1144 void Run(Solver* const s) override { dim_->PropagateAll(); }
1147 AssignedWeightedSumDimension* const dim_;
1150 AssignedWeightedSumDimension(Solver* const s, Pack* const p,
1151 const std::vector<int64_t>& weights,
1152 int bins_count, IntVar* const cost_var)
1154 vars_count_(weights.size()),
1158 first_unbound_backward_(0),
1159 sum_of_assigned_items_(0LL),
1160 sum_of_unassigned_items_(0LL),
1164 DCHECK_GT(vars_count_, 0);
1165 DCHECK_GT(bins_count_, 0);
1166 for (int i = 0; i < vars_count_; ++i) {
1169 SortIndexByWeight(&ranked_, weights_);
1170 first_unbound_backward_.SetValue(s, ranked_.size() - 1);
1173 ~AssignedWeightedSumDimension() override {}
1176 Demon* const uv = solver()->RevAlloc(new VarDemon(this));
1177 cost_var_->WhenRange(uv);
1181 cost_var_->SetRange(sum_of_assigned_items_.Value(),
1182 sum_all_weights_ - sum_of_unassigned_items_.Value());
1183 const int64_t slack_up = cost_var_->Max() - sum_of_assigned_items_.Value();
1184 const int64_t slack_down = sum_all_weights_ - cost_var_->Min();
1185 int last_unbound = first_unbound_backward_.Value();
1186 for (; last_unbound >= 0; --last_unbound) {
1187 const int var_index = ranked_[last_unbound];
1188 if (!IsAssignedStatusKnown(var_index)) {
1189 const int64_t coefficient = weights_[var_index];
1190 if (coefficient > slack_up) {
1191 SetUnassigned(var_index);
1192 } else if (coefficient > slack_down) {
1199 first_unbound_backward_.SetValue(solver(), last_unbound);
1202 void InitialPropagate(int bin_index, const std::vector<int>& forced,
1203 const std::vector<int>& undecided) override {}
1205 void EndInitialPropagate() override {}
1207 void InitialPropagateUnassigned(const std::vector<int>& assigned,
1208 const std::vector<int>& unassigned) override {
1209 for (int index = 0; index < vars_count_; ++index) {
1210 sum_all_weights_ += weights_[index];
1213 PropagateUnassigned(assigned, unassigned);
1216 void Propagate(int bin_index, const std::vector<int>& forced,
1217 const std::vector<int>& removed) override {}
1219 void PropagateUnassigned(const std::vector<int>& assigned,
1220 const std::vector<int>& unassigned) override {
1221 int64_t sum_assigned = sum_of_assigned_items_.Value();
1222 for (int index = 0; index < assigned.size(); ++index) {
1223 const int var_index = assigned[index];
1224 sum_assigned += weights_[var_index];
1227 int64_t sum_unassigned = sum_of_unassigned_items_.Value();
1228 for (int index = 0; index < unassigned.size(); ++index) {
1229 const int var_index = unassigned[index];
1230 sum_unassigned += weights_[var_index];
1233 Solver* const s = solver();
1234 sum_of_assigned_items_.SetValue(s, sum_assigned);
1235 sum_of_unassigned_items_.SetValue(s, sum_unassigned);
1239 void EndPropagate() override {}
1241 void Accept(ModelVisitor* const visitor) const override {
1242 visitor->BeginVisitExtension(
1248 visitor->EndVisitExtension(
1254 const std::vector<int64_t> weights_;
1257 Rev<int> first_unbound_backward_;
1258 Rev<int64_t> sum_of_assigned_items_;
1259 Rev<int64_t> sum_of_unassigned_items_;
1260 std::vector<int> ranked_;
1261 int64_t sum_all_weights_;
1266class CountAssignedItemsDimension : public Dimension {
1268 class VarDemon : public Demon {
1270 explicit VarDemon(CountAssignedItemsDimension* const dim) : dim_(dim) {}
1273 void Run(Solver* const s) override { dim_->PropagateAll(); }
1276 CountAssignedItemsDimension* const dim_;
1279 CountAssignedItemsDimension(Solver* const s, Pack* const p, int vars_count,
1280 int bins_count, IntVar* const cost_var)
1285 first_unbound_backward_(0),
1289 DCHECK_GT(vars_count, 0);
1290 DCHECK_GT(bins_count, 0);
1293 ~CountAssignedItemsDimension() override {}
1296 Demon* const uv = solver()->RevAlloc(new VarDemon(this));
1297 cost_var_->WhenRange(uv);
1301 cost_var_->SetRange(assigned_count_.Value(),
1302 vars_count_ - unassigned_count_.Value());
1303 if (assigned_count_.Value() == cost_var_->Max()) {
1304 UnassignAllRemainingItems();
1305 } else if (cost_var_->Min() == vars_count_ - unassigned_count_.Value()) {
1306 AssignAllRemainingItems();
1310 void InitialPropagate(int bin_index, const std::vector<int>& forced,
1311 const std::vector<int>& undecided) override {}
1313 void EndInitialPropagate() override {}
1315 void InitialPropagateUnassigned(const std::vector<int>& assigned,
1316 const std::vector<int>& unassigned) override {
1317 PropagateUnassigned(assigned, unassigned);
1320 void Propagate(int bin_index, const std::vector<int>& forced,
1321 const std::vector<int>& removed) override {}
1323 void PropagateUnassigned(const std::vector<int>& assigned,
1324 const std::vector<int>& unassigned) override {
1325 Solver* const s = solver();
1326 assigned_count_.SetValue(s, assigned_count_.Value() + assigned.size());
1327 unassigned_count_.SetValue(s,
1328 unassigned_count_.Value() + unassigned.size());
1332 void EndPropagate() override {}
1334 void Accept(ModelVisitor* const visitor) const override {
1345 Rev<int> first_unbound_backward_;
1346 Rev<int> assigned_count_;
1347 Rev<int> unassigned_count_;
1352class CountUsedBinDimension : public Dimension {
1354 class VarDemon : public Demon {
1356 explicit VarDemon(CountUsedBinDimension* const dim) : dim_(dim) {}
1359 void Run(Solver* const s) override { dim_->PropagateAll(); }
1362 CountUsedBinDimension* const dim_;
1365 CountUsedBinDimension(Solver* const s, Pack* const p, int vars_count,
1366 int bins_count, IntVar* const count_var)
1372 candidates_(bins_count_, 0),
1378 DCHECK_GT(vars_count, 0);
1379 DCHECK_GT(bins_count, 0);
1382 ~CountUsedBinDimension() override {}
1385 Demon* const uv = solver()->RevAlloc(new VarDemon(this));
1386 count_var_->WhenRange(uv);
1388 initial_max_ = bins_count_;
1392 count_var_->SetRange(card_min_.Value(), card_max_.Value());
1393 if (card_min_.Value() == count_var_->Max()) {
1394 for (int bin_index = 0; bin_index < bins_count_; ++bin_index) {
1395 if (!used_.IsSet(bin_index) && candidates_[bin_index] > 0) {
1396 RemoveAllPossibleFromBin(bin_index);
1399 } else if (card_max_.Value() == count_var_->Min()) {
1400 for (int bin_index = 0; bin_index < bins_count_; ++bin_index) {
1401 if (candidates_[bin_index] == 1) {
1402 AssignFirstPossibleToBin(bin_index);
1408 void InitialPropagate(int bin_index, const std::vector<int>& forced,
1409 const std::vector<int>& undecided) override {
1411 used_.SetToOne(solver(), bin_index);
1413 } else if (!undecided.empty()) {
1414 candidates_.SetValue(solver(), bin_index, undecided.size());
1420 void EndInitialPropagate() override {
1421 card_min_.SetValue(solver(), initial_min_);
1422 card_max_.SetValue(solver(), initial_max_);
1426 void InitialPropagateUnassigned(const std::vector<int>& assigned,
1427 const std::vector<int>& unassigned) override {
1430 void Propagate(int bin_index, const std::vector<int>& forced,
1431 const std::vector<int>& removed) override {
1432 if (!used_.IsSet(bin_index)) {
1434 used_.SetToOne(solver(), bin_index);
1435 card_min_.SetValue(solver(), card_min_.Value() + 1);
1436 } else if (!removed.empty()) {
1437 candidates_.SetValue(solver(), bin_index,
1438 candidates_.Value(bin_index) - removed.size());
1439 if (candidates_[bin_index] == 0) {
1440 card_max_.SetValue(solver(), card_max_.Value() - 1);
1446 void PropagateUnassigned(const std::vector<int>& assigned,
1447 const std::vector<int>& unassigned) override {}
1449 void EndPropagate() override { PropagateAll(); }
1451 void Accept(ModelVisitor* const visitor) const override {
1461 IntVar* const count_var_;
1463 RevArray<int> candidates_;
1473class VariableUsageDimension : public Dimension {
1475 VariableUsageDimension(Solver* const solver, Pack* const pack,
1476 const std::vector<int64_t>& capacities,
1477 const std::vector<IntVar*>& weights)
1478 : Dimension(solver, pack), capacities_(capacities), weights_(weights) {}
1480 ~VariableUsageDimension() override {}
1483 Solver* const s = solver();
1484 const int num_bins = capacities_.size();
1485 const int num_items = weights_.size();
1487 for (int bin_index = 0; bin_index < num_bins; ++bin_index) {
1488 std::vector<IntVar*> terms;
1489 for (int item_index = 0; item_index < num_items; ++item_index) {
1490 IntVar* const assign_var = AssignVar(item_index, bin_index);
1491 terms.push_back(s->MakeProd(assign_var, weights_[item_index])->Var());
1493 s->AddConstraint(s->MakeSumLessOrEqual(terms, capacities_[bin_index]));
1497 void InitialPropagate(int bin_index, const std::vector<int>& forced,
1498 const std::vector<int>& undecided) override {}
1499 void InitialPropagateUnassigned(const std::vector<int>& assigned,
1500 const std::vector<int>& unassigned) override {
1502 void EndInitialPropagate() override {}
1503 void Propagate(int bin_index, const std::vector<int>& forced,
1504 const std::vector<int>& removed) override {}
1505 void PropagateUnassigned(const std::vector<int>& assigned,
1506 const std::vector<int>& unassigned) override {}
1507 void EndPropagate() override {}
1509 std::string DebugString() const override { return "VariableUsageDimension"; }
1511 void Accept(ModelVisitor* const visitor) const override {
1512 visitor->BeginVisitExtension(
1518 visitor->EndVisitExtension(
1523 const std::vector<int64_t> capacities_;
1524 const std::vector<IntVar*> weights_;
1531 const std::vector<int64_t>& weights, const std::vector<int64_t>& bounds) {
1532 CHECK_EQ(weights.size(), vars_.size());
1533 CHECK_EQ(bounds.size(), bins_);
1536 s->RevAlloc(new DimensionLessThanConstant(s, this, weights, bounds));
1542 CHECK(weights != nullptr);
1543 CHECK_EQ(bounds.size(), bins_);
1545 Dimension* const dim = s->RevAlloc(new DimensionSumCallbackLessThanConstant(
1552 CHECK(weights != nullptr);
1553 CHECK_EQ(bounds.size(), bins_);
1555 Dimension* const dim = s->RevAlloc(new DimensionLessThanConstantCallback2(
1561 const std::vector<IntVar*>& loads) {
1562 CHECK_EQ(weights.size(), vars_.size());
1563 CHECK_EQ(loads.size(), bins_);
1566 s->RevAlloc(new DimensionWeightedSumEqVar(s, this, weights, loads));
1571 const std::vector<IntVar*>& loads) {
1572 CHECK(weights != nullptr);
1573 CHECK_EQ(loads.size(), bins_);
1575 Dimension* const dim = s->RevAlloc(new DimensionWeightedCallback2SumEqVar(
1581 const std::vector<int64_t>& weights, IntVar* const cost_var) {
1582 CHECK_EQ(weights.size(), vars_.size());
1585 new AssignedWeightedSumDimension(s, this, weights, bins_, cost_var));
1590 const std::vector<IntVar*>& usage, const std::vector<int64_t>& capacity) {
1591 CHECK_EQ(usage.size(), vars_.size());
1592 CHECK_EQ(capacity.size(), bins_);
1595 s->RevAlloc(new VariableUsageDimension(s, this, capacity, usage));
Constraint(Solver *const solver)
void SetUnassigned(int var_index)
Definition pack.cc:84
void AssignAllRemainingItems()
Definition pack.cc:98
void SetAssigned(int var_index)
Definition pack.cc:82
virtual void EndInitialPropagate()=0
bool IsPossible(int var_index, int bin_index) const
Definition pack.cc:62
virtual void EndPropagate()=0
~Dimension() override
Definition pack.cc:40
void AssignFirstPossibleToBin(int bin_index)
Definition pack.cc:94
void UnassignAllRemainingItems()
Definition pack.cc:100
void SetImpossible(int var_index, int bin_index)
Definition pack.cc:70
void AssignAllPossibleToBin(int bin_index)
Definition pack.cc:90
bool IsUndecided(int var_index, int bin_index) const
Definition pack.cc:58
std::string DebugString() const override
Definition pack.cc:53
virtual void PropagateUnassigned(const std::vector< int > &assigned, const std::vector< int > &unassigned)=0
virtual void Propagate(int bin_index, const std::vector< int > &forced, const std::vector< int > &removed)=0
virtual void InitialPropagate(int bin_index, const std::vector< int > &forced, const std::vector< int > &undecided)=0
void RemoveAllPossibleFromBin(int bin_index)
Definition pack.cc:86
IntVar * AssignVar(int var_index, int bin_index) const
Definition pack.cc:66
void Assign(int var_index, int bin_index)
Definition pack.cc:74
virtual void Accept(ModelVisitor *visitor) const =0
Solver * solver() const
Definition pack.cc:56
bool IsAssignedStatusKnown(int var_index) const
Definition pack.cc:78
virtual void InitialPropagateUnassigned(const std::vector< int > &assigned, const std::vector< int > &unassigned)=0
Dimension(Solver *const s, Pack *const pack)
Definition pack.cc:38
virtual bool Bound() const
Returns true if the min and the max of the expression are equal.
virtual void SetRange(int64_t l, int64_t u)
This method sets both the min and the max of the expression.
virtual int64_t Min() const =0
virtual int64_t Max() const =0
virtual void WhenDomain(Demon *d)=0
virtual IntVarIterator * MakeDomainIterator(bool reversible) const =0
virtual int64_t OldMax() const =0
Returns the previous max.
virtual int64_t OldMin() const =0
Returns the previous min.
static const char kVariableUsageLessConstantExtension[]
virtual void VisitIntegerArgument(const std::string &arg_name, int64_t value)
Visit integer arguments.
static const char kCoefficientsArgument[]
static const char kTargetArgument[]
static const char kSizeArgument[]
virtual void VisitIntegerVariableArrayArgument(const std::string &arg_name, const std::vector< IntVar * > &arguments)
static const char kValuesArgument[]
static const char kPack[]
static const char kUsageLessConstantExtension[]
static const char kVarsArgument[]
virtual void BeginVisitConstraint(const std::string &type_name, const Constraint *constraint)
static const char kUsageEqualVariableExtension[]
static const char kCountAssignedItemsExtension[]
Extension names:
static const char kCountUsedBinsExtension[]
static const char kWeightedSumOfAssignedEqualVariableExtension[]
virtual void EndVisitConstraint(const std::string &type_name, const Constraint *constraint)
void SetUnassigned(int var_index)
Definition pack.cc:445
void SetAssigned(int var_index)
Definition pack.cc:437
void Assign(int var_index, int bin_index)
Definition pack.cc:417
bool IsPossible(int var_index, int bin_index) const
Definition pack.cc:429
void Post() override
Definition pack.cc:128
void AddCountUsedBinDimension(IntVar *count_var)
Definition pack.cc:1599
void Propagate()
Definition pack.cc:276
Pack(Solver *s, const std::vector< IntVar * > &vars, int number_of_bins)
Definition pack.cc:109
void AddWeightedSumLessOrEqualConstantDimension(const std::vector< int64_t > &weights, const std::vector< int64_t > &bounds)
Definition pack.cc:1530
void AddSumVariableWeightsLessOrEqualConstantDimension(const std::vector< IntVar * > &usage, const std::vector< int64_t > &capacity)
Definition pack.cc:1589
void AddCountAssignedItemsDimension(IntVar *count_var)
Definition pack.cc:1606
void UnassignAllRemainingItems()
Definition pack.cc:494
void AssignAllRemainingItems()
Definition pack.cc:484
std::string DebugString() const override
Definition pack.cc:381
void AssignAllPossibleToBin(int bin_index)
Definition pack.cc:467
IntVar * AssignVar(int var_index, int bin_index) const
Definition pack.cc:433
bool IsAssignedStatusKnown(int var_index) const
Definition pack.cc:425
void AssignFirstPossibleToBin(int bin_index)
Definition pack.cc:477
void AddWeightedSumEqualVarDimension(const std::vector< int64_t > &weights, const std::vector< IntVar * > &loads)
Definition pack.cc:1560
void RemoveAllPossibleFromBin(int bin_index)
Definition pack.cc:457
void PropagateDelayed()
Definition pack.cc:155
~Pack() override
Definition pack.cc:126
void AddWeightedSumOfAssignedDimension(const std::vector< int64_t > &weights, IntVar *cost_var)
Definition pack.cc:1580
void InitialPropagate() override
Definition pack.cc:191
void SetImpossible(int var_index, int bin_index)
Definition pack.cc:409
void ClearAll()
Definition pack.cc:144
void OneDomain(int var_index)
Definition pack.cc:335
bool IsUndecided(int var_index, int bin_index) const
Definition pack.cc:405
void EnqueueDelayedDemon(Demon *const d)
virtual void PopContext()=0
virtual void PushContext(const std::string &context)=0
Matrix version of the RevBitSet class.
Demon * RegisterDemon(Demon *demon)
Adds a new demon and wraps it inside a DemonProfiler if necessary.
uint64_t fail_stamp() const
The fail_stamp() is incremented after each backtrack.
IntVar * MakeIsEqualCstVar(IntExpr *var, int64_t value)
status var of (var == value)
std::function< int64_t(int64_t)> IndexEvaluator1
Callback typedefs.
Pack * MakePack(const std::vector< IntVar * > &vars, int number_of_bins)
Definition pack.cc:1613
std::function< int64_t(int64_t, int64_t)> IndexEvaluator2
PropagationMonitor * GetPropagationMonitor() const
Returns the propagation monitor.
bool InstrumentsVariables() const
Returns whether we are tracing variables.
Demon * MakeDelayedConstraintDemon0(Solver *const s, T *const ct, void(T::*method)(), const std::string &name)
Demon * MakeConstraintDemon1(Solver *const s, T *const ct, void(T::*method)(P), const std::string &name, P param1)
trees with all degrees equal w the current value of w