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

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>
            久久久久久久999| 久久精品一二三| 国产原创一区二区| 久久成人精品视频| 欧美激情一区二区三区不卡| 亚洲狼人精品一区二区三区| 欧美午夜在线观看| 免费在线观看精品| 亚洲四色影视在线观看| 欧美高清视频一二三区| 亚洲综合色视频| 亚洲视频一区在线观看| 国产日韩精品在线| 欧美日韩中文精品| 欧美精品一线| 欧美国产日韩一区二区三区| 久久亚洲精品一区二区| 亚洲欧美在线磁力| 亚洲午夜久久久| 99在线精品观看| 亚洲精品视频在线播放| 亚洲国产精品尤物yw在线观看| 亚洲欧美区自拍先锋| 在线不卡中文字幕| 亚洲国产午夜| 国产一区二区久久精品| 国产亚洲欧洲997久久综合| 国产亚洲一区二区精品| 影音国产精品| 在线免费观看成人网| 亚洲电影在线观看| 亚洲最黄网站| 久久久久久久综合日本| 久久另类ts人妖一区二区| 欧美一区二区三区免费观看视频 | 久久视频免费观看| 亚洲一区二区精品在线| 性8sex亚洲区入口| 欧美一区二区三区日韩视频| 欧美黄色aaaa| 久久看片网站| 欧美屁股在线| 欧美成人一区二区三区在线观看| 亚洲性人人天天夜夜摸| 亚洲作爱视频| 亚洲卡通欧美制服中文| 国产精品视频一二| 久久在线观看视频| 久久久久在线| 免费观看久久久4p| 免费日韩av片| 欧美成人免费视频| 欧美午夜电影网| 欧美色大人视频| 国产精品久久国产三级国电话系列| 国产精品乱子久久久久| 欧美日本网站| 国产婷婷色一区二区三区在线 | 国产精品一区二区你懂的| 欧美精品www| 国产欧美视频一区二区| 亚洲美女区一区| 亚洲少妇自拍| 可以看av的网站久久看| 亚洲精品视频一区| 久久久最新网址| 国产日韩精品久久| 亚洲国产中文字幕在线观看| 亚洲精品免费在线播放| 亚洲影音一区| 国产精品日日摸夜夜添夜夜av| 国产日韩欧美中文在线播放| 日韩一级片网址| 欧美成人四级电影| 亚洲免费在线| 国产欧美日韩一区二区三区在线观看| 亚洲精品美女在线观看播放| 亚洲精品国产精品久久清纯直播 | 最近看过的日韩成人| 亚洲无线一线二线三线区别av| 免费国产自线拍一欧美视频| 狠狠干综合网| 亚洲性夜色噜噜噜7777| 久久夜色精品| 久久福利影视| 国产亚洲精久久久久久| 久久久91精品国产一区二区三区 | 国产日韩1区| 香蕉久久国产| 亚洲欧美电影在线观看| 国产精品大片| 亚洲欧美久久久| 中文网丁香综合网| 国产精品―色哟哟| 欧美国产综合| 国产精品久久999| 老司机午夜免费精品视频 | 久久国产免费| 欧美在线www| 一本大道久久精品懂色aⅴ| 中文亚洲字幕| 一区二区视频免费在线观看 | 99riav1国产精品视频| 国产精品久久久久久久7电影| 亚洲欧美制服中文字幕| 老司机一区二区| 香蕉乱码成人久久天堂爱免费| 欧美日韩综合网| 麻豆freexxxx性91精品| 欧美成人精品h版在线观看| 一区二区不卡在线视频 午夜欧美不卡'| 美国成人毛片| 欧美吻胸吃奶大尺度电影| 久久裸体视频| 国产精品嫩草99av在线| 欧美sm极限捆绑bd| 国产精品进线69影院| av成人黄色| 国产精品久久久久久久久久妞妞| 91久久国产综合久久| 免费看黄裸体一级大秀欧美| 亚洲系列中文字幕| 亚洲第一二三四五区| 亚洲国产精品久久| 亚洲一区二区三区中文字幕在线| 亚洲自拍偷拍一区| 久久久中精品2020中文| 99视频有精品| 亚洲午夜羞羞片| 欧美一级专区| 亚洲国产日韩美| 亚洲美女视频在线观看| 亚洲国产日韩综合一区| 这里只有精品电影| 久久青草久久| 欧美理论视频| 欧美色欧美亚洲另类七区| 亚洲尤物在线视频观看| 久久国产免费看| 欧美精品1区| 狠狠色综合网| 久久亚洲精品伦理| 久久国产精品免费一区| 国产乱码精品一区二区三区五月婷| 韩日视频一区| 亚洲欧美日韩在线播放| 欧美黑人国产人伦爽爽爽| 99国产精品国产精品久久| 欧美人牲a欧美精品| 亚洲精品视频在线看| 一区二区电影免费观看| aⅴ色国产欧美| 午夜精品久久久久久久| 亚洲欧洲在线一区| 每日更新成人在线视频| 精品福利免费观看| 亚洲在线观看| av不卡在线| 免费成人激情视频| 亚洲电影一级黄| 久久久精品国产一区二区三区| 亚洲精品网址在线观看| 久久在线播放| 狠狠色丁香久久婷婷综合丁香| 性色av一区二区三区红粉影视| 欧美成年人网| 欧美在线视频不卡| 国产日韩综合| 久久亚洲一区| 亚洲欧美日韩精品久久奇米色影视| 蜜桃伊人久久| 亚洲精品欧洲精品| 99精品欧美一区二区蜜桃免费| 欧美日韩不卡视频| 亚洲一区区二区| 亚洲精品乱码久久久久久久久| 欧美高潮视频| 亚洲电影免费观看高清完整版在线观看| 久久综合色8888| 久久亚洲精品一区| 日韩亚洲欧美一区| 亚洲免费伊人电影在线观看av| 国产性色一区二区| 亚洲第一精品在线| 国产精品免费久久久久久| 久久久久青草大香线综合精品| 欧美日韩免费一区二区三区视频| 久久夜色精品| 国产精品亚洲а∨天堂免在线| 免费在线成人| 国产精品成人免费| 99在线精品视频| 欧美一区2区视频在线观看 | 亚洲乱码日产精品bd| 欧美日韩精品二区第二页| 午夜精品一区二区三区在线| 欧美日韩免费在线视频| 蜜桃精品久久久久久久免费影院| 国产精品爽爽爽| 欧美激情一区二区|