方程组x+y=k+3,x-y=3k+1的解是一对正数,求k的取值范围