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

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;
}
//啟發(fā)式尋找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 閱讀(207) 評論(0)  編輯 收藏 引用 所屬分類: 數論

<2025年10月>
2829301234
567891011
12131415161718
19202122232425
2627282930311
2345678

導航

統(tǒng)計

公告

常用鏈接

留言簿(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>
            亚洲色图自拍| 亚洲精品一区二区在线观看| 久久精品动漫| 欧美在线亚洲| 欧美在线免费观看视频| 久久成人18免费网站| 欧美在线视频一区二区| 久久久99国产精品免费| 久久中文字幕一区| 欧美国产免费| 欧美午夜久久| 国产性做久久久久久| 在线观看亚洲| 一本久久综合亚洲鲁鲁五月天| 一本一本久久a久久精品综合麻豆| 亚洲网站在线看| 久久成人在线| 久久亚洲国产精品一区二区| 久久综合中文字幕| 欧美成人a视频| 91久久精品国产91久久性色tv| 欧美国产高清| 一区二区三区日韩精品视频| 亚洲韩国一区二区三区| 欧美日韩一区二区三区免费| 久久久久se| 免费成人av资源网| 欧美精品成人| 国产精品视频999| 国产日韩在线看片| 国产一区二区三区久久悠悠色av| 国产亚洲永久域名| 亚洲电影免费观看高清| 亚洲免费电影在线| 欧美一区二区三区免费观看| 久久露脸国产精品| 欧美激情亚洲另类| 一区二区三区欧美亚洲| 午夜精品久久| 欧美aa国产视频| 欧美天堂亚洲电影院在线观看| 欧美亚洲一区在线| 久久精品在线观看| 亚洲国产日本| 亚洲视频精品| 亚洲欧洲日产国码二区| 亚洲欧美一区二区三区久久| 久久国产精彩视频| 欧美高清视频| 亚洲一区三区电影在线观看| 久久本道综合色狠狠五月| 欧美成人精品一区二区三区| 国产精品久久久久久久久动漫 | 欧美日韩国产色综合一二三四| 欧美日韩亚洲高清一区二区| 红桃视频国产一区| 中文国产成人精品| 美女脱光内衣内裤视频久久网站| 一区二区欧美激情| 中文在线不卡视频| 麻豆成人在线| 亚洲欧美在线磁力| 欧美性大战久久久久久久蜜臀| 亚洲电影免费在线 | 亚洲电影在线免费观看| 午夜精品久久久久久久白皮肤| 欧美激情一区二区三区成人| 欧美伊人久久久久久久久影院 | 中日韩美女免费视频网站在线观看| 久久精品首页| 宅男噜噜噜66一区二区| 欧美精品一区二区在线播放| 国产日韩欧美综合精品| 亚洲在线免费视频| 亚洲精品一二区| 欧美成人精品不卡视频在线观看| 韩国久久久久| 久久亚洲图片| 久久精品国产精品亚洲综合| 国产午夜精品一区二区三区欧美| 午夜精品一区二区三区四区| 99成人在线| 欧美日本一区二区三区 | 国产精品色在线| 亚洲影视在线| 亚洲私人影院在线观看| 欧美日韩一区二区精品| 日韩视频一区二区三区在线播放| 亚洲电影在线免费观看| 免费成人高清视频| 亚洲激情女人| 日韩一区二区精品视频| 欧美日韩国产首页| 在线亚洲欧美视频| 一本色道精品久久一区二区三区 | 欧美日韩999| 亚洲精品美女免费| 亚洲欧洲在线观看| 欧美日韩国产丝袜另类| 亚洲欧美国产制服动漫| 亚洲福利视频专区| 久久免费视频一区| 亚洲毛片在线免费观看| 一区二区毛片| 国产一区二区三区直播精品电影| 美女在线一区二区| 欧美国产视频一区二区| 亚洲欧美一区二区在线观看| 久久成人18免费观看| 在线看不卡av| 亚洲美女中文字幕| 国产午夜久久| 亚洲国产欧美一区二区三区同亚洲 | 亚洲淫性视频| 久久www免费人成看片高清| 亚洲国产精品久久久久婷婷884| 亚洲精品韩国| 国内精品伊人久久久久av影院 | 亚洲黑丝在线| 国产女优一区| 亚洲三级免费观看| 国产日产欧美a一级在线| 亚洲国产精彩中文乱码av在线播放| 欧美久久久久久久久| 久久久精品午夜少妇| 欧美日韩不卡| 毛片一区二区| 国产伦精品一区| 亚洲电影免费观看高清| 国产亚洲成年网址在线观看| 亚洲激情中文1区| 激情自拍一区| 中文精品视频| 日韩视频在线你懂得| 欧美专区一区二区三区| 亚洲桃花岛网站| 美女免费视频一区| 久久久噜噜噜久久| 国产精品久久久久久妇女6080 | 一本色道久久综合精品竹菊| 尤物网精品视频| 欧美在线观看一区二区三区| 亚洲一区美女视频在线观看免费| 免费91麻豆精品国产自产在线观看| 欧美一区永久视频免费观看| 欧美日韩在线一二三| 亚洲黄页视频免费观看| 亚洲第一色在线| 欧美专区一区二区三区| 欧美成人午夜| 香蕉av777xxx色综合一区| 欧美不卡在线| 久久夜色精品国产欧美乱极品| 国产精品网红福利| 亚洲一区二区三区欧美| 亚洲欧美日韩精品| 国产精品啊啊啊| 日韩一区二区精品在线观看| 99热在这里有精品免费| 美女精品在线观看| 欧美成人午夜激情在线| 在线观看国产精品淫| 久久青青草综合| 欧美国产另类| 日韩视频免费观看高清完整版| 欧美国产精品中文字幕| 亚洲精品视频在线观看免费| 国产一区二区福利| 亚洲欧美日韩区| 欧美中文字幕在线播放| 国产午夜精品久久久| 久久久九九九九| 亚洲国产成人av好男人在线观看| 亚洲免费成人| 欧美午夜www高清视频| 一区二区三区免费网站| 欧美主播一区二区三区| 久久一区二区精品| 亚洲激情一区二区三区| 欧美激情一区二区三区不卡| 亚洲精品乱码久久久久久日本蜜臀| 在线视频亚洲欧美| 国产精品乱子久久久久| 亚洲欧美视频一区| 久久综合影音| 一本久道久久综合婷婷鲸鱼| 国产精品视频一区二区三区| 久久久久久久波多野高潮日日| 欧美电影在线观看| 亚洲一级在线| 一区免费观看| 欧美日韩 国产精品| 亚洲国产日韩在线一区模特| 99riav1国产精品视频| 国产日韩久久| 欧美激情精品久久久久久免费印度| 一区二区三区色| 欧美大片91| 欧美在线三区| 一区二区国产精品|