def R_mon(*args):
R"""Return the monomial $\prod R_{ij}$.
args is i1, j1, i2, j2, ..."""
result = DualSteenrodDense.unit()
for i in range(len(args) // 2):
s, t = args[2 * i], args[2 * i + 1]
result *= DualSteenrodDense.gen(t - s, 2**s)
for m in result.data:
return m
def R(i, j):
"""Return $R_{ij}$"""
return DualSteenrodDense.gen(j - i, 2**i)