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

poj 1811 Prime Test 數(shù)論 素?cái)?shù)測(cè)試

 代碼如下:
#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;
}
//計(jì)算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是不是合數(shù)
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;
}
//素?cái)?shù)測(cè)試
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 閱讀(213) 評(píng)論(0)  編輯 收藏 引用 所屬分類: 數(shù)論

<2012年9月>
2627282930311
2345678
9101112131415
16171819202122
23242526272829
30123456

導(dǎo)航

統(tǒng)計(jì)

公告

常用鏈接

留言簿(3)

隨筆分類

隨筆檔案

me

好友

同學(xué)

網(wǎng)友

搜索

最新評(píng)論

閱讀排行榜

評(píng)論排行榜

青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品
  • <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>
            欧美日韩在线视频观看| 91久久精品日日躁夜夜躁欧美| 欧美黄色精品| 国产伦精品一区二区三区高清| 欧美激情中文字幕乱码免费| 国产精品永久免费在线| 99在线视频精品| 最新国产の精品合集bt伙计| 欧美一区综合| 久久精品国产99国产精品| 欧美日韩在线电影| 91久久久在线| 亚洲国产成人精品女人久久久| 制服丝袜激情欧洲亚洲| 在线视频免费在线观看一区二区| 免费一区视频| 欧美夫妇交换俱乐部在线观看| 国产综合婷婷| 欧美一级片在线播放| 午夜精品婷婷| 国产精品欧美一区二区三区奶水 | 亚洲国产一区二区在线| 狠狠色综合播放一区二区| 亚洲欧美综合国产精品一区| 亚洲欧美日韩国产综合精品二区| 欧美精品三区| 亚洲精品在线视频| 在线视频日韩| 欧美手机在线视频| 亚洲午夜视频在线| 欧美在线亚洲在线| 国产精品中文字幕欧美| 亚洲免费中文| 久久噜噜亚洲综合| 伊人久久大香线| 欧美96在线丨欧| 亚洲国产日韩精品| 亚洲最新视频在线| 国产精品久久久久久久7电影 | 99精品福利视频| 亚洲欧美日韩精品在线| 国产精品家庭影院| 欧美一区二区高清在线观看| 久久免费的精品国产v∧| 精品动漫3d一区二区三区| 久久久青草婷婷精品综合日韩| 美女被久久久| 亚洲伦理在线观看| 国产精品福利网站| 午夜精品免费视频| 欧美a级片网站| 99精品久久免费看蜜臀剧情介绍| 欧美精品日韩三级| 亚洲性人人天天夜夜摸| 欧美一区1区三区3区公司| 狠狠色丁香久久婷婷综合丁香| 免费观看日韩| 亚洲视频欧美在线| 久热精品在线| 一区二区三区日韩| 国产一区二区三区网站| 欧美成人午夜| 亚洲欧美成人综合| 欧美国产日韩免费| 亚洲欧美日韩精品久久亚洲区| 在线播放国产一区中文字幕剧情欧美| 欧美国产视频一区二区| 亚洲综合精品自拍| 欧美福利电影网| 欧美一区二区免费视频| 亚洲国产欧美一区二区三区同亚洲 | 亚洲精品久久久久久久久久久久 | 亚洲一区二区三区高清 | 一二三四社区欧美黄| 国产日韩一区| 欧美高清在线一区二区| 欧美一区二区三区免费看| 亚洲精品乱码久久久久久蜜桃91| 久久成人资源| 亚洲素人一区二区| 亚洲激情第一区| 国产日韩一区二区三区在线播放| 欧美喷水视频| 麻豆成人在线观看| 性伦欧美刺激片在线观看| 日韩视频一区二区三区| 欧美国产日本| 久久综合999| 亚洲欧美日韩一区二区在线| 亚洲精品久久久一区二区三区| 国产精品久久97| 欧美日韩午夜在线| 欧美搞黄网站| 美日韩精品视频| 久久久精品一品道一区| 午夜在线不卡| 亚洲性线免费观看视频成熟| 99精品国产高清一区二区| 亚洲国产精品传媒在线观看| 久久亚洲视频| 久久一区二区三区国产精品| 欧美在线91| 欧美一区2区三区4区公司二百| 亚洲综合成人婷婷小说| 中文在线一区| 在线亚洲+欧美+日本专区| 亚洲精品久久久久中文字幕欢迎你| 在线电影院国产精品| 一区在线电影| 在线看片日韩| 在线观看亚洲精品视频| 激情久久久久久| 一区二区视频免费在线观看| 在线播放豆国产99亚洲| 激情欧美国产欧美| 伊人成年综合电影网| 在线日本成人| 最新国产成人在线观看| 99re6热只有精品免费观看| 亚洲激情视频在线观看| 日韩午夜av在线| 一区二区三区欧美激情| 亚洲女人天堂成人av在线| 亚洲欧美中文字幕| 久久成人18免费网站| 免费亚洲视频| 亚洲高清av在线| 亚洲美女在线国产| 亚洲午夜黄色| 久久aⅴ国产欧美74aaa| 老司机午夜精品视频在线观看| 农村妇女精品| 国产精品v欧美精品∨日韩| 国产精品入口麻豆原神| 伊人成人网在线看| 亚洲精品乱码| 亚洲欧美乱综合| 欧美一区二区黄色| 蜜桃av一区二区| 亚洲人成绝费网站色www| 亚洲私人黄色宅男| 久久九九全国免费精品观看| 男人插女人欧美| 欧美午夜三级| 很黄很黄激情成人| 亚洲美女视频在线观看| 香港久久久电影| 欧美www视频在线观看| 亚洲精品中文字幕有码专区| 亚洲制服丝袜在线| 裸体女人亚洲精品一区| 国产精品大片wwwwww| 永久555www成人免费| 亚洲一区二区在线免费观看视频 | 亚洲欧美清纯在线制服| 免费成人黄色| 亚洲一级电影| 蜜臀av一级做a爰片久久| 国产精品爽爽ⅴa在线观看| 亚洲高清网站| 久久9热精品视频| 91久久视频| 久久久久久久综合| 国产精品中文在线| 夜夜爽www精品| 玖玖在线精品| 亚洲一区二区三区中文字幕| 猛男gaygay欧美视频| 国产欧美日本一区视频| 一本色道久久99精品综合 | 欧美一区二区成人6969| 欧美日韩性视频在线| 亚洲电影中文字幕| 久久精品二区三区| 正在播放日韩| 欧美日韩精品免费观看视频完整| 激情成人av| 久久精品人人做人人综合| 一区二区三区欧美成人| 欧美精品色综合| 亚洲精品一区二区三区樱花| 久久免费视频在线| 亚洲综合色在线| 国产精品美女在线观看| 99国产精品久久久久久久久久| 美日韩精品免费| 久久久另类综合| 狠狠色综合色综合网络| 欧美亚洲一区二区在线| 在线一区二区视频| 欧美视频观看一区| 中文国产成人精品久久一| 亚洲人成久久| 欧美理论片在线观看| 亚洲蜜桃精久久久久久久| 亚洲人成绝费网站色www| 欧美精品一区在线发布| 亚洲精品在线免费| 亚洲日本久久| 欧美日韩午夜在线|