• <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>

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

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

            導航

            統計

            公告

            常用鏈接

            留言簿(3)

            隨筆分類

            隨筆檔案

            me

            好友

            同學

            網友

            搜索

            最新評論

            閱讀排行榜

            評論排行榜

            国内精品久久久久久久影视麻豆| 91久久精品国产成人久久| 人妻无码αv中文字幕久久琪琪布| 久久人人超碰精品CAOPOREN| 亚洲欧美日韩久久精品| 99精品国产在热久久| 无码人妻久久一区二区三区蜜桃 | 久久久久久亚洲精品不卡| 亚洲日韩中文无码久久| 久久国产香蕉视频| 国产午夜精品久久久久免费视| 欧洲国产伦久久久久久久 | 国产精品久久久天天影视| 三级片免费观看久久| 青青草国产成人久久91网| 无码日韩人妻精品久久蜜桃 | 精品久久久无码人妻中文字幕豆芽| 久久久久亚洲AV无码专区桃色| 久久久久四虎国产精品| 久久人人妻人人爽人人爽| 伊人久久大香线蕉精品不卡| 亚洲国产成人久久综合碰碰动漫3d| 久久人人爽人人爽人人AV东京热 | 亚洲精品乱码久久久久久蜜桃 | 天天爽天天狠久久久综合麻豆| 亚洲午夜久久久| 久久久久久久综合狠狠综合| 亚洲一区精品伊人久久伊人| 久久久久99这里有精品10| 午夜精品久久久久久久无码| 亚洲人成网站999久久久综合| 波多野结衣久久一区二区 | 国产精品热久久毛片| 狠狠人妻久久久久久综合| 久久伊人中文无码| 日本精品久久久久久久久免费| 亚洲色欲久久久久综合网| 伊人色综合久久天天人手人婷| 欧美午夜精品久久久久免费视| 国内精品久久久久久99蜜桃| 久久久精品一区二区三区|