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

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>
            亚洲免费福利视频| 亚洲免费中文| 欧美一二三区在线观看| 亚洲永久免费观看| 亚洲一区二区精品| 欧美一级视频| 久久精品中文字幕一区二区三区| 香蕉久久精品日日躁夜夜躁| 久久精品欧洲| 亚洲大片免费看| 亚洲国产日韩欧美在线动漫| 99re6这里只有精品视频在线观看| 日韩午夜激情| 久久av资源网| 欧美日本一区二区三区| 国产一本一道久久香蕉| 亚洲人屁股眼子交8| 亚洲欧美另类综合偷拍| 免费观看在线综合| 一区二区三区精密机械公司| 久久岛国电影| 欧美视频一二三区| 一区二区三区无毛| 亚洲男女自偷自拍图片另类| 免费欧美视频| 亚洲欧美中文另类| 欧美日韩国产综合网| 黑人一区二区| 性刺激综合网| 亚洲美女av电影| 美女精品一区| 国产真实乱子伦精品视频| 亚洲性av在线| 亚洲日韩欧美视频一区| 久久精品最新地址| 国产精品免费小视频| 日韩亚洲综合在线| 欧美成人精品| 欧美制服丝袜第一页| 欧美三级视频在线| 亚洲精品一区二区三区樱花| 久久理论片午夜琪琪电影网| 亚洲午夜免费视频| 欧美午夜电影网| 一本色道久久综合| 亚洲国产日韩一区| 欧美高清成人| 日韩视频精品在线| 欧美风情在线观看| 久久香蕉国产线看观看av| 国产日产欧美一区| 欧美二区视频| 欧美高清免费| 亚洲激情亚洲| 欧美成人精品不卡视频在线观看 | 亚洲精品五月天| 久久精品一本| 亚洲在线视频免费观看| 欧美日韩另类丝袜其他| 日韩一级网站| 亚洲美女av电影| 欧美午夜精品久久久久久人妖| 99国产精品视频免费观看一公开| 欧美激情欧美激情在线五月| 麻豆国产精品一区二区三区| 亚洲国内精品| 亚洲精品久久久久久下一站| 欧美日韩日日夜夜| 亚洲综合不卡| 欧美一区日韩一区| 在线观看欧美视频| 亚洲黄色高清| 国产精品va在线| 久久精品国产999大香线蕉| 久久国产精品第一页 | 欧美电影免费观看大全| 亚洲欧美一区二区原创| 国产精品视频自拍| 久久精品一区二区国产| 久久精品91久久久久久再现| 亚洲电影免费| 亚洲乱码国产乱码精品精天堂| 欧美日韩麻豆| 久久精品视频免费| 欧美成人午夜激情| 亚洲中午字幕| 久久久精品国产免费观看同学| 亚洲国产三级在线| 一本色道久久综合亚洲91| 在线高清一区| 麻豆精品一区二区综合av| 另类图片国产| 亚洲精品极品| 午夜精品福利视频| 亚洲精品国精品久久99热| 一区二区高清视频| 在线精品一区二区| 亚洲视频网站在线观看| 在线免费观看欧美| 一本久道久久综合中文字幕| 国产一区二区看久久| 亚洲欧洲日本mm| 国产丝袜一区二区三区| 亚洲激情第一页| 国产在线播放一区二区三区| 亚洲国产欧美国产综合一区| 国内偷自视频区视频综合| 欧美黄色大片网站| 国产视频自拍一区| 国产精品乱码| 久久久青草青青国产亚洲免观| 在线观看成人av| 这里只有精品丝袜| 亚洲国产天堂久久综合| 亚洲一区久久| 亚洲另类在线一区| 久久成人免费| 亚洲综合欧美日韩| 欧美刺激午夜性久久久久久久| 欧美中文日韩| 欧美午夜一区二区| 亚洲激情电影中文字幕| 永久免费视频成人| 亚洲欧美激情视频| 亚洲午夜av| 欧美国产一区在线| 欧美a级一区二区| 国产一区二区久久精品| 亚洲图片欧洲图片av| 亚洲国产精品一区二区久| 亚洲电影av在线| 亚洲尤物在线| 中日韩美女免费视频网站在线观看| 午夜在线观看免费一区| 亚洲一区影音先锋| 欧美激情一区二区三区成人| 久久久久久久久久看片| 国产精品午夜国产小视频| 夜夜嗨av一区二区三区网站四季av| 亚洲精品国产视频| 欧美高清一区二区| 亚洲国产高清一区| 亚洲日本精品国产第一区| 美女精品一区| 亚洲黄色在线| 亚洲视频观看| 国产精品捆绑调教| 亚洲欧美国产高清va在线播| 欧美一区二区高清| 国产欧美日韩不卡免费| 午夜精品美女久久久久av福利| 午夜视黄欧洲亚洲| 国产欧美一区二区白浆黑人| 性亚洲最疯狂xxxx高清| 麻豆亚洲精品| 亚洲精品少妇| 国产精品久久久久久久久久妞妞| 亚洲欧美日韩另类| 蜜臀av国产精品久久久久| 亚洲人成亚洲人成在线观看| 欧美日韩另类国产亚洲欧美一级| 亚洲午夜精品久久| 久久久久国产免费免费| 亚洲国产福利在线| 欧美性片在线观看| 久久本道综合色狠狠五月| 欧美激情一区二区三区全黄| 美女国产一区| 一区二区国产在线观看| 欧美高清在线播放| 亚洲一区一卡| 亚洲国产高清自拍| 午夜在线精品偷拍| 亚洲国产精品久久久久久女王| 欧美精品国产精品| 午夜久久影院| 亚洲激情在线视频| 久久精品国产精品亚洲精品| 亚洲免费观看高清在线观看| 国产日韩精品在线| 欧美日韩高清免费| 久久在线精品| 午夜一级久久| 夜夜嗨av一区二区三区网页| 欧美成在线观看| 久久九九免费视频| 亚洲图片欧美日产| 亚洲黄色免费电影| 国产一区二区三区四区五区美女 | 欧美色图一区二区三区| 久久久精品国产一区二区三区| 亚洲人成在线播放网站岛国| 久久久久久久久久码影片| 亚洲一级黄色av| 亚洲免费精品| 91久久黄色| 在线精品国产成人综合| 国产精品入口福利| 国产精品成人免费| 欧美精品 日韩|