7-35 蒙特卡罗方法求圆周率

最新推荐文章于 2023-09-23 21:44:03 发布

懒懒猪～。～

最新推荐文章于 2023-09-23 21:44:03 发布

阅读量137

点赞数

文章标签： java 算法开发语言

本文链接：https://blog.csdn.net/qq_74222774/article/details/129960631

版权

题目描述：

使用蒙特卡洛仿真方法求圆周率。（具体要求见作业指导书
2020-OO第05次作业-2指导书V1.0.pdf
）

输入格式:

从键盘输入四个实型数和一个整型数，分别为矩形左上角的横坐标、纵坐标、矩形长度、矩形宽度和投点次数，数与数之间可以用一个或多个空格或回车分隔。

输出格式:

如果矩形长度与宽度不相等（非正方形）或长宽数据非法，则输出“Wrong Format”。
如果估算出的π与Math.PI差值小于1E-4，则输出“Success”，否则输出“failed”。

输入样例：

在这里给出一组输入。例如：

0 0 1 1 20000000

输出样例：

在这里给出相应的输出。例如：

Success

答案：

import java.util.Random;
import java.util.Scanner;

public class Main {

public static void main(String[] args) {
    Scanner sc = new Scanner(System.in);
    Random r = new Random();
    double x = sc.nextDouble();
    double y = sc.nextDouble();
    double width = sc.nextDouble();
    double height = sc.nextDouble();
    int n = sc.nextInt();
    if (width == height) {
        double ox = x + width / 2.0;  //圆心x
        double oy = y - height / 2.0;  //圆心y
        double or = x - ox;
        int m = 0;
        for (int i = 0; i < n; i++) {
            double ranx = r.nextDouble() * width+x;
            double rany = r.nextDouble() * height+y-height;
            if (Math.pow(ranx - ox, 2) + Math.pow(rany - oy, 2) <= Math.pow(or, 2)) {
                m++;
            }
        }
        double pai = m * 4.0 / n;
        if (Math.abs(pai - Math.PI) <= 1e-3) {
            System.out.println("Success");
        } else {
            System.out.println("failed");
        }
    } else {
        System.out.println("Wrong Format");
    }
}

}