GTY's math problem （BestCoder Round #29 A）

题目大意：
给出四个数a,b,c,d,(1<=a,b,c,d<=1000)，比较a^b与c^d的大小；

思路分析：
直接比较b*lna与d*lnc的大小就行；

调用<math.h>中的log, double log(double x);

注意浮点数精度比较问题，浮点数不能判断相等问题，此时要用eps，大小是可以直接比较的。

 **eps**

eps缩写自epsilon，表示一个小量，但这个小量又要确保远大于浮点运算结果的不确定量。eps最常见的取值是1e-8左右。引入eps后，我们判断两浮点数a、b相等的方式如下:

定义三出口函数如下: int sgn(double a){return a < -eps ? -1 : a < eps ? 0 : 1;}

则各种判断大小的运算都应做如下修正:

传统意义	修正写法1	修正写法2
a == b	sgn(a - b) == 0	fabs(a – b) < eps
a != b	sgn(a - b) != 0	fabs(a – b) > eps
a < b	sgn(a - b) < 0	a – b < -eps
a <= b	sgn(a - b) <= 0	a – b < eps
a > b	sgn(a - b) > 0	a – b > eps
a >= b	sgn(a - b) >= 0	a – b > -eps

这样，我们才能把相差非常近的浮点数判为相等;同时把确实相差较大(差值大于eps)的数判为不相等。

代码实现：

#include<stdio.h>
#include<math.h>
const double eps = 1e-10;

int dcmp(double x) {
  if(fabs(x)<eps) return 0;
  else return x<0?-1:1;
}

int main(){
    double a,b,c,d;
    while(~scanf("%lf%lf%lf%lf",&a,&b,&c,&d)){
        double t1=b*log(a);
        double t2=d*log(c);
        int ans=dcmp(t1-t2);
        if(ans==1) printf(">\n");
        else if(ans==0) printf("=\n");
        else printf("<\n");
    }
}

用Java大数快速幂超时代码：

import java.math.*;
import java.util.*;

public class Main {

	public static void main(String[] args) {
		Scanner in = new Scanner(System.in);
		while(in.hasNext()){
			BigInteger a=in.nextBigInteger();
			int b=in.nextInt();
			BigInteger c=in.nextBigInteger();
			int d=in.nextInt();
			BigInteger t1=new BigInteger("1");
			BigInteger t2=new BigInteger("1");
			while(b!=0){
	    		if((b&1)==1) t1=t1.multiply(a);
	    		a=a.multiply(a);
	    		b>>=1;
	    	}
			while(d!=0){
	    		if((d&1)==1) t2=t2.multiply(c);
	    		c=c.multiply(c);
	    		d>>=1;
	    	}
			if(t1.compareTo(t2)>0) System.out.println(">");
			else if(t1.compareTo(t2)==0) System.out.println("=");
			else System.out.println("<");
				
		}
	}
}