def hanleyMcNeil(auc: float, n0: int, n1: int) -> float:
	"""The very good power-law variance estimate from Hanley/McNeil"""
	auc2=auc*auc
	q1=auc/(2.-auc)
	q2=2.*auc2/(1.+auc)
	return( (auc-auc2+(n1-1.)*(q1-auc2)+(n0-1.)*(q2-auc2))/n0/n1 )