Pairing operators of electrons with an arbitrary distance on a square lattice are defined in each possible wave channel, {ital s}, {ital p}, {ital d}{sub {ital x}}{sup 2}{minus}{ital y}{sup 2}, {ital d}{sub {ital x}{ital y}}, and {ital g}{sub {ital x}}{sup 3}{ital y}{minus}{ital xy}{sup 3}. Exact identities among pairing amplitudes in the Hubbard model are presented, giving powerful constraints on the nature of pairings. At half filling, some of the pairing amplitudes, including those of the second-shortest pairing distance, vanish identically in each wave channel except the {ital p} wave. There are correlations among pairing amplitudes away from half filling, as well as nonidentically vanishing ones at half filling. Even if the pairing amplitudes vanish up to a certain distance, those of larger distances can still possibly be nonzero. The order of magnitude is estimated for the pairine amplitudes, indicating that those of large distances can by no means be neglected. The scheme is also used in the extended Hubbard model including nearest-neighbor density-density interaction and that including next-nearest-neighbor hopping. In that including nearest-neighbor interaction, the features of pairing amplitudes are just the same as in the original Hubbard model at half filling, and are also similar away from half filling,more » but with the difference that the amplitude is now dependent on the nearest-neighbor interaction. In the model including next-nearest-neighbor hopping, the features are quite different, the constraint identities are more underconstrained, and the amplitude has no monotonic dependence upon the pairing distance if the next-nearest-neighbor hopping has an opposite sign to the nearest-neighbor hopping as required physically. The same scheme is also used to discuss the possible pairings in the coupled two-plane Hubbard model and in the infinite-layer Hubbard model.« less