Clear["Global`*"];Subscript[q, r] =
Q - (Subscript[p, r] - Subscript[p, e])/(1 - \[Theta]);Subscript[q,
e] = (\[Theta] Subscript[p, r] - Subscript[p,
e])/(\[Theta] (1 - \[Theta]));Subscript[\[Pi], m] =
w Subscript[q, r] +
Subscript[p, e] Subscript[q,
e] - (c - Subscript[o, 2] l) (Subscript[q, r] + Subscript[q, e]);
Subscript[\[Pi],
r] = (Subscript[p, r] - w) Subscript[q, r] +
Subscript[o, 1] l (Subscript[q, r] + Subscript[q, e]) - kl^2;
FullSimplify[
Solve[D[Subscript[\[Pi], r], Subscript[p, r]] == 0, Subscript[p, r]]]

好的谢谢～可以出结果了
但是下一步结果又出不来了，可以再请您看一下吗？
Clear["Global`*"];
Subscript[q, r] =
Q - (Subscript[p, r] - Subscript[p, e])/(1 - \[Theta]);
Subscript[q,
e] = (\[Theta] Subscript[p, r] - Subscript[p,
e])/(\[Theta] (1 - \[Theta]));
Subscript[\[Pi], m] =
w Subscript[q, r] +
Subscript[p, e] Subscript[q,
e] - (c - Subscript[o, 2] l) (Subscript[q, r] + Subscript[q, e]);
Subscript[\[Pi],
r] = (Subscript[p, r] - w) Subscript[q, r] +
Subscript[o, 1] l (Subscript[q, r] + Subscript[q, e]) - kl^2;
FullSimplify[
Solve[D[Subscript[\[Pi], r], Subscript[p, r]] == 0, Subscript[p, r]]]
FullSimplify[Solve[D[Subscript[\[Pi], r], l] == 0, l]]