110702 Carmichael Numbers


import java.util.*;

import static java.lang.Math.*;

public class Main {
	private static boolean IsPrime(long number)
	{
		long mid = (long)(sqrt((double)number));
		for (long i = 2; i <= mid; ++i)
		{
			if ((number % i) == 0)
				return false;
		}
		
		return true;
	}
	
	private static long CarmichaelMod(long a, long powers, long number)
	{
		if (powers == 1)
			return a % number;
		if (powers == 2)
			return (a * a) % number;
		
		long result = CarmichaelMod(a, powers/2, number);
		result = (result * result) % number;
		
		if ((powers % 2) != 0)
			result = (result * a) % number;
		return result;
	}
	
	private static boolean MeetCarmichaelRule(long number, long a)
	{
		return CarmichaelMod(a, number, number) == a;
	}
	
	private static boolean IsCarmichaelNumber(long number)
	{
		if (IsPrime(number))
			return false;
		
		for (long i = 2; i < number; ++i)
		{
			if (!MeetCarmichaelRule(number, i))
				return false;
		}
		
		return true;
	}
	
	private static boolean IsCarmichaelNumber(HashSet<Long> allCarmichaelNumbers, long number)
	{
		return allCarmichaelNumbers.contains(number);
	}

	private static HashSet<Long> GetAllCarmichaelNumbers(long limit)
	{
		HashSet<Long> result = new HashSet<Long>();
		for (long i = 4; i < limit; ++i)
		{
			if (IsCarmichaelNumber(i))
				result.add(i);
		}
		
		return result;
	}

	public static void main(String[] args) {
		HashSet<Long> allCarmichaelNumbers = GetAllCarmichaelNumbers(65000);
    	Scanner in = new Scanner(System.in);
    	while (true)
    	{
    		long number = in.nextLong();
    		if (number == 0)
    			return;
    		if (IsCarmichaelNumber(allCarmichaelNumbers, number))
    			System.out.println("The number " + number + " is a Carmichael number.");
    		else
    			System.out.println(number + " is normal.");
    	}
	}

}


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

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值