圆周率 π 的求法(二)

圆周率 π 的求法(二)

一个胖子

昨天写了个 demo 求圆周率,后来想想发现自己的方法有点冗余现在来重新构建下方法.

大体思路并没有改变,依然是通过扇形面积来求解圆周率 π .但是对扇形的面积求解方式进行了调整.

在上篇文章中,我们采用的是微元分割的方式对整个图形进行了分割,对所有单元的面积进行叠加求和,最后得出的扇形面积的近似值.

其实我们不难发现,我们在求解面积时,我们关注的仅仅是边界值,即在边界位置的(x,y)值.

由图像我们不难得出结论:随着随着 x 增大, y 值不断减小.所以边界方格整体向右下方移动.

因此我们在求解其边界时可以不用逐一遍历,只需在上一边界值的右下方向上的方格进行寻找即可.

例如 x=0.1 是 方格否边界 y =0.9 .当x = 0.2 时 y值直接从 0.9 开始向下寻找即可,这样计算量将会被大大的减少.时间复杂度从原来的 n² 降低为 n.

代码实现如下

 public static void main(String[] args) {
    System.out.println("请输入你想划分的网格数:");
    Scanner sc = new Scanner(System.in);
    //获取每条边被分成多少格
    int num = sc.nextInt();
    //声明扇形中所有单元格个数
    int sum = 0;
    //初始化第一次运算时的上一个 y ;
    int lastY = num -1 ;
    for (int x = 0; x < num; x++) {

        for(int y = lastY; y >= 0; y--){
            double newX = (2*x+1)*1.0/(2*num);
            double newY = (2*y+1)*1.0/(2*num);
            //判断该节点是否在扇形中
            if((newX*newX+newY*newY)<=1){
                lastY = y ;
                break ;
            }
        }
        //lastY为当前第 n-1 个方格因此在最后时应加 1 ;
        sum += lastY + 1;
    }
    //通过扇形中单元个数和总单元个数的比值来确定 圆周率 π 的大小
    System.out.println("圆周率π的大小为 π = "+ (sum*4*1.0/(num*num)));
}

---

下面的方式我并不想写的

上述方式已经很类似于 对函数 y = √￣(1 - x²) 求解的过程

因此整个面积的求解的过程可以变化为积分的求解的过程.

代码如下

 public static void main(String[] args) {
    System.out.println("请输入你想划分的网格数:");
    Scanner sc = new Scanner(System.in);
    //获取每条边被分成多少格
    int num = sc.nextInt();
    //声明扇形中所有单元格个数
    int sum = 0;
    for (int x = 0; x < num; x++) {
        double newX = (2*x+1)*1.0/(2*num);
        sum += Math.sqrt(1-newX*newX);
    }
    //通过扇形中单元个数和总单元个数的比值来确定 圆周率 π 的大小
    System.out.println("圆周率π的大小为 π = "+ (sum*4*1.0/(num*num)));
}

我不想用这种方法是因为感觉很蠢.

因为调用了Math提供的方法,而Math中直接提供了圆周率的数值

Math.PI=3.14159265358979323846

这种精度的数据我不知道需要多精细的划分才能达到.

未完待续

(最近在学大数据分布式计算,明天试着将任务拆分扔给集群)