Both the equilibrium and transport properties of the vortex matter are essentially affected by the behavior of the intervortex interaction potential. In isotropic bulk superconductors this potential is well known to be repulsive and is screened at intervortex distances R greater than the London penetration depth λ. As a result, in perfect crystals quantized Abrikosov vortices form a triangular lattice. In thin films of anisotropic superconductors this standard interaction potential behavior appears to be strongly modified because of the interplay between the long-ranged repulsion predicted in the pioneering work by J. Pearl and the attraction caused by the tilt of the vortex lines with respect to the anisotropy axes. This interplay results in a new type of vortex arrangement formed by finite-size vortex chains, i.e., vortex molecules. Tilted vortices with such unusual interaction potential form clusters with the size depending on the field tilting angle and film thickness or/and can arrange into multiquanta flux lattice. The magnetic flux through the unit cells of the corresponding flux line lattices equals to an integer number N of flux quanta. Thus, the increase in the field tilting (or varying temperature) should be accompanied by the series of the phase transitions between the vortex lattices with different N. A similar scenario should be realized in strongly anisotropic BSCCO high-T c superconductors where in tilted field a crossing lattice of Abrikosov vortices (the stacks of pancakes in this case) and Josephson vortices appears. This crossing leads to the zigzag deformation of the pancakes stacks which is responsible for the attraction interaction competing with the long-ranged Pearl’s repulsion.