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

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

<2012年10月>
30123456
78910111213
14151617181920
21222324252627
28293031123
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>
            午夜精品理论片| 久久久久这里只有精品| 亚洲免费高清视频| 欧美另类久久久品| 亚洲性图久久| 午夜精品福利一区二区蜜股av| 国产精品久久久久久一区二区三区| 亚洲性图久久| 亚洲欧洲av一区二区三区久久| 国产日韩一区二区三区| 久久三级福利| 欧美好骚综合网| 亚洲一区二区三区在线观看视频| 亚洲一级二级在线| 红桃视频成人| 日韩午夜电影| 国产午夜精品在线观看| 亚洲成色777777女色窝| 欧美紧缚bdsm在线视频| 性色av一区二区三区在线观看| 欧美在线视频在线播放完整版免费观看| 黄色一区二区在线| 亚洲麻豆一区| 影音先锋久久精品| 一区二区国产日产| 亚洲春色另类小说| 亚洲视频图片小说| 亚洲国产精品久久久久秋霞影院 | 亚洲一区在线视频| 欧美在线日韩| 亚洲午夜av在线| 久久午夜电影网| 亚洲一区日韩在线| 美女任你摸久久| 欧美一区二区三区免费视频| 麻豆国产va免费精品高清在线| 亚洲一区二区在线看| 久久天天躁狠狠躁夜夜爽蜜月| 一区二区日韩欧美| 蜜臀va亚洲va欧美va天堂| 午夜在线观看欧美| 欧美精品日韩一区| 另类亚洲自拍| 国产精品一区二区三区乱码| 亚洲激情影视| 亚洲福利国产精品| 欧美伊人久久| 午夜精品区一区二区三| 欧美日韩 国产精品| 牛牛精品成人免费视频| 国产亚洲一区二区三区在线播放| 亚洲精品一区二区在线| 亚洲精品久久久久中文字幕欢迎你| 香蕉成人久久| 欧美一区二区三区喷汁尤物| 欧美丝袜一区二区三区| 亚洲精品一级| 一本一本大道香蕉久在线精品| 美女诱惑一区| 欧美福利视频| 亚洲高清中文字幕| 卡通动漫国产精品| 免费不卡视频| 亚洲高清网站| 欧美a级片网| 最新日韩在线视频| 亚洲精品视频中文字幕| 欧美国产日韩一区| 亚洲激情在线激情| 亚洲福利在线观看| 欧美国产三级| 亚洲每日更新| 午夜精品网站| 国产视频欧美| 久久久久久**毛片大全| 欧美成人四级电影| 日韩午夜精品| 国产精品扒开腿爽爽爽视频| 这里是久久伊人| 久久精品在这里| 经典三级久久| 欧美精品18videos性欧美| 一区二区三区高清不卡| 午夜精品久久久久久久久久久久久 | 性欧美长视频| 狠狠色狠狠色综合系列| 久久久久久久一区二区三区| 欧美承认网站| 亚洲在线一区| 国产字幕视频一区二区| 欧美成黄导航| 亚洲一区在线观看免费观看电影高清 | 欧美freesex8一10精品| 亚洲精选视频免费看| 午夜精品999| 在线观看日韩av| 欧美日韩视频一区二区| 欧美夜福利tv在线| 亚洲国产欧美一区二区三区丁香婷| 一区二区日韩免费看| 国产毛片一区二区| 免费黄网站欧美| 亚洲在线视频网站| 欧美激情视频一区二区三区不卡| 亚洲午夜女主播在线直播| 激情综合在线| 国产精品你懂的在线| 久久亚洲春色中文字幕| 亚洲调教视频在线观看| 欧美激情91| 久久精品91久久香蕉加勒比| 亚洲久久在线| 精品69视频一区二区三区| 欧美日韩一区自拍| 老司机午夜免费精品视频| 亚洲一区国产精品| 亚洲欧洲日本国产| 久久综合伊人| 久久精品女人天堂| 亚洲图片欧美日产| 最新精品在线| 激情文学综合丁香| 国产美女精品在线| 国产精品成人一区二区三区吃奶| 久久在线播放| 久久精品一区二区三区中文字幕 | 久久青草欧美一区二区三区| 亚洲一区二区三区免费视频| 亚洲人成在线影院| 在线日韩电影| 狠狠做深爱婷婷久久综合一区| 国产精品毛片一区二区三区| 欧美日韩精品伦理作品在线免费观看| 久久久久国产精品一区二区| 午夜精品一区二区在线观看| 99re热这里只有精品免费视频| 欧美大片国产精品| 欧美99在线视频观看| 久久综合九九| 蜜臀av一级做a爰片久久| 久久久久久久网站| 久久全国免费视频| 久久久另类综合| 久久综合久久88| 美女诱惑一区| 欧美成人精精品一区二区频| 欧美波霸影院| 亚洲国产精品热久久| 亚洲国产精品成人一区二区| 欧美激情1区| 亚洲精品久久久一区二区三区| 亚洲高清不卡在线| 亚洲黄色成人网| 9久草视频在线视频精品| 99视频+国产日韩欧美| 亚洲一二三区视频在线观看| 亚洲女人天堂av| 久久精品一区| 欧美成人免费视频| 欧美日韩一区三区| 国产女人精品视频| 黄色工厂这里只有精品| 亚洲国产精品成人久久综合一区| 在线精品高清中文字幕| 亚洲精品久久久一区二区三区| 亚洲看片网站| 午夜精品视频网站| 麻豆精品一区二区av白丝在线| 欧美chengren| 一区二区三区精品在线| 亚洲欧美日韩精品一区二区 | 宅男噜噜噜66一区二区66| 亚洲一级二级在线| 久久精品综合网| 欧美日本一道本在线视频| 国产精品久久久久aaaa九色| 国产一区二区三区久久| 亚洲欧洲午夜| 性感少妇一区| 亚洲国产精品黑人久久久| 一本到高清视频免费精品| 欧美在线视频一区二区| 欧美经典一区二区| 国产一区观看| 亚洲视频免费观看| 蜜臀久久久99精品久久久久久| 日韩小视频在线观看| 欧美专区亚洲专区| 欧美色另类天堂2015| 亚洲电影在线播放| 亚洲欧美日韩国产综合在线| 欧美~级网站不卡| 亚洲视频自拍偷拍| 欧美.www| 影音国产精品| 久久成人18免费网站| 亚洲欧洲视频在线| 久久深夜福利| 国产日韩在线不卡| 亚洲免费视频网站|