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

            搜索

            最新評論

            閱讀排行榜

            評論排行榜

            中文字幕久久精品无码| 欧美午夜A∨大片久久| 久久久久av无码免费网| 2020国产成人久久精品| 亚洲精品NV久久久久久久久久| 色悠久久久久久久综合网| 欧美精品丝袜久久久中文字幕| 一日本道伊人久久综合影| 久久乐国产综合亚洲精品| 亚洲精品tv久久久久久久久久| 久久香综合精品久久伊人| 久久国产精品偷99| 国产精品久久久久免费a∨| 欧美黑人激情性久久| 久久高清一级毛片| 久久91精品国产91久久麻豆| 久久精品国产精品亚洲下载| 久久婷婷五月综合97色一本一本 | 成人综合久久精品色婷婷| 亚洲色欲久久久综合网| 品成人欧美大片久久国产欧美...| 一本久久a久久精品综合香蕉 | 性做久久久久久久久久久| 午夜人妻久久久久久久久| 国产视频久久| 久久精品99久久香蕉国产色戒 | 久久精品久久久久观看99水蜜桃| 日韩AV无码久久一区二区| 久久国产AVJUST麻豆| 色综合久久88色综合天天| 亚洲中文字幕久久精品无码APP | 国产精品成人无码久久久久久| 欧美精品乱码99久久蜜桃| 久久精品国产99久久久香蕉| 成人综合伊人五月婷久久| 亚洲国产欧洲综合997久久| 成人综合久久精品色婷婷| 国产精品热久久毛片| 久久青青草原国产精品免费| 久久久av波多野一区二区| 久久婷婷五月综合色高清|