• <ins id="pjuwb"></ins>
    <blockquote id="pjuwb"><pre id="pjuwb"></pre></blockquote>
    <noscript id="pjuwb"></noscript>
          <sup id="pjuwb"><pre id="pjuwb"></pre></sup>
            <dd id="pjuwb"></dd>
            <abbr id="pjuwb"></abbr>

            poj 1811 Prime Test 數論 素數測試

             代碼如下:
            #include <stdio.h>
            #include <time.h>
            #include <math.h>
            #include <stdlib.h>
            #include <algorithm>
            using namespace std;
            typedef unsigned long long LL;
            #define MAX (5000000)
            bool bPrime[MAX];
            void InitPrime()
            {
                int nMax = sqrt((double)MAX) + 1;
                bPrime[0] = bPrime[1] = true;
                for (int i = 2; i <= nMax; ++i)
                {
                    if (!bPrime[i])
                    {
                        for (int j = 2 * i; j < MAX; j += i)
                        {
                            bPrime[j] = true;
                        }
                    }
                }
            }
            LL multAndMod(LL a, LL b, LL n)
            {
                LL tmp = 0;
                while (b)
                {
                    if(b & 1)
                    {
                        tmp = (tmp + a) % n;
                    }
                    a = (a << 1) % n;
                    b >>= 1;
                }
                return tmp;
            }
            //計算a^u%n
            LL ModExp(LL a, LL u, LL n)
            {
                LL d = 1;
                a %= n;
                while (u)
                {
                    if (u & 1)
                    {
                        d = multAndMod(d, a, n);
                    }
                    a = multAndMod(a, a, n);
                    u >>= 1;
                }
                return d % n;
            }
            //判斷nN是不是合數
            bool Witness(LL a, LL nN)
            {
                LL u = nN - 1, t = 0;//將nN-1表示為u*2^t
                while (u % 2 == 0)
                {
                    t++;
                    u >>= 1;
                }
                LL x0 = ModExp(a, u, nN);//x是a^u
                LL x1;
                for (int i = 1; i <= t; ++i)
                {
                    x1 = multAndMod(x0, x0, nN);
                    if (x1 == 1 && x0 != nN - 1 && x0 != 1)
                    {
                        return true;
                    }
                    x0 = x1;
                }
                if (x1 != 1)
                {
                    return true;
                }
                return false;
            }
            //素數測試
            bool MillerRabin(LL nN)
            {
                //if (nN < MAX)return !bPrime[nN];
                const int TIME = 10;
                for (int i = 0; i < TIME; ++i)
                {
                    LL a = rand() % (nN - 1) + 1;
                    if (Witness(a, nN))
                    {
                        return false;
                    }
                }
                return true;
            }
            LL gcd(LL a, LL b)
            {
                if (a < b)swap(a, b);
                while (b)
                {
                    LL t = a;
                    a = b;
                    b = t % b;
                }
                return a;
            }
            //啟發式尋找nN的因子
            LL PollardRho(LL n, LL c)
            {
                LL i = 1, t = 2;
                LL x, y;
                LL ans;
                srand(time(NULL));  
                y = x = rand() % n;
                while(1)
                {
                    i++;
                    x = (multAndMod(x, x, n) + c) % n;
                    ans = gcd(y - x, n);
                    if(ans > 1 && ans < n)
                        return ans;
                    if(x == y)
                        return n;
                    if(t == i)
                    {
                        y = x;
                        t <<= 1;
                    }
                }
            }
            LL FindMin(LL nN, LL c)
            {
                //printf("nN:%I64u\n", nN);
                if (MillerRabin(nN) || nN <= 1)
                {
                    return nN;
                }
                LL p = nN;
                while (p >= nN) p = PollardRho(p, c--);
                if (p > 1)
                    p = FindMin(p, c);//分解p的最小因子
                if (p < nN)
                {
                    LL q = nN / p;
                    q = FindMin(q, c);//找到q的最小因子
                    p = min(p, q);
                }
                return p;
            }
            int main()
            {
                int nTest;
                srand(time(NULL));
                //InitPrime();
                scanf("%d", &nTest);
                while (nTest--)
                {
                    LL nN;
                    scanf("%I64u", &nN);
                    if (nN > 2 && nN % 2 == 0)
                    {
                        printf("2\n");
                    }
                    else if (nN == 2 || MillerRabin(nN))
                    {
                        printf("Prime\n");
                    }
                    else
                    {
                        printf("%I64u\n", FindMin(nN, 181));
                    }
                }
                return 0;
            }

            posted on 2012-09-24 12:56 yx 閱讀(188) 評論(0)  編輯 收藏 引用 所屬分類: 數論

            <2011年11月>
            303112345
            6789101112
            13141516171819
            20212223242526
            27282930123
            45678910

            導航

            統計

            公告

            常用鏈接

            留言簿(3)

            隨筆分類

            隨筆檔案

            me

            好友

            同學

            網友

            搜索

            最新評論

            閱讀排行榜

            評論排行榜

            久久精品国产精品亚洲人人| 久久久WWW成人免费精品| 久久人人爽人人爽人人片AV高清| 色妞色综合久久夜夜| 97久久精品午夜一区二区| 99久久人人爽亚洲精品美女| 久久久久亚洲精品日久生情| 久久久久久国产精品无码超碰| 精品乱码久久久久久夜夜嗨| 亚洲中文字幕无码一久久区| 99久久亚洲综合精品网站| 狠狠色婷婷久久综合频道日韩| 国产—久久香蕉国产线看观看| 久久久www免费人成精品| 国产 亚洲 欧美 另类 久久| 久久精品中文騷妇女内射| 久久亚洲精品无码观看不卡| 伊人色综合久久天天| 久久久噜噜噜久久中文福利| 久久99国产精品久久99小说| 91麻豆精品国产91久久久久久| 国产精品久久久久aaaa| 日本欧美久久久久免费播放网| 欧美激情精品久久久久久| 一级做a爰片久久毛片人呢| 国产婷婷成人久久Av免费高清| 久久精品国产日本波多野结衣| 亚洲Av无码国产情品久久| 久久精品国产只有精品66| 久久www免费人成精品香蕉| 久久免费视频网站| 久久精品国产亚洲av日韩| 久久棈精品久久久久久噜噜| 亚洲国产精品久久久天堂| 伊人久久大香线蕉综合影院首页| 久久综合精品国产一区二区三区| 曰曰摸天天摸人人看久久久| 久久不见久久见免费影院www日本| 精品无码久久久久久久动漫| 人人狠狠综合久久亚洲高清| 亚洲精品午夜国产va久久|