High order time discretization methods with the strong stability property