We have studied thermal expansion and contraction of polystyrene thin film. supported on silicon substrate using Xray reflectivity and inelastic neutron scattering techniques. In annealing experiments, we found fast and very slow contraction processes in the thin films above the glass transition temperature. The fast one was attributed to normal relaxation (annealing) process observed in bulk and responsible partly for the reported apparent negative expansivity. In addition, we found extremely slow re-expansion process in the glassy state. The ultra-slow contraction and re-expansion have been assigned to the lateral expansion and contraction of the film, respectively. Using well-annealed films, we have evaluated thermal expansivity as a function of film thickness, and found that it decreases with film thickness in a thickness range below twice of radius of gyration of a polymer chain 2R(g) in the glassy state. The decrease has been attributed to hardening of harmonic force constant due to confinement in a thin film. In order to confirm this idea we have performed inelastic neutron scattering on the thin films in meV region and found the harmonic force constant decreased with the film thickness. We also evaluated the mean square displacement with an energy resolution of 25 mu eV to find that the glass transition temperature T-g increases with decreasing the film thickness, which contradicts the results of the XR. This has been tentatively assigned to a wide distribution of relaxation time in the thin films. In addition, we have investigated the dynamic anisotropy as well as the dynamic heterogeneity in the glassy state.