The bigravity models coupled with two scalar fields are constructed. We show that a wide class of the expansion history of the universe, especially corresponding to dark energy and/or inflation, can be described by a solution of the bigravity model. We discuss the stability of the solution and give the conditions for the stability. We also explicitly construct a model which gives a stable solution. By using the stable model, for an arbitrary evolution of the universe expansion, we construct the Brans-Dicke like model which reproduces the evolution.

We study the accelerating cosmology in massive $F(R)$ bigravity via the reconstruction scheme. The consistent solution of the FRW equations is presented: it includes Big and Little Rip, quintessence, de Sitter and decelerating universes described by the physical $g$ metric while the corresponding solution of the universe described by the reference $f$ metric is also found. It is demonstrated that in general the cosmological singularities of $g$ metric are not always manifested as cosmological singularities of the reference $f$ metric. We give two consistent ways to describe the Big and Little Rip, quintessence, de Sitter and decelerating universes. In one of the consistent solutions, the two metrics $g$ and $f$ coincide with each other, which may indicate the connection with the convenient single metric background formulation. For another solution, where two metrics $g$ and $f$ do not coincide with each other, here always appears super-luminal mode.

Explicit formulae of the equations in the generalized Galileon models are given. We also develop the formulation of the reconstruction. By using the formulation, we can explicitly construct an action which reproduces an arbitrary development of the expansion of the universe. The conditions how the reconstructed solution becomes stable and therefore it becomes an attractor solution are also given. Working in the static and spherically symmetric space-time, we investigate how the Vainshtein mechanism works in the generalized Galileon model and the correction to the Newton law becomes small. It is also shown that any spherically symmetric and static geometry can be realized by properly choosing the form of the action, which may tell that the solution could have fourth hair corresponding to the scalar field.