青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品

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 閱讀(208) 評論(0)  編輯 收藏 引用 所屬分類: 數論

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

導航

統計

公告

常用鏈接

留言簿(3)

隨筆分類

隨筆檔案

me

好友

同學

網友

搜索

最新評論

閱讀排行榜

評論排行榜

青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品
  • <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>
            夜夜夜精品看看| 午夜精品偷拍| 国产精品视频观看| 国产精品视频久久久| 欧美日韩中文字幕综合视频| 欧美日韩一区三区| 国产欧美日韩亚洲精品| 国产自产在线视频一区| **欧美日韩vr在线| 日韩一区二区免费高清| 午夜视频久久久| 欧美成年人视频网站欧美| 最近看过的日韩成人| 日韩视频一区二区三区在线播放| 亚洲手机成人高清视频| 久久久美女艺术照精彩视频福利播放| 欧美成人影音| 国产一区欧美日韩| 亚洲网站在线| 免费在线亚洲| 亚洲一区二区日本| 免费看成人av| 国产精品毛片a∨一区二区三区|国| 伊人成人在线| 欧美一区免费视频| 亚洲精品乱码久久久久久日本蜜臀 | 欧美成黄导航| 国产精品日韩久久久| 亚洲精选视频在线| 久久久综合激的五月天| 日韩亚洲欧美精品| 免费在线播放第一区高清av| 国产精品青草久久久久福利99| 亚洲激情网站| 久热精品视频在线观看| 亚洲主播在线| 欧美精品一区二区三区高清aⅴ| 国外视频精品毛片| 欧美一区免费| 中文一区二区在线观看| 欧美日韩精品免费| 亚洲毛片在线观看.| 麻豆精品一区二区av白丝在线| 亚洲一区二区三区四区五区午夜 | 国产酒店精品激情| 亚洲日韩成人| 欧美v亚洲v综合ⅴ国产v| 欧美一区网站| 国产日产亚洲精品| 欧美中文字幕在线视频| 亚洲一区二区三区四区五区午夜| 欧美日韩成人网| 亚洲日本乱码在线观看| 欧美激情在线| 欧美成人一区二区三区在线观看 | 欧美日韩人人澡狠狠躁视频| 亚洲黄一区二区| 欧美高清视频在线观看| 鲁大师成人一区二区三区| 亚洲第一综合天堂另类专| 毛片精品免费在线观看| 久久理论片午夜琪琪电影网| 亚洲成人在线视频网站| 欧美成人综合| 欧美国产三区| 亚洲天堂av在线免费| 中文精品视频一区二区在线观看| 国产精品区二区三区日本| 久久国产精品高清| 久久乐国产精品| 亚洲精品偷拍| 亚洲午夜成aⅴ人片| 国产农村妇女毛片精品久久莱园子| 欧美一级黄色录像| 久久久精品国产99久久精品芒果| 亚洲国产高清自拍| 亚洲国产欧美在线人成| 欧美视频二区36p| 欧美一区二粉嫩精品国产一线天| 欧美一区日韩一区| 亚洲精品乱码久久久久久| 在线一区免费观看| 国精产品99永久一区一区| 欧美+亚洲+精品+三区| 欧美黑人国产人伦爽爽爽| 亚洲自啪免费| 久久综合色综合88| 亚洲视屏一区| 久久久蜜桃精品| 亚洲视屏在线播放| 久久激情综合网| 在线综合亚洲欧美在线视频| 欧美在线free| 一区二区三区欧美激情| 欧美在线观看网址综合| 日韩午夜在线| 久久成人精品无人区| 在线性视频日韩欧美| 久久动漫亚洲| 鲁大师影院一区二区三区| a4yy欧美一区二区三区| 亚洲一区三区视频在线观看| 在线观看日韩一区| 亚洲视频网站在线观看| 亚洲黄色免费网站| 亚洲综合精品一区二区| 亚洲区第一页| 欧美在线综合视频| 亚洲视频碰碰| 模特精品在线| 美女免费视频一区| 国产免费亚洲高清| 亚洲精品美女久久7777777| 狠狠色香婷婷久久亚洲精品| 中文网丁香综合网| 日韩一级黄色大片| 欧美成人性网| 欧美风情在线观看| 在线成人免费观看| 欧美中文字幕久久| 欧美一区二区久久久| 欧美日韩一区二区在线| 亚洲精品国产精品国自产观看| 在线播放国产一区中文字幕剧情欧美 | 欧美国产视频在线观看| 韩日精品视频一区| 欧美一区二区视频97| 篠田优中文在线播放第一区| 欧美三级午夜理伦三级中视频| 亚洲国产成人一区| 亚洲人成网站精品片在线观看| 久久亚洲电影| 欧美国产精品日韩| 亚洲国产欧美在线人成| 久久亚洲精选| 欧美激情第9页| 99精品欧美一区二区三区| 欧美日产在线观看| 一区二区三区四区五区精品视频| 亚洲一区二区在线免费观看视频| 欧美色精品天天在线观看视频| 一本久久青青| 久久成年人视频| 国内外成人免费激情在线视频| 久久精品导航| 欧美黄免费看| 制服诱惑一区二区| 国产精品视频网址| 久久久免费精品视频| 欧美插天视频在线播放| 日韩一级大片在线| 国产精品日本精品| 欧美专区福利在线| 欧美激情一区二区三区高清视频| 一本色道久久综合亚洲二区三区 | 欧美国产精品va在线观看| 亚洲精品视频中文字幕| 亚洲小说区图片区| 国产欧美视频一区二区| 久久婷婷久久一区二区三区| 91久久精品国产| 久久精品国产一区二区电影 | 亚洲精品久久久久久一区二区 | 午夜欧美精品| 免费看黄裸体一级大秀欧美| 亚洲久久视频| 国产视频在线一区二区| 久久人人爽人人爽爽久久| 亚洲人成久久| 久久久久国产精品麻豆ai换脸| 亚洲黄色成人| 国产精品揄拍一区二区| 乱人伦精品视频在线观看| 亚洲图中文字幕| 欧美 亚欧 日韩视频在线| 亚洲欧美激情在线视频| 亚洲第一福利在线观看| 欧美日韩一区在线观看| 久久精品视频在线| 在线亚洲一区| 亚洲韩国日本中文字幕| 乱中年女人伦av一区二区| 亚洲一区高清| 亚洲精品国产精品乱码不99按摩| 国产欧美日韩三级| 欧美日韩综合在线| 欧美高清在线一区| 久久久在线视频| 欧美在线播放一区| 亚洲女性喷水在线观看一区| 亚洲久久一区二区| 亚洲第一区色| 男女激情视频一区| 久久精品理论片| 欧美一区二区黄色| 亚洲女爱视频在线| 亚洲午夜精品一区二区三区他趣 | av不卡免费看| 亚洲伦理自拍| 日韩午夜精品|