ns = seq(1, 25)
ms = seq(1, 21)

four_plus_one = 4*ns+1
five_minus_one = 5*ms-1

print(intersect(four_plus_one, five_minus_one))