You'd need to be a little more careful than that. The coefficients do matter:
a*b+c*d
is invariant under more permutations than
a*b-c*d
is.
A brain-dead way is to just iterate through all permutations, which is entirely reasonable for 4 variables and infeasible by the time you hit 14 variables or so:
sage: P.<a,b,c,d>=QQ[]
sage: G=SymmetricGroup([a,b,c,d])
sage: f=a*b+c*d
sage: [g for g in G if f(g(a),g(b),g(c),g(d)) == f]
[(), (c,d), (a,b), (a,b)(c,d), (a,c)(b,d), (a,c,b,d), (a,d,b,c), (a,d)(b,c)]