常用小算法

标签：

1.快速模取幂

要求 a^b % c的时候，当然可以直接循环然后每一步都mod c，这样的复杂度就是O(b)了，但是快速模取幂算法可以降低复杂度。

描述如下：

   数论计算中经常出现的一种运算就是求一个数的幂a^b对另外一个数n个模的运算，即计算:

a^b mod n (a,b,n是正整数)

    由于计算机只能表示有限位的整数，所以编程时模取幂的运算要注意值的大小范围，当a^b的值超过整数范围时,mod运算便无法进行。

    如何解决这个问题，我们引出一个能计算a^b mod n的值的有用算法——反复平方法，首先我们必须明确：

d=a^b mod n=(…((((a mod n)*a)mod n)*a)mod n…*a)mod n    {共b个a}

    由此可以引出一个迭代式

         d:=a;

         for i:=2 to b do

            d:=d mod n*a;

         d:=d mod n;

    时间复杂度为O（b），当b很大时，效率很低。我们可以将b转换为二进制数<b_k,b_k-1,...,b₁,b₀>,然后从最低位b₀开始，由右至左逐位扫描，每次迭代时，用到下面两个恒等式：

a^2c mod n =(a^c)² mod n      b_i=0            

a^2c+1 mod n =a*(a^c)² mod n   b_i=1 (0<=c<=b)

    其中c为b的二进制数的后缀（b_i-1...b₀）对应的十进制数，当c成倍增加时，算法保持d=a^c mod n不变，直至c=b。

    //快速模取幂算法
    private static long mod_pow(long x, long n, long mod) {
        long res = 1;
        while(n > 0) {
            if((n & 1) != 0) res = res * x % mod;
            x = x * x % mod;
            n >>= 1;
        }
        return res;
    }

2.判断素数

    public boolean isPrime(int n) {
        if (n < 2)
            return false;
        int k = (int)Math.sqrt(n);
        for (int i = 2; i <= k; i++) {
            if (n % i == 0)
                return false;
        }
        return true;
    }

常用小算法

标签：

原文地址：http://www.cnblogs.com/fisherinbox/p/5497491.html