86 unsigned int _norders;
96 virtual void help(
void);
98 bool parse(
int& argc,
char* argv[]);
101 unsigned int size(
void)
const {
return _size; }
109 int*
loss(
void)
const {
return _loss; }
126 bool operator() (
int i,
int j) {
198 capacities(opt.capacities()), ncapacities(opt.ncapacities()),
199 maxcapacity(opt.maxcapacity()), loss(opt.loss()),
200 ncolors(opt.ncolors()), orders(opt.orders()),
201 norders(opt.
size()), nslabs(opt.
size()),
203 slab(*this, norders, 0,nslabs-1),
204 slabload(*this, nslabs, 0,45),
205 slabcost(*this, nslabs, 0, Int::Limits::
max),
206 total_cost(*this, 0, Int::Limits::
max)
210 for (
unsigned int i = 0;
i < norders; ++
i) {
212 for (
int j = nslabs; j--; ) {
213 boolslab[j +
i*nslabs] = tmp[j] =
BoolVar(*
this, 0, 1);
219 for (
unsigned int s = 0; s < nslabs; ++s) {
222 for (
int i = norders;
i--; ) {
224 x[
i] = boolslab[s +
i*nslabs];
230 for (
unsigned int i = norders;
i-- ; )
237 for (
int c = ncolors;
c--; ) {
239 for (
int o = norders; o--; ) {
244 for (
unsigned int s = 0; s < nslabs; ++s) {
246 for (
int c = ncolors;
c--; ) {
250 for (
int o = norders; o--; ) {
252 hasc[pos++] = boolslab[s + o*nslabs];
254 assert(pos == nofcolor[c]);
266 for (
int s = nslabs; s--; ) {
267 element(*
this, l, slabload[s], slabcost[s]);
272 if (opt.
symmetry() == SYMMETRY_BRANCHING) {
275 }
else if (opt.
symmetry() == SYMMETRY_NONE) {
288 for (
unsigned int i = 0 ;
i < norders ;
i++)
292 for (
unsigned int i = 0 ;
i < norders ;
i++) {
293 sorted_orders[
i] = slab[indices[
i]];
302 os <<
"What slab=" << slab << std::endl;
303 os <<
"Slab load=" << slabload << std::endl;
304 os <<
"Slab cost=" << slabcost << std::endl;
305 os <<
"Total cost=" << total_cost << std::endl;
308 bool unassigned =
false;
309 for (
int i = nslabs;
i--; ) {
314 if (slabload[
i].
min()>0) ++nslabsused;
315 if (slabcost[
i].
min()>0) ++nslabscost;
318 os <<
"Number of slabs used=" << nslabsused
319 <<
", slabs with cost=" << nslabscost
327 capacities(s.capacities), ncapacities(s.ncapacities),
328 maxcapacity(s.maxcapacity), loss(s.loss),
329 ncolors(s.ncolors), orders(s.orders),
330 norders(s.norders), nslabs(s.nslabs) {
374 e << alternatives() << pos << val;
383 :
Brancher(home, b), start(b.start) {
390 for (
unsigned int i = start;
i < sm.
norders; ++
i)
405 unsigned int pos = start;
406 for (
unsigned int i = start;
i<sm.
norders; ++
i) {
419 unsigned int val = sm.
slab[
pos].min();
421 unsigned int firstzero = 0;
426 return new Choice(*
this, (val<firstzero) ? 2 : 1, pos, val);
429 unsigned int alt;
int pos, val;
430 e >> alt >> pos >> val;
431 return new Choice(*
this, alt, pos, val);
448 std::ostream& o)
const {
450 o <<
"slab[" << c.
pos <<
"] " 451 << ((a == 0) ?
"=" :
"!=")
464 return sizeof(*this);
480 if (!opt.
parse(argc,argv))
482 Script::run<SteelMill,BAB,SteelMillOptions>(
opt);
490 std::cerr <<
"\t(string), optional" << std::endl
491 <<
"\t\tBenchmark to load." << std::endl
492 <<
"\t\tIf none is given, the standard CSPLib instance is used." 494 std::cerr <<
"\t(unsigned int), optional" << std::endl
495 <<
"\t\tNumber of orders to use, in the interval [0..norders]." 497 <<
"\t\tIf none is given, all orders are used." << std::endl;
505 std::cerr <<
"Too many arguments given, max two allowed (given={";
506 for (
int i = 1;
i < argc; ++
i) {
507 std::cerr <<
"\"" << argv[
i] <<
"\"";
508 if (
i < argc-1) std::cerr <<
",";
510 std::cerr <<
"})." << std::endl;
516 for (
int i = strlen(argv[argc-1]);
i-- && issize; )
517 issize &= (isdigit(argv[argc-1][
i]) != 0);
519 _size = atoi(argv[argc-1]);
521 std::ifstream instance(argv[argc-1]);
522 if (instance.fail()) {
523 std::cerr <<
"Argument \"" << argv[argc-1]
524 <<
"\" is neither an integer nor a readable file" 529 instance >> _ncapacities;
530 _capacities =
new int[_ncapacities];
532 for (
int i = 0;
i < _ncapacities; ++
i) {
533 instance >> _capacities[
i];
534 _maxcapacity =
std::max(_maxcapacity, _capacities[
i]);
536 instance >> _ncolors >> _norders;
537 _orders =
new int[_norders][2];
538 for (
unsigned int i = 0;
i < _norders; ++
i) {
547 _loss =
new int[_maxcapacity+1];
550 for (
int c = 1;
c < _maxcapacity; ++
c) {
551 if (
c > _capacities[currcap]) ++currcap;
552 _loss[
c] = _capacities[currcap] -
c;
560 if (_size == 0 || _size > _norders) {
561 std::cerr <<
"Size must be between 1 and " << _norders << std::endl;
int ncapacities
Number of capacities.
SteelMillBranch(Home home)
Construct brancher.
ModEvent nq(Space &home, int n)
Restrict domain values to be different from n.
static IntArgs create(int n, int start, int inc=1)
Allocate array with n elements such that for all .
IntVarArray slab
Slab assigned to order i.
int val
Value of variable.
order_t orders
The orders.
unsigned int csplib_maxcapacity
Maximum capacity.
int * loss
Loss for all sizes.
int * capacities
Capacities.
void post(Home home, Term *t, int n, FloatRelType frt, FloatVal c)
Post propagator for linear constraint over floats.
void branch(Home home, const FloatVarArgs &x, FloatVarBranch vars, FloatValBranch vals, FloatBranchFilter bf, FloatVarValPrint vvp)
Branch over x with variable selection vars and value selection vals.
void channel(Home home, FloatVar x0, IntVar x1)
Post propagator for channeling a float and an integer variable .
int size(void) const
Return size of array (number of elements)
const FloatNum max
Largest allowed float value.
void update(Space &home, VarArray< Var > &a)
Update array to be a clone of array a.
virtual bool status(const Space &home) const
Check status of brancher, return true if alternatives left.
ModEvent eq(Space &home, int n)
Restrict domain values to be equal to n.
unsigned int norders
Number of orders.
order_t orders(void) const
Return orders.
void max(Home home, FloatVar x0, FloatVar x1, FloatVar x2)
Post propagator for .
GECODE_FLATZINC_EXPORT FlatZincSpace * parse(const std::string &fileName, Printer &p, std::ostream &err=std::cerr, FlatZincSpace *fzs=NULL, Rnd &rnd=defrnd)
Parse FlatZinc file fileName into fzs and return it.
int ncolors
Number of colors.
int ncapacities(void) const
Return number of capacities.
bool assigned(void) const
Test if all variables are assigned.
bool pos(const View &x)
Test whether x is postive.
virtual void help(void)
Print help text.
int maxcapacity
Maximum capacity.
Collection of symmetries.
void linear(Home home, const FloatVarArgs &x, FloatRelType frt, FloatVal c)
Post propagator for .
Breaking symmetries with symmetry.
SortByWeight(order_t _orders)
Initialize orders.
int start
Cache of first unassigned value.
Multi _c(Gecode::IntArgs(3, 1, 2, 3))
const int order_weight
Weight-position in order-array elements.
unsigned int csplib_norders
Number of orders.
IntVarArray slabcost
Cost of slab j.
Parametric base-class for scripts.
Base-class for both propagators and branchers.
SteelMillOptions(const char *n)
Initialize options for example with name n.
virtual void print(std::ostream &os) const
Print solution.
virtual IntVar cost(void) const
Return solution cost.
Gecode::FloatVal c(-8, 8)
virtual Space * copy(void)
Copy during cloning.
Gecode::IntArgs i(4, 1, 2, 3, 4)
Base-class for branchers.
int maxcapacity(void) const
Return maximum of capacities.
void quicksort(Type *l, Type *r, Less &less)
Standard quick sort.
int n
Number of negative literals for node type.
static void post(Home home)
Post brancher.
Execution has resulted in failure.
int csplib_orders[][2]
Orders.
virtual void archive(Archive &e) const
Archive into e.
virtual ExecStatus commit(Space &home, const Gecode::Choice &_c, unsigned int a)
Perform commit for choice _c and alternative a.
virtual Actor * copy(Space &home)
Copy brancher.
int norders(void) const
Return number of orders.
IntValBranch INT_VAL_MIN(void)
Select smallest value.
unsigned int size(I &i)
Size of all ranges of range iterator i.
SymmetryHandle ValueSymmetry(const IntArgs &vs)
Values in v are interchangeable.
IntVar total_cost
Total cost.
int * loss(void) const
Return loss values.
virtual Choice * choice(const Space &, Archive &e)
Return choice from e.
void update(Space &home, VarImpVar< VarImp > &y)
Update this variable to be a clone of variable y.
int pos
Position of variable.
int csplib_capacities[]
Constants for CSPLib instance of the Steel Mill Slab Design Problem.
const int order_color
Color-position in order-array elements.
virtual Gecode::Choice * choice(Space &home)
Return choice.
int(* order_t)[2]
Order-specifications.
unsigned int nslabs
Number of slabs.
Passing integer variables.
Passing integer arguments.
Passing Boolean variables.
int csplib_loss[]
Loss for all sizes.
Boolean integer variables.
Choice(const Brancher &b, unsigned int a, int pos0, int val0)
void parse(int &argc, char *argv[])
Parse options from arguments argv (number is argc)
bool parse(int &argc, char *argv[])
Parse options from arguments argv (number is argc)
void min(Home home, FloatVar x0, FloatVar x1, FloatVar x2)
Post propagator for .
struct Gecode::@585::NNF::@62::@63 b
For binary nodes (and, or, eqv)
IntVarBranch INT_VAR_SIZE_MIN(BranchTbl tbl)
Select variable with smallest domain size.
Post propagator for f(x \diamond_{\mathit{op}} y) \sim_r z \f$ void rel(Home home
Integer view for integer variables.
SteelMill(const SteelMillOptions &opt)
Actual model.
virtual void archive(Archive &e) const
Archive into e.
struct Gecode::@585::NNF::@62::@64 a
For atomic nodes.
Choice for performing commit
virtual size_t dispose(Space &)
Delete brancher and return its size.
SteelMillOptions for examples with size option and an additional optional file name parameter...
Use LDSB for symmetry breaking.
bool assigned(View x, int v)
Whether x is assigned to value v.
void symmetry(int v)
Set default symmetry value.
void rel(Home home, FloatVar x0, FloatRelType frt, FloatVal n)
Propagates .
IntVarBranch INT_VAR_MAX_MIN(BranchTbl tbl)
Select variable with smallest max.
unsigned int size(void) const
Return size.
void solutions(unsigned int n)
Set default number of solutions to search for.
Post propagator for SetVar x
int ncolors(void) const
Return number of colors.
int * capacities(void) const
Return capacities.
Gecode toplevel namespace
SteelMillBranch(Space &home, SteelMillBranch &b)
Copy constructor.
Home class for posting propagators
unsigned int csplib_ncolors
Number of colors.
IntVarArray slabload
Load of slab j.
unsigned int csplib_ncapacities
Number of capacities.
Example: Steel-mill slab design problem
int main(int argc, char *argv[])
Main-function.
bool me_failed(ModEvent me)
Check whether modification event me is failed.
SteelMill(SteelMill &s)
Constructor for cloning s.
void element(Home home, IntSharedArray c, IntVar x0, IntVar x1, IntPropLevel)
Post domain consistent propagator for .
virtual void print(const Space &, const Gecode::Choice &_c, unsigned int a, std::ostream &o) const
Print explanation.
Custom brancher for steel mill slab design.
virtual void help(void)
Print help text.