BC开方

function bcroot($幂, $指数 = '3', $scale = 0){
	#开方

	if($幂 == '0'){
		return '0';
	}elseif($幂 == '1'){
		return '1';
	}

	$底扩位 = strval($scale +1);
	$幂扩位 = bcmul($指数,$底扩位);

	$幂 = bcmul($幂,bcpow('10',$幂扩位));

	$小数点位 = strpos($幂,'.');
	if($小数点位 !== FALSE){
		$幂 = substr($幂,0,$小数点位);
	}
	
	$幂长 = strlen($幂);
	$底数长 = ceil($幂长 / intval($指数));

	$底数 = '';
	for($i=0;$i<$底数长;$i++){
		for($j=0;$j<10;$j++){
			$新数 = $底数.strval($j);
			$新数 = str_pad($新数,$底数长,'0');
			$新幂 = bcpow($新数,$指数);
			$比较 = bccomp($幂,$新幂);
			if($比较 == 1){
				#幂大
				if($j == 9){
					$底数 .= '9';
				}
			}elseif($比较 == -1){
				#新幂大
				$底数 .= strval($j -1);
				break;
			}else{
				#相等
				$底数 = $新数;
				break 2;
			}
		}
	}
	
	$底数 = bcdiv($底数,bcpow('10',$底扩位),$scale);
	return $底数;
}

 

转载于:https://my.oschina.net/u/2356887/blog/2248893