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

            bon

              C++博客 :: 首頁 :: 聯系 :: 聚合  :: 管理
              46 Posts :: 0 Stories :: 12 Comments :: 0 Trackbacks

            常用鏈接

            留言簿(2)

            我參與的團隊

            搜索

            •  

            最新評論

            閱讀排行榜

            評論排行榜

            PKU 2992
            給出n跟k,求有多少個約數。
            思路如下:
            要知道一個數n約數的個數,首先要知道n的質因子分解:

                                             

            接著就可以按下式來算所有約數的個數了:

                                             

            其次,要知道一個數n的質因子分解,要打素數表,然后對每個小于n的素數p試除,看n含有多少個p的因子。

            題目是要求n!/(k!(n-k)!的約數個數,如果直接算出n!/(k!(n-k)!會很麻煩,可以分別對n!, k!, (n-k)!做如上分解,然后對于同一個底,用n在的指數減去k對應的指數及(n-k)對應的指數,就得到n!/(k!(n-k)!在底的指數了。
            可以用遞推的方法來求n!的質因子分解,具體見代碼。
             1 // pku 2992 質因子分解
             2 #include <iostream>
             3 #include <math.h>
             4 using namespace std;
             5 
             6 int p[90];        // 素數表
             7 int pn;
             8 int e[432][90];    // 對應素數的指數,質因子分解
             9 int n,k;
            10 
            11 // 篩法復雜度為O(n^2),下面的試除法只要sigma(sqrt(i)),2<=i<=n
            12 void prime()
            13 {
            14     int i,j;
            15     pn=0;
            16     p[pn++]=2;
            17     for(i=3;i<=431;i++){
            18         int up=(int)sqrt(i);
            19         int flag=0;
            20         for(j=2;j<=up && !flag;j++){
            21             if(i%j==0) flag=1;
            22         }
            23         if(flag==0) p[pn++]=i;
            24     }
            25     //printf("cnt=%d\n",pn);
            26     //for(i=0;i<pn;i++) printf("%d ",p[i]);
            27     return;
            28 }
            29 
            30 void fact()
            31 {
            32     int i,j;
            33     for(i=2;i<=431;i++){
            34         int a=i;
            35         //printf("%d:",a);
            36         for(j=0;j<pn;j++){
            37             while(a>1 && a%p[j]==0) {e[i][j]++;a/=p[j];}
            38             //printf("%d ",e[i][j]);
            39         }
            40         //printf("\n");
            41     }
            42     for(i=3;i<=431;i++){
            43         for(j=0;j<pn;j++) e[i][j]+=e[i-1][j];
            44     }
            45 }
            46 
            47 /*
            48 void divide(int a)
            49 {
            50     int i;
            51     for(i=0;i<pn;i++){
            52         while(a%p[i]==0) {e[e]--;a/=p[i];}
            53     }
            54 }
            55 */
            56 
            57 int main()
            58 {
            59     //freopen("in.txt","r",stdin);
            60     int i,j;
            61     prime();
            62     fact();
            63     while(scanf("%d%d",&n,&k)!=EOF){
            64     //for(n=0;n<=431;n++){
            65         //for(k=0;k<n;k++){
            66             __int64 sum=1;
            67             for(i=0;i<pn;i++) sum*=(1+(e[n][i]-e[k][i]-e[n-k][i]));
            68             printf("%I64d\n",sum);
            69         //}
            70         //printf("1\n");
            71     }
            72     return 1;
            73 }
            posted on 2008-03-14 13:57 bon 閱讀(297) 評論(0)  編輯 收藏 引用
            Google PageRank 
Checker - Page Rank Calculator
            伊人久久精品无码av一区| 狠狠人妻久久久久久综合蜜桃| 国内精品久久久久影院亚洲| 午夜精品久久久久久毛片| 久久狠狠高潮亚洲精品| 亚洲精品视频久久久| 国产精品无码久久久久久| 久久久久久av无码免费看大片 | 亚洲国产精品一区二区久久| 久久久亚洲精品蜜桃臀 | 久久99久久成人免费播放| 亚洲午夜精品久久久久久浪潮| 青青草原综合久久| 久久久久久亚洲Av无码精品专口 | 久久午夜伦鲁片免费无码| 国产日韩久久久精品影院首页 | 一本色道久久88综合日韩精品| 91精品观看91久久久久久| 亚洲精品乱码久久久久久中文字幕| 欧美精品福利视频一区二区三区久久久精品 | 久久亚洲精品无码VA大香大香| 久久96国产精品久久久| 久久精品国产亚洲AV蜜臀色欲 | 久久国产精品一区| 国产精品视频久久久| 亚洲国产另类久久久精品小说| 一本色道久久综合| 合区精品久久久中文字幕一区| 久久人人爽人人爽人人片AV麻豆| 久久亚洲国产精品一区二区| 97久久天天综合色天天综合色hd| 亚洲午夜久久久久久久久电影网| 久久国产亚洲精品| 亚洲乱码日产精品a级毛片久久| 久久久久国产| 亚洲欧美一级久久精品| 亚洲国产成人久久精品99| 香蕉久久永久视频| 伊人久久五月天| 亚洲第一极品精品无码久久| 亚洲AV成人无码久久精品老人 |