The correct answer is 4x the power. Let me give an example to answer your comment (on Tesla23's great answer) on whether an antenna pattern changes with the number of antennas.
Assume the antennas in your image are separated by distance d and that both antennas are isotropic radiators so that radiated field for each INDIVIDUAL antenna is \$ E_{iso} = \dfrac{e^{-jkr}}{4\pi |r|} \$ where \$r \$ is a vector pointing from the antenna to any location we wish to know the fields at. In your case, the fields at the receiving antenna are of primary interest.
The image above is the scenario you described, but now I have added two vectors, one showing the path of fields from TX1 to the RX \$ r_1 \$ and one showing the path from TX2 to the RX \$ r_2 \$. Furthermore, assume that each antenna is excited with power \$ P_{TX} \$. The field associated with this power will be \$ E_{TX} \propto \sqrt{P_{TX}} \$. Now then, the fields at the receiver due to both antennas is \$ E_{RX} = E_{TX1}+E_{TX2} = E_{TX}\dfrac{e^{-jkr_1}}{4\pi |r_1|}+E_{TX}\dfrac{e^{-jkr_2}}{4\pi |r_2|} \$
Since in the given situation \$r_1 = r_2 \$ then the equation is simplified as
\$ E_{RX} = 2 E_{TX}\dfrac{e^{-jkr_1}}{4\pi |r_1|} \$
This is to say that the received fields at the receiver are double what they would have been for one transmitter. Since \$ E_{TX} \propto \sqrt{P_{TX}} \$, the power received is quadrupled due to the doubling of the field intensity. You are right to wonder where this power comes from, and that is what I will answer next.
An ideal isotropic antenna radiates fields in all directions equally. So no matter where I place my receive antenna, the expected gain would be constant for a SINGLE isotropic element. But now consider the case where there are two. Assume now for simplicity that the separation between antennas is \$ d = \lambda/2 \$. Now also assume that I placed an antenna to the left of TX1, as shown below.
Now the fields at the RX are again given by
\$ E_{RX} = E_{TX1}+E_{TX2} = E_{TX}\dfrac{e^{-jkr_1}}{4\pi |r_1|}+E_{TX}\dfrac{e^{-jkr_2}}{4\pi |r_2|} \$
In the case above \$ r_2 = r_1 + d \$. if I assume that \$ |r_1| \approx |r_2| \$ which is equivalent to \$ |r_1| >> d \$ then the total fields at the receiver are given by
\$ E_{RX} = E_{TX1}+E_{TX2} = E_{TX}\dfrac{e^{-jkr_1}}{4\pi |r_1|}+E_{TX}\dfrac{e^{-jk(r_1+d)}}{4\pi |r_1|} \$
Which is equivalent to
\$ E_{RX} = E_{TX}\dfrac{e^{-jk(r_1)}}{4\pi |r_1|}(1+e^{-jk(d)}) \$
Using the definition of k
\$ E_{RX} = E_{TX}\dfrac{e^{-jk(r_1)}}{4\pi |r_1|}(1-1) = 0 \$
So the fields for this RX location would be 0. This is the essence of where the power is coming from. The fields add up constructively in some locations (RX location 1) and destructively in some locations (RX location 2). So if I integrated along the entire polar sphere, the total power would still be the same, so power is still conserved, but because of these field interactions at a specific location, the power can change drastically (as demonstrated).
Antenna Theory by Constantine Balanis has a whole chapter dedicated to this subject. Phased arrays are a very interesting field I would encourage you to look into. Finally, one disclaimer I should make about this simple analysis is that I assumed the power radiated by each antenna did not change by their proximity to one another. Generally, this is not true, and antennas will "couple." This coupling changes how each antenna radiates its power and adds another complexity to antenna analysis.