The four-pomeron vertex is studied in the perturbative QCD. Its dominating terms of the leading (zeroth and first) orders in the coupling constant and subdominant in the number of colors are constructed. The vertex consists of two terms, one with a derivative in rapidity ∂y and the other with the BFKL interaction between pomerons. The corresponding part of the action and equations of motion are found. The iterative solution of the latter is possible only for rapidities smaller than 2 and quite large coupling constant αs, of the order or greater than unity, when the quadruple pomeron interaction is relatively small. Also iteration of the part with ∂y is unstable in the infrared region and compels to introduce an infrared cut. The variational approach with simple trying functions allows to find the minimum of the action at αs of the order 0.2 and rapidities up to 25. Numerical estimates for O–O collisions show that actually the influence of the quadruple pomeron interaction turns out to be rather small.