We analyze the ␣ decay between ground states along N − Z chains in deformed heavy and superheavy nuclei, by using the pairing approach. We show that the derivative of the preformation amplitude is practically a constant along any ␣ chain, while that of the outgoing wave function changes exponentially upon the Coulomb parameter. This leads to the breakdown of the continuity equation and therefore to wrong decay widths. The behavior cannot be explained within the standard shell model. We significantly correct this deficiency by considering an ␣-cluster factor in the preformation amplitude, depending exponentially upon the Coulomb parameter. Thus, four-body correlations, connected with the radial shape of the preformation factor, are directly evidenced by the ␣-decay systematics. Moreover, this procedure, in principle, fully determines the Q value and is an important development in the ␣-decay theory. It also allows us to analyze the relative ␣-clustering structure of the emitter. It turns out that the isotopes close to the region N Ͼ 126 and superheavy nuclei have a stronger clustering behavior. For superheavy region an additional dependence upon the number of interacting ␣ particles is necessary.