Correct option is (a) Rd || {Rs*(1+gm * ro) + ro}
Easy explanation: The impedance looking into the drain of M1 would be similar to that of a CS stage with source degeneration. That implies the impedance is Rs*(1+gm * ro) + ro. But this impedance is parallel to Rd. Hence, the overall output impedance becomes Rd || {Rs*(1+gm* ro) + ro}.