与a(3^e)=3^e相乘,a(p^e)=((p-1)*e+p)*p^(e-1)如果p mod 3=1,a(p^e)=p^(2*地板(e/2)),如果p mod 3=2。-弗拉·约德维奇2003年9月27日
枫木
A087694号:=proc(n)选项记得;当地pf,p,f,e;如果n=1,那么1就1;如果nops(pf)=ifactors(n)[2];如果nops(pf)=1,那么f:=op(1,pf);p:=op(1,f);e:=op(2,f);如果p=3,那么n;elif p p mod 3=1,然后(p-1)*e+p)*p^(e-1);else p ^(2*地板(e/2));结束if结束,else;else mul(procnname(op(1,p)^op(2,p)),p=pfpf pf pf pf pf pf p=pf(e mod 3,p)3=1,p;结束if;end if;end proc: