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

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>
            欧美性猛交xxxx免费看久久久| 欧美影院在线播放| 乱人伦精品视频在线观看| 亚洲网站视频| 亚洲福利视频免费观看| 国产欧美三级| 国产精品日韩高清| 欧美在线视频观看免费网站| 中日韩美女免费视频网址在线观看| 亚洲成人在线视频网站| 欧美一区在线视频| 亚洲欧美日韩区| 午夜欧美不卡精品aaaaa| 亚洲精品国产精品乱码不99按摩| 国产亚洲精品久| 欧美日韩一级大片网址| 欧美日韩成人免费| 欧美视频观看一区| 国产精品福利片| 国产欧美日本在线| 国产精品视屏| 韩国成人理伦片免费播放| 欧美猛交免费看| 欧美精品在线观看一区二区| 欧美精品一区二区三区很污很色的 | 国产欧美亚洲日本| 欧美午夜不卡在线观看免费| 欧美日韩视频在线第一区| 久久人人爽人人爽| 欧美国产欧美亚州国产日韩mv天天看完整 | 一区精品在线| 亚洲欧洲视频| 亚洲一区精品电影| 一区二区三欧美| 久久精品久久综合| 国产麻豆午夜三级精品| 伊人精品成人久久综合软件| 一本色道久久综合| 久久久国产午夜精品| 日韩一区二区免费看| 久久久久成人精品免费播放动漫| 欧美丝袜一区二区| 最新成人在线| 久久亚洲图片| 亚洲欧美日韩中文播放| 久久综合精品国产一区二区三区| 国产精品国产三级欧美二区| 亚洲国产精品久久久久婷婷老年| 欧美自拍偷拍午夜视频| 亚洲日本精品国产第一区| 欧美在线观看日本一区| 国产精品扒开腿做爽爽爽视频| 亚洲精品日产精品乱码不卡| 老妇喷水一区二区三区| 午夜久久久久| 国产精品女人毛片| 亚洲午夜高清视频| 亚洲国产视频一区| 免费欧美电影| 亚洲精品一区二区三区樱花| 久久久久久久综合狠狠综合| 毛片一区二区三区| 精品成人a区在线观看| 性欧美1819sex性高清| 亚洲色在线视频| 欧美性猛交xxxx免费看久久久| 一区二区高清视频在线观看| 亚洲精品国产精品国产自| 欧美a级理论片| 日韩天堂在线视频| 99xxxx成人网| 国产精品xnxxcom| 亚洲在线中文字幕| 亚洲免费在线电影| 激情欧美一区| 欧美二区不卡| 欧美视频一区二区三区在线观看| 亚洲一区二区三区在线看 | 一区二区三区蜜桃网| 亚洲国产导航| 国产精品久久久久高潮| 这里只有精品电影| 亚洲一区二区三区777| 国产精品美女诱惑| 久久裸体视频| 欧美韩国一区| 亚洲欧美在线免费| 欧美在线视频免费| 亚洲区欧美区| 一二三区精品福利视频| 国产欧美视频一区二区三区| 欧美成人日本| 国产精品久久久久毛片大屁完整版 | 亚洲综合第一| 久久激情视频| 亚洲一区二区三区在线| 久久国产精品久久久久久电车| 亚洲黄页视频免费观看| av成人国产| 在线看片第一页欧美| 亚洲毛片在线免费观看| 激情欧美日韩一区| 一区二区高清视频在线观看| 国产日韩欧美亚洲| 99riav1国产精品视频| 在线观看视频一区二区| 99在线精品视频在线观看| 国内偷自视频区视频综合| 亚洲精品免费看| 国模叶桐国产精品一区| 一区二区三区视频在线| 亚洲欧洲精品一区| 性高湖久久久久久久久| 一区二区三区久久网| 久久综合成人精品亚洲另类欧美| 亚洲综合视频1区| 亚洲一区二区欧美日韩| 日韩一级精品| 国产一区二区三区免费观看| 最新国产の精品合集bt伙计| 国产欧美精品一区二区色综合| 亚洲精品美女| 亚洲黄色免费| 久久激情综合网| 午夜免费电影一区在线观看| 欧美久久久久免费| 亚洲第一色在线| 亚洲成色精品| 久久伊人一区二区| 久久中文字幕导航| 国产欧美欧美| 午夜激情一区| 久久久999成人| 国内精品久久久久影院色| 欧美一区观看| 久久久久久一区二区三区| 国产一区 二区 三区一级| 亚洲一区二区在线免费观看视频 | 国产精品99久久久久久宅男| 欧美成人精品在线| 亚洲成色777777在线观看影院| 精品成人一区二区| 久久精品一区二区三区不卡| 久久国产精品一区二区| 国产日韩精品一区二区三区在线| 亚洲在线国产日韩欧美| 久久久精品国产一区二区三区| 国产一区二区三区日韩欧美| 久久av一区二区| 欧美大片一区二区三区| 亚洲人成亚洲人成在线观看图片| 欧美激情久久久久| 日韩午夜av在线| 欧美一区二区三区四区高清 | 午夜精品网站| 国产精品一区免费观看| 久久精品视频一| 亚洲人成77777在线观看网| 日韩系列在线| 国产精品激情偷乱一区二区∴| 先锋影音国产一区| 另类激情亚洲| 99精品欧美一区二区三区综合在线| 欧美日韩中文精品| 欧美一区二区视频97| 欧美激情精品久久久| 亚洲一区二区网站| 国语自产精品视频在线看抢先版结局| 久久久久久尹人网香蕉| 亚洲精品一品区二品区三品区| 午夜精品短视频| 在线日韩av| 国产精品女人毛片| 美女诱惑一区| 亚洲一区二区精品在线观看| 女主播福利一区| 性做久久久久久免费观看欧美| 亚洲成色www久久网站| 99精品国产在热久久婷婷| 欧美一级欧美一级在线播放| 国产人成一区二区三区影院| 蜜臀91精品一区二区三区| 99在线精品观看| 欧美α欧美αv大片| 亚洲亚洲精品在线观看| 精品成人在线| 国产精品久久久久久久第一福利| 亚洲毛片在线免费观看| 国产精品免费网站| 欧美成人一品| 久久精品一区二区三区四区| 亚洲网站在线播放| 亚洲国产日韩欧美一区二区三区| 久久av一区二区三区| 一区二区三区回区在观看免费视频| 国产一区二区三区久久精品| 国产精品成人一区二区三区吃奶| 六月婷婷一区| 久久大逼视频| 久久国产精品久久w女人spa|