# polynomial ring variable limit workaround

On my machine, the maximum number of variables allowed in a polynomial ring is 32767. For instance,

S = PolynomialRing(QQ, 32767,'s')

works, but

S = PolynomialRing(QQ, 32768,'s')

does not.

Is there a workaround with enough functionality to allow me to do something like compute a primary decomposition of an ideal?