package cn.edu.ylanhds.fabl.mathematical.formula;
/**
* 一些不怎么常用的数学公式的JAVA代码
*
* @author fablen
*
*/
public class CommonDivisor {
public static void main(String [] args) {
CommonDivisor commonDivisor = new CommonDivisor();
System.out.println(commonDivisor.gcd(1125, 855));
System.out.println(commonDivisor.abs(9));
System.out.println(commonDivisor.abs(9.5));
System.out.println(commonDivisor.isPrime(8));
System.out.println(commonDivisor.sqrt(8.0));
System.out.println(commonDivisor.hypotenuse(3.0,4.0));
System.out.println(commonDivisor.H(2));
}
/**
* 计算两个非负整数P和q的最大公约数;若q是0,则最大公约数为p.否则,
* 将p除以q得到余数r,p和q的最大公约数即为q和r的最大公约数
* @param p
* @param q
* @return
*/
public static int gcd(int p, int q) {
if (q == 0) {
return p;
}
int r = p % q;
return gcd(q, r);
}
/**
* 计算一个整数的绝对值
* @param x
* @return
*/
public static int abs(int x) {
if(x<0) {
return -x;
}else {
return x;
}
}
/**
* 计算一个浮点数的绝对值
* @param x
* @return
*/
public static double abs(double x) {
if(x< 0.0) {
return -x;
}else {
return x;
}
}
/**
* 判断一个数是否是素数
* @param N
* @return
*/
public static boolean isPrime(int N) {
if (N < 2) {
return false;
}
for (int i = 2; i*i <=N ; i++) {
if (N % i == 0) {
return false;
}
}
return true;
}
/**
* 计算平方根(牛顿迭代法)
* @param c
* @return
*/
public static double sqrt(double c) {
if (c < 0) {
return Double.NaN;
}
double err = 1e-15;
double t = c;
while (Math.abs(t-c/t) > err * t) {
t = (c/t + t)/2.0;
}
return t;
}
/**
* 计算直角三角形的斜边
* @param a 边长a
* @param b 边长b
* @return
*/
public static double hypotenuse(double a, double b) {
return Math.sqrt(a*a + b*b);
}
/**
* 计算调和级数
* @param N
* @return
*/
public static double H(int N){
double sum = 0.0;
for(int i = 1; i<= N;i++) {
sum += 1.0/i;
}
return sum;
}
}
本文是学习《算法》一书中的代码的实现笔记。