It is well known from the seminal paper by Fujita  for 1 < p < p0, and Hayakawa  for the critical case p = p0, that all the solutions u ≥ 0 of the semilinear heat equation ut =δu + |u|p-1 u in ℝN × ℝ+, in the range 1<p≤p0 = 1 + N 2, (0.1) with arbitrary initial data u0(x) ≥ 0, ≢ 0, blow-up in finite time, while for p > p0 there exists a class of sufficiently "small" global in time solutions. This fundamental result from the 1960-70s (see also  for related contributions), was a cornerstone of further active blow-up research. Nowadays, similar Fujita-type critical exponents p0, as important characteristics of stability, unstability, and blow-up of solutions, have been calculated for various nonlinear PDEs. The above blow-up conclusion does not include solutions of changing sign, so some of them may remain global even for p ≤ p0. Our goal is a thorough description of blow-up and global in time oscillatory solutions in the subcritical range in (0.1) on the basis of various analytic methods including nonlinear capacity, variational, category, fibering, and invariant manifold techniques. Two countable sets of global solutions of changing sign are shown to exist. Most of them are not radially symmetric in any dimension N ≥ 2 (previously, only radial such solutions in ℝN or in the unit ball B1 ⊆ ℝN were mostly studied). A countable sequence of critical exponents, at which the whole set of global solutions changes its structure, is detected: pl = 1 + N 2 +l, 1 = 0, 1, 2, See [47, 48] for earlier interesting contributions on sign changing solutions.