def elementary(n, k):
R = PolynomialRing(ZZ, ['x%d' % (i+1) for i in range(n)])
f = SymmetricFunctions(QQ).e()[k]
return R(f.expand(n, alphabet=[str(g) for g in R.gens()]))
elementary(7, 3) # more variables than this table holds