A Multistage scenario Cluster Lagrangian Decomposition (MCLD) approach for obtaining strong lower bounds on the solution value of large sized instances of the multistage stochastic mixed 0–1 problem is presented.

MCLD frequently gives the solution value of the original stochastic mixed 0–1 problem by using merged or pure scenario cluster submodels either without Lagrange multipliers usage or considering any of the four Lagrange multipliers schemes that we have selected to work with.

It is easy to prove that the performance of the MCLD procedure outperforms the traditional Lagrangian Decomposition scheme based on single scenarios in both the bound׳s quality and elapsed time.

The efficiency of the four updating schemes, as contrasted in the testbed we have experimented with, is very similar in quality of the MCLD bound. However, a careful implementation of the Subgradient Method makes it more robust than the other three updating schemes. Notice that the SM based MCLD approach obtains the solution value in the original problem in thirteen out of the fourteen instances included in the testbed.

The proposed approach requires an elapsed time that is one order of magnitude smaller than the time required by CPLEX, in the worst case.