ZOJ-1337

import java.util.Scanner;
import java.util.TreeSet;

public class Main
{
	static TreeSet<Integer> primes = new TreeSet<Integer>();

	public static void main(String[] args)
	{
		generatePrimes();

		Scanner sc = new Scanner(System.in);

		while (true)
		{
			int n = sc.nextInt();

			if (n == 0)
				break;

			int[] array = new int[n];

			int total = 0;
			int common = 0;

			for (int i = 0; i < n; i++)
			{
				array[i] = sc.nextInt();
				for (int j = 0; j < i; j++)
				{
					total++;

					if (hasCommonFactors(array[i], array[j]))
						common++;
				}
			}

			if (total - common == 0)
				System.out.println("No estimate for this data set.");
			else
				System.out.format("%.6f\n", Math.sqrt(6.0 * total / (total - common)));
		}

	}

	static boolean hasCommonFactors(int a, int b)
	{
		if (a == 1 || b == 1)
			return false;
		if (a == 0 || b == 0)
			return true;
		if (a == b)
			return true;

		int max = a > b ? a : b;
		int min = a + b - max;

		if (primes.contains(min))
		{
			if (primes.contains(max))
				return false;
			else if (max % min == 0)
				return true;
			return false;
		}
		else if (primes.contains(max))
			return false;
		else
		{
			for (Integer n : primes)
			{
				if (n > min)
					return false;

				if (min % n == 0 && max % n == 0)
					return true;
			}
			return false;
		}
	}

	static void generatePrimes()
	{
		primes.add(2);
		primes.add(5);

		for (int i = 3; i < 32768; i += 2)
		{
			int tail = i % 10;

			if (tail == 1 || tail == 3 || tail == 7 || tail == 9)
			{
				boolean isPrime = true;
				for (Integer n : primes)
				{
					if (i % n == 0)
					{
						isPrime = false;
						break;
					}
				}

				if (isPrime)
					primes.add(i);
			}
		}
	}
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值