设 $y^2 \mathrm{~d} x+(2 x y+1) \mathrm{d} y$ 是函数 $f(x, y)$ 的全微分, 其中 $f(0,0)=0$, 求 $f(x, y)$, 并计算曲面积分 $I=\iint_{\Sigma} z f(x, y) \mathrm{d} S$, 其中 $\Sigma$ 是椎面 $z=\sqrt{x^2+y^2}$ 被柱面 $x^2+(y-1)^2=1$ 所截下的有限部分.
$\text{A.}$
$\text{B.}$
$\text{C.}$
$\text{D.}$