line = y == 2 x - 1;
curve = Sin[x + 2 y] == 2 x - y;
tangentSystem = Flatten[
{#, (* equations *)
Dt@#, (* differentials *)
Dt[x] != 0}] &@ {line, curve}
Reduce[tangentSystem,
{Dt[x], Dt[y]}] (* if solution(s), then tangent;
if False, then not tangent;
if Reduce[...], then the problem has not been solved *)
tangentPoints = {x, y} /. Solve[ (* bound to find explicit solutions *)
Flatten[{tangentSystem, -10 < y < 10, -10 < x < 10}],
{Dt[x], Dt[y], x, y}]