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

poj 1284 Primitive Roots

   這個題是求原根的個數。所謂原根,意思是給定一個數n,存在數g,g^j能夠產生乘法群Zn*中所有的數字。即g^j = {x|x與n互質,
1<=x<n}。如果n是奇素數p(大于2的素數),那么滿足g^j={1,2,...,p-1}。
   這個題目要求求原根的個數。由費馬定理由,對任意1<=x<p,即Zp*中的數字,都由x^(p-1) = 1 % p。從費馬定理可以看出,
再往下計算就開始循環了。那么有,x^i%p(1<=i<p) = {1, 2, 3,...,p-1},意思是能夠生成Zp*中的所有數字。
   根據上面的那個式子可以得到,x^i%(p-1)(1<=i<p) = {0, 1, 2,...,p-2}。 如果由gcd(x,p-1) = 1,那么必然存在某個x^i,
使得x^i*x = (p-1)%p。
   因此可以得到,原根的個數是p-1的乘法群中元素的個數,也就是歐拉函數(p-1)。

代碼如下:
   
#include <stdio.h>
#include <math.h>
#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;
            }
        }
    }
}
bool IsPrime(int nN)
{
    if (nN < MAX)return !bPrime[nN];
    int nMax = sqrt((double)nN) + 1;
    for (int i = 2; i <= nMax; ++i)
    {
        if (nN % i == 0)
            return false;
    }
    return true;
}
int main()
{
    int nN;
    InitPrime();
    while (scanf("%d", &nN) == 1)
    {
        nN--;
        int nAns = 1;
        if (IsPrime(nN))
        {
            nAns = nN - 1;
        }
        else
        {
            for (int i = 2; i <= nN; ++i)
            {
                if (nN % i == 0)
                {
                    nAns *= i - 1;
                    nN /= i;
                    while (nN % i == 0)
                    {
                        nAns *= i;
                        nN /= i;
                    }
                    if (IsPrime(nN))
                    {
                        nAns *= nN - 1;
                        break;
                    }
                }
            }
        }
        printf("%d\n", nAns);
    }
    return 0;
}

posted on 2012-07-30 22:36 yx 閱讀(930) 評論(0)  編輯 收藏 引用 所屬分類: 數論

<2012年7月>
24252627282930
1234567
891011121314
15161718192021
22232425262728
2930311234

導航

統計

公告

常用鏈接

留言簿(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>
            中文网丁香综合网| 国产精品美女久久久久久久| 国产精品久久久久久久久久久久 | 另类欧美日韩国产在线| 一区二区精品在线| 欧美精品成人91久久久久久久| 黄色日韩在线| 性色av香蕉一区二区| 久久在线视频在线| 久久激情婷婷| 韩日欧美一区二区| 久久另类ts人妖一区二区| 欧美一区二区三区在线| 中文国产成人精品久久一| 欧美日韩喷水| 亚洲在线播放| 亚洲自啪免费| 国产香蕉97碰碰久久人人| 欧美视频三区在线播放| 欧美大片免费观看在线观看网站推荐| 一区在线观看视频| 免费久久99精品国产自| 麻豆免费精品视频| 亚洲伦伦在线| 99国产精品视频免费观看| 欧美亚洲第一区| 午夜综合激情| 久久都是精品| 亚洲人成欧美中文字幕| 亚洲精品综合精品自拍| 欧美性色视频在线| 久久精品中文| 欧美不卡福利| 亚洲欧美日韩一区二区在线| 亚洲欧美另类在线观看| 午夜精品久久一牛影视| 国产性做久久久久久| 男女激情久久| 欧美日韩国产综合久久| 亚洲欧美另类在线观看| 久久精品1区| 亚洲免费精彩视频| 亚洲男同1069视频| 亚洲国产综合在线看不卡| 亚洲日韩欧美视频一区| 精品成人一区二区三区| 亚洲女优在线| 久久国产精品亚洲va麻豆| 91久久国产综合久久| 亚洲色图制服丝袜| 1000部国产精品成人观看| 99精品久久久| 18成人免费观看视频| 一区二区三区鲁丝不卡| 在线观看视频亚洲| 亚洲一区二区视频在线| 亚洲欧洲在线观看| 亚洲欧美日韩高清| 亚洲美女免费视频| 久久岛国电影| 午夜精品久久一牛影视| 欧美电影免费网站| 亚洲欧美国产高清| 久久久7777| 亚洲欧美激情一区| 久久乐国产精品| 香蕉免费一区二区三区在线观看| 久久亚洲精品中文字幕冲田杏梨| 午夜精品电影| 91久久综合| 久久精品国产99精品国产亚洲性色| 99这里只有久久精品视频| 久久亚洲国产精品一区二区| 欧美一级片在线播放| 欧美日韩在线影院| 最新国产の精品合集bt伙计| 国产乱码精品一区二区三区av| 最新亚洲视频| 亚洲日本免费| 久久夜色精品| 美女亚洲精品| 精品盗摄一区二区三区| 午夜免费日韩视频| 亚洲欧美日本日韩| 亚洲资源av| 亚洲一区国产| 欧美视频一区在线| 日韩视频在线永久播放| 日韩视频一区二区在线观看| 男人的天堂成人在线| 欧美成人精品h版在线观看| 欧美三区免费完整视频在线观看| 国产一区二区精品在线观看| 亚洲男人的天堂在线| 午夜伦理片一区| 国产麻豆日韩欧美久久| 午夜欧美大尺度福利影院在线看| 亚洲欧美另类国产| 国产精品无码永久免费888| 亚洲尤物在线| 久久久精品动漫| 国产自产精品| 噜噜噜躁狠狠躁狠狠精品视频| 欧美a级一区二区| 亚洲人成啪啪网站| 欧美日韩国产色视频| 国产精品99久久99久久久二8| 午夜精品久久久99热福利| 久久国产欧美精品| 老司机午夜精品视频在线观看| 亚洲国产精品久久精品怡红院| 亚洲性视频网站| 欧美日韩一区二区在线观看| 一区二区免费在线播放| 午夜精品久久久久| 亚洲欧美国产日韩中文字幕| 国产精品影院在线观看| 欧美在线免费视屏| 蜜桃视频一区| 日韩视频在线免费观看| 欧美午夜精品久久久久久浪潮| 亚洲尤物视频网| 免费成人激情视频| 99这里只有精品| 国产精品视频一区二区高潮| 欧美在线www| 亚洲经典三级| 欧美在线视频网站| 亚洲人成在线影院| 国产精品一区二区久激情瑜伽| 久久免费视频在线观看| 亚洲欧美一区二区激情| 狠狠色综合一区二区| 欧美日韩国产精品成人| 欧美中文在线观看国产| 亚洲精品免费观看| 久久婷婷影院| 亚洲欧美清纯在线制服| 亚洲人人精品| 狠狠色综合网| 欧美影院在线播放| 国产综合精品| 欧美深夜影院| 欧美激情视频一区二区三区免费 | 国产欧美日韩在线| 欧美黑人多人双交| 久久国产一区二区三区| 一区二区日韩精品| 亚洲高清不卡在线观看| 久久久噜噜噜久久久| 亚洲女女女同性video| 亚洲欧洲三级| 欧美中日韩免费视频| 在线亚洲欧美视频| 亚洲国产精品视频一区| 久热精品视频| 久久久99爱| 欧美在线观看网站| 亚洲欧美日本精品| 亚洲欧美国产不卡| 久久久青草婷婷精品综合日韩| 免费成人你懂的| 久热精品视频在线免费观看| 性欧美videos另类喷潮| 亚洲制服丝袜在线| 一本一本久久a久久精品牛牛影视| 在线观看欧美日本| 狠狠色丁香久久婷婷综合丁香| 国产精品免费一区二区三区观看 | 亚洲高清在线观看一区| 国产亚洲日本欧美韩国| 国产精品无码专区在线观看| 欧美性猛交xxxx免费看久久久| 欧美日本亚洲韩国国产| 欧美极品aⅴ影院| 免费成人激情视频| 欧美成人国产| 免费成人av在线看| 亚洲精选成人| 亚洲老板91色精品久久| 亚洲三级免费| 亚洲美洲欧洲综合国产一区| 亚洲精品美女久久久久| 日韩亚洲欧美中文三级| 午夜精彩国产免费不卡不顿大片| 欧美亚州韩日在线看免费版国语版| 欧美欧美午夜aⅴ在线观看| 欧美激情精品久久久久久免费印度 | 亚洲毛片av| 一级日韩一区在线观看| 亚洲视频网在线直播| 亚洲欧美在线播放| 久久噜噜亚洲综合| 欧美精品首页| 国产精品久久影院| 一区在线视频观看| 日韩一区二区精品| 欧美一级视频| 欧美激情成人在线| 亚洲伦理在线观看|