Find all polynomials $P$ with real coefficients such that for all $x, y ,z \in R$,
$$P(x)+P(y)+P(z)+P(x+y+z)=P(x+y)+P(y+z)+P(z+x)$$
We are given the functional equation for a polynomial \( P \) with real coefficients:
\[ P(x) + P(y) + P(z) + P(x+y+z) = P(x+y) + P(y+z) + P(z+x) \]
for all \( x, y, z \in \mathbb{R} \).

1. **Setting \( x = y = z \):**
   \[
   P(x) + P(x) + P(x) + P(3x) = P(2x) + P(2x) + P(2x)
   \]
   Simplifying, we get:
   \[
   3P(x) + P(3x) = 3P(2x)
   \]

2. **Comparing leading coefficients:**
   Let \( P(x) = a_d x^d + \cdots \) where \( a_d \) is the leading coefficient and \( d \) is the degree of the polynomial. The leading term of \( 3P(x) \) is \( 3a_d x^d \), and the leading term of \( P(3x) \) is \( a_d (3x)^d = a_d \cdot 3^d x^d \). Similarly, the leading term of \( 3P(2x) \) is \( 3a_d \cdot 2^d x^d \).

   Equating the leading terms, we get:
   \[
   3a_d + a_d \cdot 3^d = 3a_d \cdot 2^d
   \]
   Simplifying, we have:
   \[
   3 + 3^d = 3 \cdot 2^d
   \]
   Dividing both sides by 3:
   \[
   1 + 3^{d-1} = 2^d
   \]

3. **Solving for \( d \):**
   We need to find \( d \) such that \( 1 + 3^{d-1} = 2^d \). Testing small values of \( d \):
   - For \( d = 1 \):
     \[
     1 + 3^0 = 2^1 \implies 1 + 1 = 2 \quad \text{(True)}
     \]
   - For \( d = 2 \):
     \[
     1 + 3^1 = 2^2 \implies 1 + 3 = 4 \quad \text{(True)}
     \]
   - For \( d \geq 3 \):
     \[
     3^{d-1} > 2^d - 1 \quad \text{(False for all \( d \geq 3 \))}
     \]

   Therefore, \( d \) can only be 1 or 2.

4. **Case \( d = 1 \):**
   Let \( P(x) = ax + b \). Substituting into the original equation:
   \[
   a(x) + b + a(y) + b + a(z) + b + a(x+y+z) + b = a(x+y) + b + a(y+z) + b + a(z+x) + b
   \]
   Simplifying, we get:
   \[
   ax + ay + az + a(x+y+z) + 4b = ax + ay + az + a(x+y) + a(y+z) + a(z+x) + 3b
   \]
   This simplifies to:
   \[
   a(x+y+z) + b = a(x+y+z) + b
   \]
   This is always true, so \( P(x) = ax + b \) is a solution.

5. **Case \( d = 2 \):**
   Let \( P(x) = ax^2 + bx + c \). Substituting into the original equation:
   \[
   a(x^2) + bx + c + a(y^2) + by + c + a(z^2) + bz + c + a((x+y+z)^2) + b(x+y+z) + c = a((x+y)^2) + b(x+y) + c + a((y+z)^2) + b(y+z) + c + a((z+x)^2) + b(z+x) + c
   \]
   Simplifying, we get:
   \[
   ax^2 + ay^2 + az^2 + a(x^2 + y^2 + z^2 + 2xy + 2yz + 2zx) + b(x+y+z) + 4c = a(x^2 + 2xy + y^2) + a(y^2 + 2yz + z^2) + a(z^2 + 2zx + x^2) + b(x+y) + b(y+z) + b(z+x) + 3c
   \]
   Simplifying further, we get:
   \[
   2ax^2 + 2ay^2 + 2az^2 + 2a(xy + yz + zx) + b(x+y+z) + 4c = 2ax^2 + 2ay^2 + 2az^2 + 2a(xy + yz + zx) + b(x+y+z) + 3c
   \]
   This simplifies to:
   \[
   4c = 3c
   \]
   Therefore, \( c = 0 \). So, \( P(x) = ax^2 + bx \) is a solution.

Conclusion:
The polynomials that satisfy the given functional equation are of the form \( P(x) = ax + b \) or \( P(x) = ax^2 + bx \).

The final answer is \( \boxed{ P(x) = ax + b } \) or \( P(x) = ax^2 + bx \).