COIN-OR::LEMON - Graph Library

Changeset 581:aa1804409f29 in lemon-1.2


Ignore:
Timestamp:
04/14/09 10:35:38 (15 years ago)
Author:
Peter Kovacs <kpeter@…>
Branch:
default
Phase:
public
Message:

Exploit that the standard maps are reference maps (#190)

Files:
9 edited

Legend:

Unmodified
Added
Removed
  • lemon/circulation.h

    r559 r581  
    454454
    455455      for(NodeIt n(_g);n!=INVALID;++n) {
    456         _excess->set(n, (*_delta)[n]);
     456        (*_excess)[n] = (*_delta)[n];
    457457      }
    458458
    459459      for (ArcIt e(_g);e!=INVALID;++e) {
    460460        _flow->set(e, (*_lo)[e]);
    461         _excess->set(_g.target(e), (*_excess)[_g.target(e)] + (*_flow)[e]);
    462         _excess->set(_g.source(e), (*_excess)[_g.source(e)] - (*_flow)[e]);
     461        (*_excess)[_g.target(e)] += (*_flow)[e];
     462        (*_excess)[_g.source(e)] -= (*_flow)[e];
    463463      }
    464464
     
    483483
    484484      for(NodeIt n(_g);n!=INVALID;++n) {
    485         _excess->set(n, (*_delta)[n]);
     485        (*_excess)[n] = (*_delta)[n];
    486486      }
    487487
     
    489489        if (!_tol.positive((*_excess)[_g.target(e)] + (*_up)[e])) {
    490490          _flow->set(e, (*_up)[e]);
    491           _excess->set(_g.target(e), (*_excess)[_g.target(e)] + (*_up)[e]);
    492           _excess->set(_g.source(e), (*_excess)[_g.source(e)] - (*_up)[e]);
     491          (*_excess)[_g.target(e)] += (*_up)[e];
     492          (*_excess)[_g.source(e)] -= (*_up)[e];
    493493        } else if (_tol.positive((*_excess)[_g.target(e)] + (*_lo)[e])) {
    494494          _flow->set(e, (*_lo)[e]);
    495           _excess->set(_g.target(e), (*_excess)[_g.target(e)] + (*_lo)[e]);
    496           _excess->set(_g.source(e), (*_excess)[_g.source(e)] - (*_lo)[e]);
     495          (*_excess)[_g.target(e)] += (*_lo)[e];
     496          (*_excess)[_g.source(e)] -= (*_lo)[e];
    497497        } else {
    498498          Value fc = -(*_excess)[_g.target(e)];
    499499          _flow->set(e, fc);
    500           _excess->set(_g.target(e), 0);
    501           _excess->set(_g.source(e), (*_excess)[_g.source(e)] - fc);
     500          (*_excess)[_g.target(e)] = 0;
     501          (*_excess)[_g.source(e)] -= fc;
    502502        }
    503503      }
     
    538538            if(!_tol.less(fc, exc)) {
    539539              _flow->set(e, (*_flow)[e] + exc);
    540               _excess->set(v, (*_excess)[v] + exc);
     540              (*_excess)[v] += exc;
    541541              if(!_level->active(v) && _tol.positive((*_excess)[v]))
    542542                _level->activate(v);
    543               _excess->set(act,0);
     543              (*_excess)[act] = 0;
    544544              _level->deactivate(act);
    545545              goto next_l;
     
    547547            else {
    548548              _flow->set(e, (*_up)[e]);
    549               _excess->set(v, (*_excess)[v] + fc);
     549              (*_excess)[v] += fc;
    550550              if(!_level->active(v) && _tol.positive((*_excess)[v]))
    551551                _level->activate(v);
     
    562562            if(!_tol.less(fc, exc)) {
    563563              _flow->set(e, (*_flow)[e] - exc);
    564               _excess->set(v, (*_excess)[v] + exc);
     564              (*_excess)[v] += exc;
    565565              if(!_level->active(v) && _tol.positive((*_excess)[v]))
    566566                _level->activate(v);
    567               _excess->set(act,0);
     567              (*_excess)[act] = 0;
    568568              _level->deactivate(act);
    569569              goto next_l;
     
    571571            else {
    572572              _flow->set(e, (*_lo)[e]);
    573               _excess->set(v, (*_excess)[v] + fc);
     573              (*_excess)[v] += fc;
    574574              if(!_level->active(v) && _tol.positive((*_excess)[v]))
    575575                _level->activate(v);
     
    580580        }
    581581
    582         _excess->set(act, exc);
     582        (*_excess)[act] = exc;
    583583        if(!_tol.positive(exc)) _level->deactivate(act);
    584584        else if(mlevel==_node_num) {
  • lemon/core.h

    r559 r581  
    13161316    virtual void clear() {
    13171317      for(NodeIt n(_g);n!=INVALID;++n) {
    1318         _head.set(n, INVALID);
     1318        _head[n] = INVALID;
    13191319      }
    13201320    }
     
    13231323      Node s = _g.source(arc);
    13241324      Node t = _g.target(arc);
    1325       _left.set(arc, INVALID);
    1326       _right.set(arc, INVALID);
     1325      _left[arc] = INVALID;
     1326      _right[arc] = INVALID;
    13271327
    13281328      Arc e = _head[s];
    13291329      if (e == INVALID) {
    1330         _head.set(s, arc);
    1331         _parent.set(arc, INVALID);
     1330        _head[s] = arc;
     1331        _parent[arc] = INVALID;
    13321332        return;
    13331333      }
     
    13351335        if (t < _g.target(e)) {
    13361336          if (_left[e] == INVALID) {
    1337             _left.set(e, arc);
    1338             _parent.set(arc, e);
     1337            _left[e] = arc;
     1338            _parent[arc] = e;
    13391339            splay(arc);
    13401340            return;
     
    13441344        } else {
    13451345          if (_right[e] == INVALID) {
    1346             _right.set(e, arc);
    1347             _parent.set(arc, e);
     1346            _right[e] = arc;
     1347            _parent[arc] = e;
    13481348            splay(arc);
    13491349            return;
     
    13581358      if (_left[arc] == INVALID) {
    13591359        if (_right[arc] != INVALID) {
    1360           _parent.set(_right[arc], _parent[arc]);
     1360          _parent[_right[arc]] = _parent[arc];
    13611361        }
    13621362        if (_parent[arc] != INVALID) {
    13631363          if (_left[_parent[arc]] == arc) {
    1364             _left.set(_parent[arc], _right[arc]);
     1364            _left[_parent[arc]] = _right[arc];
    13651365          } else {
    1366             _right.set(_parent[arc], _right[arc]);
     1366            _right[_parent[arc]] = _right[arc];
    13671367          }
    13681368        } else {
    1369           _head.set(_g.source(arc), _right[arc]);
     1369          _head[_g.source(arc)] = _right[arc];
    13701370        }
    13711371      } else if (_right[arc] == INVALID) {
    1372         _parent.set(_left[arc], _parent[arc]);
     1372        _parent[_left[arc]] = _parent[arc];
    13731373        if (_parent[arc] != INVALID) {
    13741374          if (_left[_parent[arc]] == arc) {
    1375             _left.set(_parent[arc], _left[arc]);
     1375            _left[_parent[arc]] = _left[arc];
    13761376          } else {
    1377             _right.set(_parent[arc], _left[arc]);
     1377            _right[_parent[arc]] = _left[arc];
    13781378          }
    13791379        } else {
    1380           _head.set(_g.source(arc), _left[arc]);
     1380          _head[_g.source(arc)] = _left[arc];
    13811381        }
    13821382      } else {
     
    13881388          }
    13891389          Arc s = _parent[e];
    1390           _right.set(_parent[e], _left[e]);
     1390          _right[_parent[e]] = _left[e];
    13911391          if (_left[e] != INVALID) {
    1392             _parent.set(_left[e], _parent[e]);
     1392            _parent[_left[e]] = _parent[e];
    13931393          }
    13941394
    1395           _left.set(e, _left[arc]);
    1396           _parent.set(_left[arc], e);
    1397           _right.set(e, _right[arc]);
    1398           _parent.set(_right[arc], e);
    1399 
    1400           _parent.set(e, _parent[arc]);
     1395          _left[e] = _left[arc];
     1396          _parent[_left[arc]] = e;
     1397          _right[e] = _right[arc];
     1398          _parent[_right[arc]] = e;
     1399
     1400          _parent[e] = _parent[arc];
    14011401          if (_parent[arc] != INVALID) {
    14021402            if (_left[_parent[arc]] == arc) {
    1403               _left.set(_parent[arc], e);
     1403              _left[_parent[arc]] = e;
    14041404            } else {
    1405               _right.set(_parent[arc], e);
     1405              _right[_parent[arc]] = e;
    14061406            }
    14071407          }
    14081408          splay(s);
    14091409        } else {
    1410           _right.set(e, _right[arc]);
    1411           _parent.set(_right[arc], e);
    1412           _parent.set(e, _parent[arc]);
     1410          _right[e] = _right[arc];
     1411          _parent[_right[arc]] = e;
     1412          _parent[e] = _parent[arc];
    14131413
    14141414          if (_parent[arc] != INVALID) {
    14151415            if (_left[_parent[arc]] == arc) {
    1416               _left.set(_parent[arc], e);
     1416              _left[_parent[arc]] = e;
    14171417            } else {
    1418               _right.set(_parent[arc], e);
     1418              _right[_parent[arc]] = e;
    14191419            }
    14201420          } else {
    1421             _head.set(_g.source(arc), e);
     1421            _head[_g.source(arc)] = e;
    14221422          }
    14231423        }
     
    14311431      if (a < m) {
    14321432        Arc left = refreshRec(v,a,m-1);
    1433         _left.set(me, left);
    1434         _parent.set(left, me);
     1433        _left[me] = left;
     1434        _parent[left] = me;
    14351435      } else {
    1436         _left.set(me, INVALID);
     1436        _left[me] = INVALID;
    14371437      }
    14381438      if (m < b) {
    14391439        Arc right = refreshRec(v,m+1,b);
    1440         _right.set(me, right);
    1441         _parent.set(right, me);
     1440        _right[me] = right;
     1441        _parent[right] = me;
    14421442      } else {
    1443         _right.set(me, INVALID);
     1443        _right[me] = INVALID;
    14441444      }
    14451445      return me;
     
    14531453          std::sort(v.begin(),v.end(),ArcLess(_g));
    14541454          Arc head = refreshRec(v,0,v.size()-1);
    1455           _head.set(n, head);
    1456           _parent.set(head, INVALID);
    1457         }
    1458         else _head.set(n, INVALID);
     1455          _head[n] = head;
     1456          _parent[head] = INVALID;
     1457        }
     1458        else _head[n] = INVALID;
    14591459      }
    14601460    }
     
    14621462    void zig(Arc v) {
    14631463      Arc w = _parent[v];
    1464       _parent.set(v, _parent[w]);
    1465       _parent.set(w, v);
    1466       _left.set(w, _right[v]);
    1467       _right.set(v, w);
     1464      _parent[v] = _parent[w];
     1465      _parent[w] = v;
     1466      _left[w] = _right[v];
     1467      _right[v] = w;
    14681468      if (_parent[v] != INVALID) {
    14691469        if (_right[_parent[v]] == w) {
    1470           _right.set(_parent[v], v);
     1470          _right[_parent[v]] = v;
    14711471        } else {
    1472           _left.set(_parent[v], v);
     1472          _left[_parent[v]] = v;
    14731473        }
    14741474      }
    14751475      if (_left[w] != INVALID){
    1476         _parent.set(_left[w], w);
     1476        _parent[_left[w]] = w;
    14771477      }
    14781478    }
     
    14801480    void zag(Arc v) {
    14811481      Arc w = _parent[v];
    1482       _parent.set(v, _parent[w]);
    1483       _parent.set(w, v);
    1484       _right.set(w, _left[v]);
    1485       _left.set(v, w);
     1482      _parent[v] = _parent[w];
     1483      _parent[w] = v;
     1484      _right[w] = _left[v];
     1485      _left[v] = w;
    14861486      if (_parent[v] != INVALID){
    14871487        if (_left[_parent[v]] == w) {
    1488           _left.set(_parent[v], v);
     1488          _left[_parent[v]] = v;
    14891489        } else {
    1490           _right.set(_parent[v], v);
     1490          _right[_parent[v]] = v;
    14911491        }
    14921492      }
    14931493      if (_right[w] != INVALID){
    1494         _parent.set(_right[w], w);
     1494        _parent[_right[w]] = w;
    14951495      }
    14961496    }
  • lemon/elevator.h

    r559 r581  
    7777    void copy(Item i, Vit p)
    7878    {
    79       _where.set(*p=i,p);
     79      _where[*p=i] = p;
    8080    }
    8181    void copy(Vit s, Vit p)
     
    8585          Item i=*s;
    8686          *p=i;
    87           _where.set(i,p);
     87          _where[i] = p;
    8888        }
    8989    }
     
    9292      Item ti=*i;
    9393      Vit ct = _where[ti];
    94       _where.set(ti,_where[*i=*j]);
    95       _where.set(*j,ct);
     94      _where[ti] = _where[*i=*j];
     95      _where[*j] = ct;
    9696      *j=ti;
    9797    }
     
    227227    {
    228228      Item it = *_last_active[_highest_active];
    229       _level.set(it,_level[it]+1);
     229      ++_level[it];
    230230      swap(_last_active[_highest_active]--,_last_active[_highest_active+1]);
    231231      --_first[++_highest_active];
     
    250250        }
    251251      copy(li,_first[new_level]);
    252       _level.set(li,new_level);
     252      _level[li] = new_level;
    253253      _highest_active=new_level;
    254254    }
     
    270270      copy(li,_first[_max_level]);
    271271      --_last_active[_max_level];
    272       _level.set(li,_max_level);
     272      _level[li] = _max_level;
    273273
    274274      while(_highest_active>=0 &&
     
    300300    {
    301301      Item it =*_last_active[level];
    302       _level.set(it,_level[it]+1);
     302      ++_level[it];
    303303      swap(_last_active[level]--, --_first[level+1]);
    304304      if (level+1>_highest_active) ++_highest_active;
     
    320320        }
    321321      copy(ai,_first[new_level]);
    322       _level.set(ai,new_level);
     322      _level[ai] = new_level;
    323323      if (new_level>_highest_active) _highest_active=new_level;
    324324    }
     
    340340      copy(ai,_first[_max_level]);
    341341      --_last_active[_max_level];
    342       _level.set(ai,_max_level);
     342      _level[ai] = _max_level;
    343343
    344344      if (_highest_active==level) {
     
    371371        }
    372372      copy(i,_first[new_level]);
    373       _level.set(i,new_level);
     373      _level[i] = new_level;
    374374      if(new_level>_highest_active) _highest_active=new_level;
    375375    }
     
    383383    ///\pre The item is on the top level.
    384384    void dirtyTopButOne(Item i) {
    385       _level.set(i,_max_level - 1);
     385      _level[i] = _max_level - 1;
    386386    }
    387387
     
    395395      const Vit tl=_first[_max_level];
    396396      for(Vit i=f;i!=tl;++i)
    397         _level.set(*i,_max_level);
     397        _level[*i] = _max_level;
    398398      for(int i=l;i<=_max_level;i++)
    399399        {
     
    434434        {
    435435          *n=i;
    436           _where.set(i,n);
    437           _level.set(i,_max_level);
     436          _where[i] = n;
     437          _level[i] = _max_level;
    438438          ++n;
    439439        }
     
    444444    {
    445445      swap(_where[i],_init_num);
    446       _level.set(i,_init_lev);
     446      _level[i] = _init_lev;
    447447      ++_init_num;
    448448    }
     
    552552    ///\pre Item \c i shouldn't be active before.
    553553    void activate(Item i) {
    554       _active.set(i, true);
     554      _active[i] = true;
    555555
    556556      int level = _level[i];
     
    561561      if (_prev[i] == INVALID || _active[_prev[i]]) return;
    562562      //unlace
    563       _next.set(_prev[i], _next[i]);
     563      _next[_prev[i]] = _next[i];
    564564      if (_next[i] != INVALID) {
    565         _prev.set(_next[i], _prev[i]);
     565        _prev[_next[i]] = _prev[i];
    566566      } else {
    567567        _last[level] = _prev[i];
    568568      }
    569569      //lace
    570       _next.set(i, _first[level]);
    571       _prev.set(_first[level], i);
    572       _prev.set(i, INVALID);
     570      _next[i] = _first[level];
     571      _prev[_first[level]] = i;
     572      _prev[i] = INVALID;
    573573      _first[level] = i;
    574574
     
    580580    ///\pre Item \c i must be active before.
    581581    void deactivate(Item i) {
    582       _active.set(i, false);
     582      _active[i] = false;
    583583      int level = _level[i];
    584584
     
    587587
    588588      //unlace
    589       _prev.set(_next[i], _prev[i]);
     589      _prev[_next[i]] = _prev[i];
    590590      if (_prev[i] != INVALID) {
    591         _next.set(_prev[i], _next[i]);
     591        _next[_prev[i]] = _next[i];
    592592      } else {
    593593        _first[_level[i]] = _next[i];
    594594      }
    595595      //lace
    596       _prev.set(i, _last[level]);
    597       _next.set(_last[level], i);
    598       _next.set(i, INVALID);
     596      _prev[i] = _last[level];
     597      _next[_last[level]] = i;
     598      _next[i] = INVALID;
    599599      _last[level] = i;
    600600
     
    686686      Item i = _first[_highest_active];
    687687      if (_next[i] != INVALID) {
    688         _prev.set(_next[i], INVALID);
     688        _prev[_next[i]] = INVALID;
    689689        _first[_highest_active] = _next[i];
    690690      } else {
     
    692692        _last[_highest_active] = INVALID;
    693693      }
    694       _level.set(i, ++_highest_active);
     694      _level[i] = ++_highest_active;
    695695      if (_first[_highest_active] == INVALID) {
    696696        _first[_highest_active] = i;
    697697        _last[_highest_active] = i;
    698         _prev.set(i, INVALID);
    699         _next.set(i, INVALID);
    700       } else {
    701         _prev.set(_first[_highest_active], i);
    702         _next.set(i, _first[_highest_active]);
     698        _prev[i] = INVALID;
     699        _next[i] = INVALID;
     700      } else {
     701        _prev[_first[_highest_active]] = i;
     702        _next[i] = _first[_highest_active];
    703703        _first[_highest_active] = i;
    704704      }
     
    715715      Item i = _first[_highest_active];
    716716      if (_next[i] != INVALID) {
    717         _prev.set(_next[i], INVALID);
     717        _prev[_next[i]] = INVALID;
    718718        _first[_highest_active] = _next[i];
    719719      } else {
     
    721721        _last[_highest_active] = INVALID;
    722722      }
    723       _level.set(i, _highest_active = new_level);
     723      _level[i] = _highest_active = new_level;
    724724      if (_first[_highest_active] == INVALID) {
    725725        _first[_highest_active] = _last[_highest_active] = i;
    726         _prev.set(i, INVALID);
    727         _next.set(i, INVALID);
    728       } else {
    729         _prev.set(_first[_highest_active], i);
    730         _next.set(i, _first[_highest_active]);
     726        _prev[i] = INVALID;
     727        _next[i] = INVALID;
     728      } else {
     729        _prev[_first[_highest_active]] = i;
     730        _next[i] = _first[_highest_active];
    731731        _first[_highest_active] = i;
    732732      }
     
    739739    void liftHighestActiveToTop() {
    740740      Item i = _first[_highest_active];
    741       _level.set(i, _max_level);
     741      _level[i] = _max_level;
    742742      if (_next[i] != INVALID) {
    743         _prev.set(_next[i], INVALID);
     743        _prev[_next[i]] = INVALID;
    744744        _first[_highest_active] = _next[i];
    745745      } else {
     
    775775      Item i = _first[l];
    776776      if (_next[i] != INVALID) {
    777         _prev.set(_next[i], INVALID);
     777        _prev[_next[i]] = INVALID;
    778778        _first[l] = _next[i];
    779779      } else {
     
    781781        _last[l] = INVALID;
    782782      }
    783       _level.set(i, ++l);
     783      _level[i] = ++l;
    784784      if (_first[l] == INVALID) {
    785785        _first[l] = _last[l] = i;
    786         _prev.set(i, INVALID);
    787         _next.set(i, INVALID);
    788       } else {
    789         _prev.set(_first[l], i);
    790         _next.set(i, _first[l]);
     786        _prev[i] = INVALID;
     787        _next[i] = INVALID;
     788      } else {
     789        _prev[_first[l]] = i;
     790        _next[i] = _first[l];
    791791        _first[l] = i;
    792792      }
     
    804804      Item i = _first[l];
    805805      if (_next[i] != INVALID) {
    806         _prev.set(_next[i], INVALID);
     806        _prev[_next[i]] = INVALID;
    807807        _first[l] = _next[i];
    808808      } else {
     
    810810        _last[l] = INVALID;
    811811      }
    812       _level.set(i, l = new_level);
     812      _level[i] = l = new_level;
    813813      if (_first[l] == INVALID) {
    814814        _first[l] = _last[l] = i;
    815         _prev.set(i, INVALID);
    816         _next.set(i, INVALID);
    817       } else {
    818         _prev.set(_first[l], i);
    819         _next.set(i, _first[l]);
     815        _prev[i] = INVALID;
     816        _next[i] = INVALID;
     817      } else {
     818        _prev[_first[l]] = i;
     819        _next[i] = _first[l];
    820820        _first[l] = i;
    821821      }
     
    833833      Item i = _first[l];
    834834      if (_next[i] != INVALID) {
    835         _prev.set(_next[i], INVALID);
     835        _prev[_next[i]] = INVALID;
    836836        _first[l] = _next[i];
    837837      } else {
     
    839839        _last[l] = INVALID;
    840840      }
    841       _level.set(i, _max_level);
     841      _level[i] = _max_level;
    842842      if (l == _highest_active) {
    843843        while (_highest_active >= 0 && activeFree(_highest_active))
     
    857857    void lift(Item i, int new_level) {
    858858      if (_next[i] != INVALID) {
    859         _prev.set(_next[i], _prev[i]);
     859        _prev[_next[i]] = _prev[i];
    860860      } else {
    861861        _last[new_level] = _prev[i];
    862862      }
    863863      if (_prev[i] != INVALID) {
    864         _next.set(_prev[i], _next[i]);
     864        _next[_prev[i]] = _next[i];
    865865      } else {
    866866        _first[new_level] = _next[i];
    867867      }
    868       _level.set(i, new_level);
     868      _level[i] = new_level;
    869869      if (_first[new_level] == INVALID) {
    870870        _first[new_level] = _last[new_level] = i;
    871         _prev.set(i, INVALID);
    872         _next.set(i, INVALID);
    873       } else {
    874         _prev.set(_first[new_level], i);
    875         _next.set(i, _first[new_level]);
     871        _prev[i] = INVALID;
     872        _next[i] = INVALID;
     873      } else {
     874        _prev[_first[new_level]] = i;
     875        _next[i] = _first[new_level];
    876876        _first[new_level] = i;
    877877      }
     
    889889    ///\pre The item is on the top level.
    890890    void dirtyTopButOne(Item i) {
    891       _level.set(i, _max_level - 1);
     891      _level[i] = _max_level - 1;
    892892    }
    893893
     
    900900        Item n = _first[i];
    901901        while (n != INVALID) {
    902           _level.set(n, _max_level);
     902          _level[n] = _max_level;
    903903          n = _next[n];
    904904        }
     
    938938      for(typename ItemSetTraits<GR,Item>::ItemIt i(_graph);
    939939          i != INVALID; ++i) {
    940         _level.set(i, _max_level);
    941         _active.set(i, false);
     940        _level[i] = _max_level;
     941        _active[i] = false;
    942942      }
    943943    }
     
    945945    ///Add an item to the current level.
    946946    void initAddItem(Item i) {
    947       _level.set(i, _init_level);
     947      _level[i] = _init_level;
    948948      if (_last[_init_level] == INVALID) {
    949949        _first[_init_level] = i;
    950950        _last[_init_level] = i;
    951         _prev.set(i, INVALID);
    952         _next.set(i, INVALID);
    953       } else {
    954         _prev.set(i, _last[_init_level]);
    955         _next.set(i, INVALID);
    956         _next.set(_last[_init_level], i);
     951        _prev[i] = INVALID;
     952        _next[i] = INVALID;
     953      } else {
     954        _prev[i] = _last[_init_level];
     955        _next[i] = INVALID;
     956        _next[_last[_init_level]] = i;
    957957        _last[_init_level] = i;
    958958      }
  • lemon/gomory_hu.h

    r546 r581  
    144144      _root = NodeIt(_graph);
    145145      for (NodeIt n(_graph); n != INVALID; ++n) {
    146         _pred->set(n, _root);
    147         _order->set(n, -1);
    148       }
    149       _pred->set(_root, INVALID);
    150       _weight->set(_root, std::numeric_limits<Value>::max());
     146        (*_pred)[n] = _root;
     147        (*_order)[n] = -1;
     148      }
     149      (*_pred)[_root] = INVALID;
     150      (*_weight)[_root] = std::numeric_limits<Value>::max();
    151151    }
    152152
     
    165165        fa.runMinCut();
    166166
    167         _weight->set(n, fa.flowValue());
     167        (*_weight)[n] = fa.flowValue();
    168168
    169169        for (NodeIt nn(_graph); nn != INVALID; ++nn) {
    170170          if (nn != n && fa.minCut(nn) && (*_pred)[nn] == pn) {
    171             _pred->set(nn, n);
     171            (*_pred)[nn] = n;
    172172          }
    173173        }
    174174        if ((*_pred)[pn] != INVALID && fa.minCut((*_pred)[pn])) {
    175           _pred->set(n, (*_pred)[pn]);
    176           _pred->set(pn, n);
    177           _weight->set(n, (*_weight)[pn]);
    178           _weight->set(pn, fa.flowValue());     
    179         }
    180       }
    181 
    182       _order->set(_root, 0);
     175          (*_pred)[n] = (*_pred)[pn];
     176          (*_pred)[pn] = n;
     177          (*_weight)[n] = (*_weight)[pn];
     178          (*_weight)[pn] = fa.flowValue();
     179        }
     180      }
     181
     182      (*_order)[_root] = 0;
    183183      int index = 1;
    184184
     
    191191        }
    192192        while (!st.empty()) {
    193           _order->set(st.back(), index++);
     193          (*_order)[st.back()] = index++;
    194194          st.pop_back();
    195195        }
     
    310310
    311311      typename Graph::template NodeMap<bool> reached(_graph, false);
    312       reached.set(_root, true);
     312      reached[_root] = true;
    313313      cutMap.set(_root, !s_root);
    314       reached.set(rn, true);
     314      reached[rn] = true;
    315315      cutMap.set(rn, s_root);
    316316
  • lemon/hao_orlin.h

    r559 r581  
    162162
    163163    void activate(const Node& i) {
    164       _active->set(i, true);
     164      (*_active)[i] = true;
    165165
    166166      int bucket = (*_bucket)[i];
     
    168168      if ((*_prev)[i] == INVALID || (*_active)[(*_prev)[i]]) return;
    169169      //unlace
    170       _next->set((*_prev)[i], (*_next)[i]);
     170      (*_next)[(*_prev)[i]] = (*_next)[i];
    171171      if ((*_next)[i] != INVALID) {
    172         _prev->set((*_next)[i], (*_prev)[i]);
     172        (*_prev)[(*_next)[i]] = (*_prev)[i];
    173173      } else {
    174174        _last[bucket] = (*_prev)[i];
    175175      }
    176176      //lace
    177       _next->set(i, _first[bucket]);
    178       _prev->set(_first[bucket], i);
    179       _prev->set(i, INVALID);
     177      (*_next)[i] = _first[bucket];
     178      (*_prev)[_first[bucket]] = i;
     179      (*_prev)[i] = INVALID;
    180180      _first[bucket] = i;
    181181    }
    182182
    183183    void deactivate(const Node& i) {
    184       _active->set(i, false);
     184      (*_active)[i] = false;
    185185      int bucket = (*_bucket)[i];
    186186
     
    188188
    189189      //unlace
    190       _prev->set((*_next)[i], (*_prev)[i]);
     190      (*_prev)[(*_next)[i]] = (*_prev)[i];
    191191      if ((*_prev)[i] != INVALID) {
    192         _next->set((*_prev)[i], (*_next)[i]);
     192        (*_next)[(*_prev)[i]] = (*_next)[i];
    193193      } else {
    194194        _first[bucket] = (*_next)[i];
    195195      }
    196196      //lace
    197       _prev->set(i, _last[bucket]);
    198       _next->set(_last[bucket], i);
    199       _next->set(i, INVALID);
     197      (*_prev)[i] = _last[bucket];
     198      (*_next)[_last[bucket]] = i;
     199      (*_next)[i] = INVALID;
    200200      _last[bucket] = i;
    201201    }
     
    204204      (*_bucket)[i] = bucket;
    205205      if (_last[bucket] != INVALID) {
    206         _prev->set(i, _last[bucket]);
    207         _next->set(_last[bucket], i);
    208         _next->set(i, INVALID);
     206        (*_prev)[i] = _last[bucket];
     207        (*_next)[_last[bucket]] = i;
     208        (*_next)[i] = INVALID;
    209209        _last[bucket] = i;
    210210      } else {
    211         _prev->set(i, INVALID);
     211        (*_prev)[i] = INVALID;
    212212        _first[bucket] = i;
    213         _next->set(i, INVALID);
     213        (*_next)[i] = INVALID;
    214214        _last[bucket] = i;
    215215      }
     
    219219
    220220      for (NodeIt n(_graph); n != INVALID; ++n) {
    221         _excess->set(n, 0);
     221        (*_excess)[n] = 0;
    222222      }
    223223
    224224      for (ArcIt a(_graph); a != INVALID; ++a) {
    225         _flow->set(a, 0);
     225        (*_flow)[a] = 0;
    226226      }
    227227
     
    233233        typename Digraph::template NodeMap<bool> reached(_graph, false);
    234234
    235         reached.set(_source, true);
     235        reached[_source] = true;
    236236        bool first_set = true;
    237237
     
    241241
    242242          queue[qlast++] = t;
    243           reached.set(t, true);
     243          reached[t] = true;
    244244
    245245          while (qfirst != qlast) {
     
    258258              Node u = _graph.source(a);
    259259              if (!reached[u] && _tolerance.positive((*_capacity)[a])) {
    260                 reached.set(u, true);
     260                reached[u] = true;
    261261                queue[qlast++] = u;
    262262              }
     
    267267
    268268        ++bucket_num;
    269         _bucket->set(_source, 0);
     269        (*_bucket)[_source] = 0;
    270270        _dormant[0] = true;
    271271      }
    272       _source_set->set(_source, true);
     272      (*_source_set)[_source] = true;
    273273
    274274      Node target = _last[_sets.back().back()];
     
    277277          if (_tolerance.positive((*_capacity)[a])) {
    278278            Node u = _graph.target(a);
    279             _flow->set(a, (*_capacity)[a]);
    280             _excess->set(u, (*_excess)[u] + (*_capacity)[a]);
     279            (*_flow)[a] = (*_capacity)[a];
     280            (*_excess)[u] += (*_capacity)[a];
    281281            if (!(*_active)[u] && u != _source) {
    282282              activate(u);
     
    319319              }
    320320              if (!_tolerance.less(rem, excess)) {
    321                 _flow->set(a, (*_flow)[a] + excess);
    322                 _excess->set(v, (*_excess)[v] + excess);
     321                (*_flow)[a] += excess;
     322                (*_excess)[v] += excess;
    323323                excess = 0;
    324324                goto no_more_push;
    325325              } else {
    326326                excess -= rem;
    327                 _excess->set(v, (*_excess)[v] + rem);
    328                 _flow->set(a, (*_capacity)[a]);
     327                (*_excess)[v] += rem;
     328                (*_flow)[a] = (*_capacity)[a];
    329329              }
    330330            } else if (next_bucket > (*_bucket)[v]) {
     
    343343              }
    344344              if (!_tolerance.less(rem, excess)) {
    345                 _flow->set(a, (*_flow)[a] - excess);
    346                 _excess->set(v, (*_excess)[v] + excess);
     345                (*_flow)[a] -= excess;
     346                (*_excess)[v] += excess;
    347347                excess = 0;
    348348                goto no_more_push;
    349349              } else {
    350350                excess -= rem;
    351                 _excess->set(v, (*_excess)[v] + rem);
    352                 _flow->set(a, 0);
     351                (*_excess)[v] += rem;
     352                (*_flow)[a] = 0;
    353353              }
    354354            } else if (next_bucket > (*_bucket)[v]) {
     
    359359        no_more_push:
    360360
    361           _excess->set(n, excess);
     361          (*_excess)[n] = excess;
    362362
    363363          if (excess != 0) {
     
    377377            } else if (next_bucket == _node_num) {
    378378              _first[(*_bucket)[n]] = (*_next)[n];
    379               _prev->set((*_next)[n], INVALID);
     379              (*_prev)[(*_next)[n]] = INVALID;
    380380
    381381              std::list<std::list<int> >::iterator new_set =
     
    383383
    384384              new_set->push_front(bucket_num);
    385               _bucket->set(n, bucket_num);
     385              (*_bucket)[n] = bucket_num;
    386386              _first[bucket_num] = _last[bucket_num] = n;
    387               _next->set(n, INVALID);
    388               _prev->set(n, INVALID);
     387              (*_next)[n] = INVALID;
     388              (*_prev)[n] = INVALID;
    389389              _dormant[bucket_num] = true;
    390390              ++bucket_num;
     
    396396            } else {
    397397              _first[*_highest] = (*_next)[n];
    398               _prev->set((*_next)[n], INVALID);
     398              (*_prev)[(*_next)[n]] = INVALID;
    399399
    400400              while (next_bucket != *_highest) {
     
    410410              --_highest;
    411411
    412               _bucket->set(n, *_highest);
    413               _next->set(n, _first[*_highest]);
     412              (*_bucket)[n] = *_highest;
     413              (*_next)[n] = _first[*_highest];
    414414              if (_first[*_highest] != INVALID) {
    415                 _prev->set(_first[*_highest], n);
     415                (*_prev)[_first[*_highest]] = n;
    416416              } else {
    417417                _last[*_highest] = n;
     
    435435          _min_cut = (*_excess)[target];
    436436          for (NodeIt i(_graph); i != INVALID; ++i) {
    437             _min_cut_map->set(i, true);
     437            (*_min_cut_map)[i] = true;
    438438          }
    439439          for (std::list<int>::iterator it = _sets.back().begin();
     
    441441            Node n = _first[*it];
    442442            while (n != INVALID) {
    443               _min_cut_map->set(n, false);
     443              (*_min_cut_map)[n] = false;
    444444              n = (*_next)[n];
    445445            }
     
    454454              new_target = (*_prev)[target];
    455455            } else {
    456               _prev->set((*_next)[target], (*_prev)[target]);
     456              (*_prev)[(*_next)[target]] = (*_prev)[target];
    457457              new_target = (*_next)[target];
    458458            }
     
    460460              _first[(*_bucket)[target]] = (*_next)[target];
    461461            } else {
    462               _next->set((*_prev)[target], (*_next)[target]);
     462              (*_next)[(*_prev)[target]] = (*_next)[target];
    463463            }
    464464          } else {
     
    476476          }
    477477
    478           _bucket->set(target, 0);
    479 
    480           _source_set->set(target, true);
     478          (*_bucket)[target] = 0;
     479
     480          (*_source_set)[target] = true;
    481481          for (OutArcIt a(_graph, target); a != INVALID; ++a) {
    482482            Value rem = (*_capacity)[a] - (*_flow)[a];
     
    486486              activate(v);
    487487            }
    488             _excess->set(v, (*_excess)[v] + rem);
    489             _flow->set(a, (*_capacity)[a]);
     488            (*_excess)[v] += rem;
     489            (*_flow)[a] = (*_capacity)[a];
    490490          }
    491491
     
    497497              activate(v);
    498498            }
    499             _excess->set(v, (*_excess)[v] + rem);
    500             _flow->set(a, 0);
     499            (*_excess)[v] += rem;
     500            (*_flow)[a] = 0;
    501501          }
    502502
     
    518518
    519519      for (NodeIt n(_graph); n != INVALID; ++n) {
    520         _excess->set(n, 0);
     520        (*_excess)[n] = 0;
    521521      }
    522522
    523523      for (ArcIt a(_graph); a != INVALID; ++a) {
    524         _flow->set(a, 0);
     524        (*_flow)[a] = 0;
    525525      }
    526526
     
    532532        typename Digraph::template NodeMap<bool> reached(_graph, false);
    533533
    534         reached.set(_source, true);
     534        reached[_source] = true;
    535535
    536536        bool first_set = true;
     
    541541
    542542          queue[qlast++] = t;
    543           reached.set(t, true);
     543          reached[t] = true;
    544544
    545545          while (qfirst != qlast) {
     
    558558              Node u = _graph.target(a);
    559559              if (!reached[u] && _tolerance.positive((*_capacity)[a])) {
    560                 reached.set(u, true);
     560                reached[u] = true;
    561561                queue[qlast++] = u;
    562562              }
     
    567567
    568568        ++bucket_num;
    569         _bucket->set(_source, 0);
     569        (*_bucket)[_source] = 0;
    570570        _dormant[0] = true;
    571571      }
    572       _source_set->set(_source, true);
     572      (*_source_set)[_source] = true;
    573573
    574574      Node target = _last[_sets.back().back()];
     
    577577          if (_tolerance.positive((*_capacity)[a])) {
    578578            Node u = _graph.source(a);
    579             _flow->set(a, (*_capacity)[a]);
    580             _excess->set(u, (*_excess)[u] + (*_capacity)[a]);
     579            (*_flow)[a] = (*_capacity)[a];
     580            (*_excess)[u] += (*_capacity)[a];
    581581            if (!(*_active)[u] && u != _source) {
    582582              activate(u);
     
    619619              }
    620620              if (!_tolerance.less(rem, excess)) {
    621                 _flow->set(a, (*_flow)[a] + excess);
    622                 _excess->set(v, (*_excess)[v] + excess);
     621                (*_flow)[a] += excess;
     622                (*_excess)[v] += excess;
    623623                excess = 0;
    624624                goto no_more_push;
    625625              } else {
    626626                excess -= rem;
    627                 _excess->set(v, (*_excess)[v] + rem);
    628                 _flow->set(a, (*_capacity)[a]);
     627                (*_excess)[v] += rem;
     628                (*_flow)[a] = (*_capacity)[a];
    629629              }
    630630            } else if (next_bucket > (*_bucket)[v]) {
     
    643643              }
    644644              if (!_tolerance.less(rem, excess)) {
    645                 _flow->set(a, (*_flow)[a] - excess);
    646                 _excess->set(v, (*_excess)[v] + excess);
     645                (*_flow)[a] -= excess;
     646                (*_excess)[v] += excess;
    647647                excess = 0;
    648648                goto no_more_push;
    649649              } else {
    650650                excess -= rem;
    651                 _excess->set(v, (*_excess)[v] + rem);
    652                 _flow->set(a, 0);
     651                (*_excess)[v] += rem;
     652                (*_flow)[a] = 0;
    653653              }
    654654            } else if (next_bucket > (*_bucket)[v]) {
     
    659659        no_more_push:
    660660
    661           _excess->set(n, excess);
     661          (*_excess)[n] = excess;
    662662
    663663          if (excess != 0) {
     
    677677            } else if (next_bucket == _node_num) {
    678678              _first[(*_bucket)[n]] = (*_next)[n];
    679               _prev->set((*_next)[n], INVALID);
     679              (*_prev)[(*_next)[n]] = INVALID;
    680680
    681681              std::list<std::list<int> >::iterator new_set =
     
    683683
    684684              new_set->push_front(bucket_num);
    685               _bucket->set(n, bucket_num);
     685              (*_bucket)[n] = bucket_num;
    686686              _first[bucket_num] = _last[bucket_num] = n;
    687               _next->set(n, INVALID);
    688               _prev->set(n, INVALID);
     687              (*_next)[n] = INVALID;
     688              (*_prev)[n] = INVALID;
    689689              _dormant[bucket_num] = true;
    690690              ++bucket_num;
     
    696696            } else {
    697697              _first[*_highest] = (*_next)[n];
    698               _prev->set((*_next)[n], INVALID);
     698              (*_prev)[(*_next)[n]] = INVALID;
    699699
    700700              while (next_bucket != *_highest) {
     
    709709              --_highest;
    710710
    711               _bucket->set(n, *_highest);
    712               _next->set(n, _first[*_highest]);
     711              (*_bucket)[n] = *_highest;
     712              (*_next)[n] = _first[*_highest];
    713713              if (_first[*_highest] != INVALID) {
    714                 _prev->set(_first[*_highest], n);
     714                (*_prev)[_first[*_highest]] = n;
    715715              } else {
    716716                _last[*_highest] = n;
     
    734734          _min_cut = (*_excess)[target];
    735735          for (NodeIt i(_graph); i != INVALID; ++i) {
    736             _min_cut_map->set(i, false);
     736            (*_min_cut_map)[i] = false;
    737737          }
    738738          for (std::list<int>::iterator it = _sets.back().begin();
     
    740740            Node n = _first[*it];
    741741            while (n != INVALID) {
    742               _min_cut_map->set(n, true);
     742              (*_min_cut_map)[n] = true;
    743743              n = (*_next)[n];
    744744            }
     
    753753              new_target = (*_prev)[target];
    754754            } else {
    755               _prev->set((*_next)[target], (*_prev)[target]);
     755              (*_prev)[(*_next)[target]] = (*_prev)[target];
    756756              new_target = (*_next)[target];
    757757            }
     
    759759              _first[(*_bucket)[target]] = (*_next)[target];
    760760            } else {
    761               _next->set((*_prev)[target], (*_next)[target]);
     761              (*_next)[(*_prev)[target]] = (*_next)[target];
    762762            }
    763763          } else {
     
    775775          }
    776776
    777           _bucket->set(target, 0);
    778 
    779           _source_set->set(target, true);
     777          (*_bucket)[target] = 0;
     778
     779          (*_source_set)[target] = true;
    780780          for (InArcIt a(_graph, target); a != INVALID; ++a) {
    781781            Value rem = (*_capacity)[a] - (*_flow)[a];
     
    785785              activate(v);
    786786            }
    787             _excess->set(v, (*_excess)[v] + rem);
    788             _flow->set(a, (*_capacity)[a]);
     787            (*_excess)[v] += rem;
     788            (*_flow)[a] = (*_capacity)[a];
    789789          }
    790790
     
    796796              activate(v);
    797797            }
    798             _excess->set(v, (*_excess)[v] + rem);
    799             _flow->set(a, 0);
     798            (*_excess)[v] += rem;
     799            (*_flow)[a] = 0;
    800800          }
    801801
  • lemon/max_matching.h

    r559 r581  
    283283
    284284        while (base != nca) {
    285           _ear->set(node, arc);
     285          (*_ear)[node] = arc;
    286286
    287287          Node n = node;
     
    290290            Arc a = (*_ear)[n];
    291291            n = _graph.target(a);
    292             _ear->set(n, _graph.oppositeArc(a));
     292            (*_ear)[n] = _graph.oppositeArc(a);
    293293          }
    294294          node = _graph.target((*_matching)[base]);
     
    296296          _tree_set->erase(node);
    297297          _blossom_set->insert(node, _blossom_set->find(base));
    298           _status->set(node, EVEN);
     298          (*_status)[node] = EVEN;
    299299          _node_queue[_last++] = node;
    300300          arc = _graph.oppositeArc((*_ear)[node]);
     
    305305      }
    306306
    307       _blossom_rep->set(_blossom_set->find(nca), nca);
     307      (*_blossom_rep)[_blossom_set->find(nca)] = nca;
    308308
    309309      {
     
    314314
    315315        while (base != nca) {
    316           _ear->set(node, arc);
     316          (*_ear)[node] = arc;
    317317
    318318          Node n = node;
     
    321321            Arc a = (*_ear)[n];
    322322            n = _graph.target(a);
    323             _ear->set(n, _graph.oppositeArc(a));
     323            (*_ear)[n] = _graph.oppositeArc(a);
    324324          }
    325325          node = _graph.target((*_matching)[base]);
     
    327327          _tree_set->erase(node);
    328328          _blossom_set->insert(node, _blossom_set->find(base));
    329           _status->set(node, EVEN);
     329          (*_status)[node] = EVEN;
    330330          _node_queue[_last++] = node;
    331331          arc = _graph.oppositeArc((*_ear)[node]);
     
    336336      }
    337337
    338       _blossom_rep->set(_blossom_set->find(nca), nca);
     338      (*_blossom_rep)[_blossom_set->find(nca)] = nca;
    339339    }
    340340
     
    345345      Node odd = _graph.target(a);
    346346
    347       _ear->set(odd, _graph.oppositeArc(a));
     347      (*_ear)[odd] = _graph.oppositeArc(a);
    348348      Node even = _graph.target((*_matching)[odd]);
    349       _blossom_rep->set(_blossom_set->insert(even), even);
    350       _status->set(odd, ODD);
    351       _status->set(even, EVEN);
     349      (*_blossom_rep)[_blossom_set->insert(even)] = even;
     350      (*_status)[odd] = ODD;
     351      (*_status)[even] = EVEN;
    352352      int tree = _tree_set->find((*_blossom_rep)[_blossom_set->find(base)]);
    353353      _tree_set->insert(odd, tree);
     
    363363      int tree = _tree_set->find((*_blossom_rep)[_blossom_set->find(even)]);
    364364
    365       _matching->set(odd, _graph.oppositeArc(a));
    366       _status->set(odd, MATCHED);
     365      (*_matching)[odd] = _graph.oppositeArc(a);
     366      (*_status)[odd] = MATCHED;
    367367
    368368      Arc arc = (*_matching)[even];
    369       _matching->set(even, a);
     369      (*_matching)[even] = a;
    370370
    371371      while (arc != INVALID) {
     
    373373        arc = (*_ear)[odd];
    374374        even = _graph.target(arc);
    375         _matching->set(odd, arc);
     375        (*_matching)[odd] = arc;
    376376        arc = (*_matching)[even];
    377         _matching->set(even, _graph.oppositeArc((*_matching)[odd]));
     377        (*_matching)[even] = _graph.oppositeArc((*_matching)[odd]);
    378378      }
    379379
     
    381381           it != INVALID; ++it) {
    382382        if ((*_status)[it] == ODD) {
    383           _status->set(it, MATCHED);
     383          (*_status)[it] = MATCHED;
    384384        } else {
    385385          int blossom = _blossom_set->find(it);
    386386          for (typename BlossomSet::ItemIt jt(*_blossom_set, blossom);
    387387               jt != INVALID; ++jt) {
    388             _status->set(jt, MATCHED);
     388            (*_status)[jt] = MATCHED;
    389389          }
    390390          _blossom_set->eraseClass(blossom);
     
    428428      createStructures();
    429429      for(NodeIt n(_graph); n != INVALID; ++n) {
    430         _matching->set(n, INVALID);
    431         _status->set(n, UNMATCHED);
     430        (*_matching)[n] = INVALID;
     431        (*_status)[n] = UNMATCHED;
    432432      }
    433433    }
     
    439439      createStructures();
    440440      for (NodeIt n(_graph); n != INVALID; ++n) {
    441         _matching->set(n, INVALID);
    442         _status->set(n, UNMATCHED);
     441        (*_matching)[n] = INVALID;
     442        (*_status)[n] = UNMATCHED;
    443443      }
    444444      for (NodeIt n(_graph); n != INVALID; ++n) {
     
    447447            Node v = _graph.target(a);
    448448            if ((*_matching)[v] == INVALID && v != n) {
    449               _matching->set(n, a);
    450               _status->set(n, MATCHED);
    451               _matching->set(v, _graph.oppositeArc(a));
    452               _status->set(v, MATCHED);
     449              (*_matching)[n] = a;
     450              (*_status)[n] = MATCHED;
     451              (*_matching)[v] = _graph.oppositeArc(a);
     452              (*_status)[v] = MATCHED;
    453453              break;
    454454            }
     
    470470
    471471      for (NodeIt n(_graph); n != INVALID; ++n) {
    472         _matching->set(n, INVALID);
    473         _status->set(n, UNMATCHED);
     472        (*_matching)[n] = INVALID;
     473        (*_status)[n] = UNMATCHED;
    474474      }
    475475      for(EdgeIt e(_graph); e!=INVALID; ++e) {
     
    478478          Node u = _graph.u(e);
    479479          if ((*_matching)[u] != INVALID) return false;
    480           _matching->set(u, _graph.direct(e, true));
    481           _status->set(u, MATCHED);
     480          (*_matching)[u] = _graph.direct(e, true);
     481          (*_status)[u] = MATCHED;
    482482
    483483          Node v = _graph.v(e);
    484484          if ((*_matching)[v] != INVALID) return false;
    485           _matching->set(v, _graph.direct(e, false));
    486           _status->set(v, MATCHED);
     485          (*_matching)[v] = _graph.direct(e, false);
     486          (*_status)[v] = MATCHED;
    487487        }
    488488      }
     
    498498          (*_blossom_rep)[_blossom_set->insert(n)] = n;
    499499          _tree_set->insert(n);
    500           _status->set(n, EVEN);
     500          (*_status)[n] = EVEN;
    501501          processSparse(n);
    502502        }
     
    513513          (*_blossom_rep)[_blossom_set->insert(n)] = n;
    514514          _tree_set->insert(n);
    515           _status->set(n, EVEN);
     515          (*_status)[n] = EVEN;
    516516          processDense(n);
    517517        }
     
    15491549        Value pot = (*_node_data)[bi].pot;
    15501550
    1551         _matching->set(base, matching);
     1551        (*_matching)[base] = matching;
    15521552        _blossom_node_list.push_back(base);
    1553         _node_potential->set(base, pot);
     1553        (*_node_potential)[base] = pot;
    15541554      } else {
    15551555
     
    16451645
    16461646      for (ArcIt e(_graph); e != INVALID; ++e) {
    1647         _node_heap_index->set(e, BinHeap<Value, IntArcMap>::PRE_HEAP);
     1647        (*_node_heap_index)[e] = BinHeap<Value, IntArcMap>::PRE_HEAP;
    16481648      }
    16491649      for (NodeIt n(_graph); n != INVALID; ++n) {
    1650         _delta1_index->set(n, _delta1->PRE_HEAP);
     1650        (*_delta1_index)[n] = _delta1->PRE_HEAP;
    16511651      }
    16521652      for (EdgeIt e(_graph); e != INVALID; ++e) {
    1653         _delta3_index->set(e, _delta3->PRE_HEAP);
     1653        (*_delta3_index)[e] = _delta3->PRE_HEAP;
    16541654      }
    16551655      for (int i = 0; i < _blossom_num; ++i) {
    1656         _delta2_index->set(i, _delta2->PRE_HEAP);
    1657         _delta4_index->set(i, _delta4->PRE_HEAP);
     1656        (*_delta2_index)[i] = _delta2->PRE_HEAP;
     1657        (*_delta4_index)[i] = _delta4->PRE_HEAP;
    16581658      }
    16591659
     
    16671667          }
    16681668        }
    1669         _node_index->set(n, index);
     1669        (*_node_index)[n] = index;
    16701670        (*_node_data)[index].pot = max;
    16711671        _delta1->push(n, max);
     
    27422742        Value pot = (*_node_data)[bi].pot;
    27432743
    2744         _matching->set(base, matching);
     2744        (*_matching)[base] = matching;
    27452745        _blossom_node_list.push_back(base);
    2746         _node_potential->set(base, pot);
     2746        (*_node_potential)[base] = pot;
    27472747      } else {
    27482748
     
    28322832
    28332833      for (ArcIt e(_graph); e != INVALID; ++e) {
    2834         _node_heap_index->set(e, BinHeap<Value, IntArcMap>::PRE_HEAP);
     2834        (*_node_heap_index)[e] = BinHeap<Value, IntArcMap>::PRE_HEAP;
    28352835      }
    28362836      for (EdgeIt e(_graph); e != INVALID; ++e) {
    2837         _delta3_index->set(e, _delta3->PRE_HEAP);
     2837        (*_delta3_index)[e] = _delta3->PRE_HEAP;
    28382838      }
    28392839      for (int i = 0; i < _blossom_num; ++i) {
    2840         _delta2_index->set(i, _delta2->PRE_HEAP);
    2841         _delta4_index->set(i, _delta4->PRE_HEAP);
     2840        (*_delta2_index)[i] = _delta2->PRE_HEAP;
     2841        (*_delta4_index)[i] = _delta4->PRE_HEAP;
    28422842      }
    28432843
     
    28512851          }
    28522852        }
    2853         _node_index->set(n, index);
     2853        (*_node_index)[n] = index;
    28542854        (*_node_data)[index].pot = max;
    28552855        int blossom =
  • lemon/min_cost_arborescence.h

    r559 r581  
    294294        }
    295295      }
    296       _arc_order->set(minimum.arc, _dual_variables.size());
     296      (*_arc_order)[minimum.arc] = _dual_variables.size();
    297297      DualVariable var(_dual_node_list.size() - 1,
    298298                       _dual_node_list.size(), minimum.value);
     
    336336        }
    337337      }
    338       _arc_order->set(minimum.arc, _dual_variables.size());
     338      (*_arc_order)[minimum.arc] = _dual_variables.size();
    339339      DualVariable var(node_bottom, _dual_node_list.size(), minimum.value);
    340340      _dual_variables.push_back(var);
     
    365365        Node source = _heap->top();
    366366        _heap->pop();
    367         _node_order->set(source, -1);
     367        (*_node_order)[source] = -1;
    368368        for (OutArcIt it(*_digraph, source); it != INVALID; ++it) {
    369369          if ((*_arc_order)[it] < 0) continue;
     
    651651      for (NodeIt it(*_digraph); it != INVALID; ++it) {
    652652        (*_cost_arcs)[it].arc = INVALID;
    653         _node_order->set(it, -3);
    654         _heap_cross_ref->set(it, Heap::PRE_HEAP);
     653        (*_node_order)[it] = -3;
     654        (*_heap_cross_ref)[it] = Heap::PRE_HEAP;
    655655        _pred->set(it, INVALID);
    656656      }
    657657      for (ArcIt it(*_digraph); it != INVALID; ++it) {
    658658        _arborescence->set(it, false);
    659         _arc_order->set(it, -1);
     659        (*_arc_order)[it] = -1;
    660660      }
    661661      _dual_node_list.clear();
  • lemon/preflow.h

    r559 r581  
    405405      _phase = true;
    406406      for (NodeIt n(_graph); n != INVALID; ++n) {
    407         _excess->set(n, 0);
     407        (*_excess)[n] = 0;
    408408      }
    409409
     
    418418
    419419      std::vector<Node> queue;
    420       reached.set(_source, true);
     420      reached[_source] = true;
    421421
    422422      queue.push_back(_target);
    423       reached.set(_target, true);
     423      reached[_target] = true;
    424424      while (!queue.empty()) {
    425425        _level->initNewLevel();
     
    430430            Node u = _graph.source(e);
    431431            if (!reached[u] && _tolerance.positive((*_capacity)[e])) {
    432               reached.set(u, true);
     432              reached[u] = true;
    433433              _level->initAddItem(u);
    434434              nqueue.push_back(u);
     
    445445          if ((*_level)[u] == _level->maxLevel()) continue;
    446446          _flow->set(e, (*_capacity)[e]);
    447           _excess->set(u, (*_excess)[u] + (*_capacity)[e]);
     447          (*_excess)[u] += (*_capacity)[e];
    448448          if (u != _target && !_level->active(u)) {
    449449            _level->activate(u);
     
    479479        }
    480480        if (excess < 0 && n != _source) return false;
    481         _excess->set(n, excess);
     481        (*_excess)[n] = excess;
    482482      }
    483483
     
    488488
    489489      std::vector<Node> queue;
    490       reached.set(_source, true);
     490      reached[_source] = true;
    491491
    492492      queue.push_back(_target);
    493       reached.set(_target, true);
     493      reached[_target] = true;
    494494      while (!queue.empty()) {
    495495        _level->initNewLevel();
     
    501501            if (!reached[u] &&
    502502                _tolerance.positive((*_capacity)[e] - (*_flow)[e])) {
    503               reached.set(u, true);
     503              reached[u] = true;
    504504              _level->initAddItem(u);
    505505              nqueue.push_back(u);
     
    509509            Node v = _graph.target(e);
    510510            if (!reached[v] && _tolerance.positive((*_flow)[e])) {
    511               reached.set(v, true);
     511              reached[v] = true;
    512512              _level->initAddItem(v);
    513513              nqueue.push_back(v);
     
    525525          if ((*_level)[u] == _level->maxLevel()) continue;
    526526          _flow->set(e, (*_capacity)[e]);
    527           _excess->set(u, (*_excess)[u] + rem);
     527          (*_excess)[u] += rem;
    528528          if (u != _target && !_level->active(u)) {
    529529            _level->activate(u);
     
    537537          if ((*_level)[v] == _level->maxLevel()) continue;
    538538          _flow->set(e, 0);
    539           _excess->set(v, (*_excess)[v] + rem);
     539          (*_excess)[v] += rem;
    540540          if (v != _target && !_level->active(v)) {
    541541            _level->activate(v);
     
    578578              if (!_tolerance.less(rem, excess)) {
    579579                _flow->set(e, (*_flow)[e] + excess);
    580                 _excess->set(v, (*_excess)[v] + excess);
     580                (*_excess)[v] += excess;
    581581                excess = 0;
    582582                goto no_more_push_1;
    583583              } else {
    584584                excess -= rem;
    585                 _excess->set(v, (*_excess)[v] + rem);
     585                (*_excess)[v] += rem;
    586586                _flow->set(e, (*_capacity)[e]);
    587587              }
     
    601601              if (!_tolerance.less(rem, excess)) {
    602602                _flow->set(e, (*_flow)[e] - excess);
    603                 _excess->set(v, (*_excess)[v] + excess);
     603                (*_excess)[v] += excess;
    604604                excess = 0;
    605605                goto no_more_push_1;
    606606              } else {
    607607                excess -= rem;
    608                 _excess->set(v, (*_excess)[v] + rem);
     608                (*_excess)[v] += rem;
    609609                _flow->set(e, 0);
    610610              }
     
    616616        no_more_push_1:
    617617
    618           _excess->set(n, excess);
     618          (*_excess)[n] = excess;
    619619
    620620          if (excess != 0) {
     
    651651              if (!_tolerance.less(rem, excess)) {
    652652                _flow->set(e, (*_flow)[e] + excess);
    653                 _excess->set(v, (*_excess)[v] + excess);
     653                (*_excess)[v] += excess;
    654654                excess = 0;
    655655                goto no_more_push_2;
    656656              } else {
    657657                excess -= rem;
    658                 _excess->set(v, (*_excess)[v] + rem);
     658                (*_excess)[v] += rem;
    659659                _flow->set(e, (*_capacity)[e]);
    660660              }
     
    674674              if (!_tolerance.less(rem, excess)) {
    675675                _flow->set(e, (*_flow)[e] - excess);
    676                 _excess->set(v, (*_excess)[v] + excess);
     676                (*_excess)[v] += excess;
    677677                excess = 0;
    678678                goto no_more_push_2;
    679679              } else {
    680680                excess -= rem;
    681                 _excess->set(v, (*_excess)[v] + rem);
     681                (*_excess)[v] += rem;
    682682                _flow->set(e, 0);
    683683              }
     
    689689        no_more_push_2:
    690690
    691           _excess->set(n, excess);
     691          (*_excess)[n] = excess;
    692692
    693693          if (excess != 0) {
     
    732732      typename Digraph::template NodeMap<bool> reached(_graph);
    733733      for (NodeIt n(_graph); n != INVALID; ++n) {
    734         reached.set(n, (*_level)[n] < _level->maxLevel());
     734        reached[n] = (*_level)[n] < _level->maxLevel();
    735735      }
    736736
     
    740740      std::vector<Node> queue;
    741741      queue.push_back(_source);
    742       reached.set(_source, true);
     742      reached[_source] = true;
    743743
    744744      while (!queue.empty()) {
     
    750750            Node v = _graph.target(e);
    751751            if (!reached[v] && _tolerance.positive((*_flow)[e])) {
    752               reached.set(v, true);
     752              reached[v] = true;
    753753              _level->initAddItem(v);
    754754              nqueue.push_back(v);
     
    759759            if (!reached[u] &&
    760760                _tolerance.positive((*_capacity)[e] - (*_flow)[e])) {
    761               reached.set(u, true);
     761              reached[u] = true;
    762762              _level->initAddItem(u);
    763763              nqueue.push_back(u);
     
    793793            if (!_tolerance.less(rem, excess)) {
    794794              _flow->set(e, (*_flow)[e] + excess);
    795               _excess->set(v, (*_excess)[v] + excess);
     795              (*_excess)[v] += excess;
    796796              excess = 0;
    797797              goto no_more_push;
    798798            } else {
    799799              excess -= rem;
    800               _excess->set(v, (*_excess)[v] + rem);
     800              (*_excess)[v] += rem;
    801801              _flow->set(e, (*_capacity)[e]);
    802802            }
     
    816816            if (!_tolerance.less(rem, excess)) {
    817817              _flow->set(e, (*_flow)[e] - excess);
    818               _excess->set(v, (*_excess)[v] + excess);
     818              (*_excess)[v] += excess;
    819819              excess = 0;
    820820              goto no_more_push;
    821821            } else {
    822822              excess -= rem;
    823               _excess->set(v, (*_excess)[v] + rem);
     823              (*_excess)[v] += rem;
    824824              _flow->set(e, 0);
    825825            }
     
    831831      no_more_push:
    832832
    833         _excess->set(n, excess);
     833        (*_excess)[n] = excess;
    834834
    835835        if (excess != 0) {
  • test/kruskal_test.cc

    r440 r581  
    100100        "Total cost should be 10");
    101101
    102   edge_cost_map.set(e1, -10);
    103   edge_cost_map.set(e2, -9);
    104   edge_cost_map.set(e3, -8);
    105   edge_cost_map.set(e4, -7);
    106   edge_cost_map.set(e5, -6);
    107   edge_cost_map.set(e6, -5);
    108   edge_cost_map.set(e7, -4);
    109   edge_cost_map.set(e8, -3);
    110   edge_cost_map.set(e9, -2);
    111   edge_cost_map.set(e10, -1);
     102  edge_cost_map[e1] = -10;
     103  edge_cost_map[e2] = -9;
     104  edge_cost_map[e3] = -8;
     105  edge_cost_map[e4] = -7;
     106  edge_cost_map[e5] = -6;
     107  edge_cost_map[e6] = -5;
     108  edge_cost_map[e7] = -4;
     109  edge_cost_map[e8] = -3;
     110  edge_cost_map[e9] = -2;
     111  edge_cost_map[e10] = -1;
    112112
    113113  vector<Edge> tree_edge_vec(5);
Note: See TracChangeset for help on using the changeset viewer.