
/**  Prime Factorial Decomposition **/
/*A: K McIsaac*/
/*S: University of Western Australia, Nedlands 6009, Australia*/
/*D: November 1985*/
/*K: Nfac, Prime, Factors */

/*: Pfac[$num]
	Gives the prime factor decomposition of $num in a 
	2-D format*/

_Pfacp[Pr][$a,$b] :: Fmt[{{0,0},{1,1}},$a,$b]
_Pfacp[Pr][0,$a] : 0
_Pfacp[Pr][1,$a] : 1
_Pfacp[Pr][$a,1] : $a

Pfac_:Tier;Pfac_:Ldist;
Pfac[$a _= Sign[$a]= -1] :: Pfac[-1] Pfac[-$a]
Pfac[$a] :: Ap[' Mult,Map[Ap[Pfacp,$1],Nfac[ $a]]]

Nfac[0]  : {{0,1}}
Nfac[1]  : {{1,1}}
Bigfac[B[0]]  : {{0,1}}
Bigfac[B[1]]  : {{1,1}}
Pfac[-1] : Pfacp[-1,1]

/*E:
SMP 1.5.0
Fri Nov  8 11:31:05 1985


#I[1]::  <XPfac

#I[2]::  Pfac[60]

	  2
#O[2]:*  2  3 5

#I[3]::  Pfac[-60]

	     2
#O[3]:*  -1 2  3 5

#I[4]::  Pfac[0]

#O[4]:*  0

#I[5]::  Pfac[6]

#O[5]:*  2 3

#I[6]::  <end>
*/
