def smoothness(n):
 if (n == 1):
  return (1, 1)
 facs = factorint(n)
 return (max(facs), max(((m ** facs[m]) for m in facs)))

def smoothness ( n ) : if ( n == 1 ) : return ( 1 , 1 ) facs = factorint ( n ) return ( max ( facs ) , max ( ( ( m ** facs [ m ] ) for m in facs ) ) )
return the b-smooth and b-power smooth values of n .