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

               這個題一看就知道是求歐拉函數(shù)。歐拉函數(shù)描述的正式題意。歐拉函數(shù)的理解可以按照算法導(dǎo)論上面的說法,對0-N-1進行篩選素數(shù)。
            那么公式n∏(1-1/p),其中p是n的素數(shù)因子,就可以得到直觀的理解了。但是計算的時候,會將這個式子變形下,得到另外一個形式。
               如圖所示:

               但是這個題,需要考慮下,有可能n是個大素數(shù),直接進行因子分解的話會超時的。怎么辦了,只能在分解的時候判斷n是不是已經(jīng)成為
            素數(shù)了,如果是素數(shù),答案再乘以n-1就行了。為了加快判斷,我用5mb的空間搞了個素數(shù)表,大于5000000的數(shù)字只能循環(huán)判斷了。

             代碼如下,注意求歐拉函數(shù)的代碼部分:
            #include <stdio.h>
            #include <math.h>
            #define MAX (5000000)
            bool bPrime[MAX];//false表示素數(shù)

            void InitPrime()
            {
                bPrime[0] = bPrime[1] = true;
                int nMax = sqrt((double)MAX) + 1;
                for (int i = 2; i <= nMax; ++i)
                {
                    if (!bPrime[i])
                    for (int j = i * 2; j < MAX; j += i)
                    {
                        bPrime[j] = true;
                    }
                }
            }

            bool IsPrime(int nN)
            {
                if (nN < MAX)
                {
                    return !bPrime[nN];
                }
                else
                {
                    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), nN)
                {
                    if (nN == 1){printf("0\n");continue;}
                    int nAns = 1;
                    for (int i = 2; i <= nN; ++i)
                    {
                        if (IsPrime(nN))
                        {
                            nAns *= nN - 1;
                            break;
                        }
                        if (nN % i == 0)
                        {
                            nAns *= i - 1;
                            nN /= i;
                            while (nN % i == 0)
                            {
                                nAns *= i;
                                nN /= i;
                            }
                        }
                    }
                    printf("%d\n", nAns);
                }
                
                return 0;
            }

            posted on 2012-07-26 21:36 yx 閱讀(728) 評論(0)  編輯 收藏 引用 所屬分類: 數(shù)論

            <2012年8月>
            2930311234
            567891011
            12131415161718
            19202122232425
            2627282930311
            2345678

            導(dǎo)航

            統(tǒng)計

            公告

            常用鏈接

            留言簿(3)

            隨筆分類

            隨筆檔案

            me

            好友

            同學(xué)

            網(wǎng)友

            搜索

            最新評論

            閱讀排行榜

            評論排行榜

            中文字幕久久精品| 国产免费福利体检区久久| 国产精品日韩深夜福利久久| 亚洲AV日韩精品久久久久| 亚洲国产精品久久久天堂| 777久久精品一区二区三区无码| 麻豆AV一区二区三区久久| 色偷偷久久一区二区三区| 亚洲国产另类久久久精品小说| 亚洲级αV无码毛片久久精品| 要久久爱在线免费观看| 青青草国产97免久久费观看| AV色综合久久天堂AV色综合在| 国产美女久久精品香蕉69| 国产亚洲精品美女久久久| 久久久久免费看成人影片| 久久亚洲精品视频| 久久免费看黄a级毛片| 亚洲国产一成人久久精品| 久久99国产亚洲高清观看首页| 大香伊人久久精品一区二区| 欧美精品久久久久久久自慰| 亚洲午夜久久久久妓女影院| 69久久夜色精品国产69| 无码国内精品久久人妻麻豆按摩| 久久精品免费网站网| 99久久中文字幕| 久久综合久久综合亚洲| 久久久综合香蕉尹人综合网| 久久国产精品久久国产精品| 成人亚洲欧美久久久久| 99久久精品免费看国产一区二区三区 | 色综合久久天天综线观看| 久久久久久国产精品免费免费| 九九99精品久久久久久| 色综合久久久久综合99| 日产精品久久久一区二区| 精品久久久久久无码中文字幕| 久久97精品久久久久久久不卡| 新狼窝色AV性久久久久久| 日韩久久无码免费毛片软件|