The optimization problem in this paper is targeted at large-scale hydrothermal power systems. The thermal part of the system is a multi-area power pool with tie-line constraints, and the hydro part is a set of cascaded hydrostations. The objective is to minimize the operation cost of the thermal subsystem. This is an integer nonlinear optimization process with a large number of variables and constraint. In order to obtain the optimal solution in a reasonable time, we decompose the problem into thermal and hydro subproblems. The coordinator between these subproblems is the system Lagrange multiplier. For the thermal subproblem, in a multi-area power pool, it is necessary to coordinate the area generations for reducing the operation cost without violating tie limits. For the hydro subsystem, network flow concepts are adopted to coordinate water usage over the entire study time span, and the reduced gradient method is used to overcome the linear characteristic of the network flow method in order to obtain the optimal solution. In this study, load forecasting errors and forced outages of generating units are incorporated in system reliability requirements. Three case studies for the proposed method are presented.