Как группировать термины по экспоненте в выражении Mathematica

мне нужно найти метод преобразования выражения типа

a^(1+m+n) b^(2+2m - 2n)

на

(a b^2)^m (a/b^2)^n (a b^2),

то есть группировать члены с одинаковым показателем. Я попытался использовать Collect [] и т. д., Но ничего не могу заставить работать.

какие предложения?

спасибо, Тома

7
wolfram-mathematica

2 ответов

используя Log в сочетании с CoefficientRules:

exp = a^(1 + m + n) b^(2 + 2 m - 2 n);

Times @@ (Exp[#[[2]]]^(Times @@ ({n, m}^#[[1]])) & /@ 
   CoefficientRules[PowerExpand[Log[exp]], {n, m}])

выход:

a (a/b^2)^n b^2 (a b^2)^m
10

Вы можете сделать это, например:

log[x_*y_] := log[x] + log[y];
log[x_^y_] := y*log[x];
log1 /: a_*log1[b_] := log1[b^a];
log1 /: Plus[x__log1] := log1[Times @@ Map[First, {x}]];
exp[HoldPattern[Plus[x__]]] := Times @@ Map[exp, {x}];
exp[log1[x_]] := x

и затем:

In[58]:= exp[Collect[Expand[log[a^(1+m+n) b^(2+2m-2n)]],{m,n}]]/.log->log1

Out[58]= a (a/b^2)^n b^2 (a b^2)^m
7