Now asked on MO here
Let $f_n $ be the $ n $-th Faulhaber polynomial. In this question, I observed the following divisibility properties:
- If $ n $ is even, then $ f_n $ is divisible by $ f_2 $.
- If $ n > 1 $ is odd, then $ f_n $ is divisible by $ f_3 $.
However, it seems that Faulhaber polynomials do not share any additional common factors (ignoring the denominator ). More precisely, for $ n \neq m $ and $n,m>1$ , we have
$$
\gcd(f_n, f_m) =
\begin{cases}
f_3, & \text{if both } n \text{ and } m \text{ are odd,} \\
f_2, & \text{if both } n \text{ and } m \text{ are even,} \\
f_1, & \text{otherwise.}
\end{cases}
$$
Beyond these cases, no additional common factors seem to exist.
Added: As @sil mentioned in the comments $\gcd(f_9,f_{10})= n (n+1) \left(n^2+n-1\right)$
However for all $n,m\le100$ that is the only counter example that I found.
I used this ugly Mathematica code to generate the gcds
For[m = 1, m < 101, m++,
For[j = 1, j < m + 1, j++,
If[ m != j,
Print[m, j,
Factor[PolynomialGCD[Sum[k^m, {k, 1, n}],
Sum[k^j, {k, 1, n}] ]]] , ]]]
Added: Using the code in this answer I have checked that for $n,m<1000$ there are no other counter examples.
Here are the first 14 Faulhaber's Polynomials: $$f_1=\frac{n(n+1)}{2}$$ $$ f_{2}=\frac16 n (1+n) (1+2 n)$$ $$f_{3}=\frac14 n^2 (1+n)^2$$ $$f_{4}=\frac1{30} n (1+n) (1+2 n) (-1+3 n+3 n^2)$$ $$f_{5}=\frac1{12} n^2 (1+n)^2 (-1+2 n+2 n^2)$$ $$f_{6}=\frac1{42} n (1+n) (1+2 n) (1-3 n+6 n^3+3 n^4)$$ $$f_{7}=\frac1{24} n^2 (1+n)^2 (2-4 n-n^2+6 n^3+3 n^4)$$ $$f_{8}=\frac1{90} n (1+n) (1+2 n) (-3+9 n-n^2-15 n^3+5 n^4+15 n^5+5 n^6)$$ $$f_{9}=\frac1{20} n^2 (1+n)^2 (-1+n+n^2) (3-3 n-n^2+4 n^3+2 n^4)$$ $$f_{10}=\frac1{66} n (1+n) (1+2 n) (-1+n+n^2) (-5+10 n+3 n^2-11 n^3+2 n^4+9 n^5+3 n^6)$$ $$f_{11}=\frac1{24} n^2 (1+n)^2 (10-20 n-3 n^2+26 n^3-5 n^4-16 n^5+4 n^6+8 n^7+2 n^8)$$ $$f_{12}=\frac1{2730}n (1+n) (1+2 n) (-691+2073 n-287 n^2-3285 n^3+1420 n^4+2310 n^5-1190 n^6-1050 n^7+525 n^8+525 n^9+105 n^{10})$$ $$f_{13}=\frac1{420} n^2 (1+n)^2 (-691+1382 n+202 n^2-1786 n^3+367 n^4+1052 n^5-326 n^6-400 n^7+125 n^8+150 n^9+30 n^{10})$$ $$f_{14}=\frac1{90} n (1+n) (1+2 n) (105-315 n+44 n^2+498 n^3-217 n^4-345 n^5+182 n^6+144 n^7-81 n^8-45 n^9+24 n^{10}+18 n^{11}+3 n^{12})$$