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

            misschuer

            常用鏈接

            統計

            積分與排名

            百事通

            最新評論

            nuaa 1284 1/N! = 1/X + 1/Y

            #include <iostream>
            #define MAX 10001
            using namespace std;
            bool prime[ MAX ] = {false};
            int r[ MAX ];
            int ff[ MAX ] , kf;
            bool v[ MAX ];
            int rank[ MAX ];
            __int64 ans[MAX];
            int k;

            void factor (int n)
            {
                
            int i;
                
            for (i = 2 ; i <= n ; ++ i)
                
            {
                    
            while (n % i == 0)
                    
            {
                        n 
            /= i;
                        r[ i ] 
            ++;
                    }

                    
            if (n == 1)
                        
            break;
                    
            if (prime[ n ])
                    
            {
                        r[ n ] 
            ++;
                        
            break;
                    }

                }

            }


            int main()
            {
                
            int i , j;
                
            int n , g;
                __int64   tt;
                kf 
            = 0;
                
            for (i = 2 ; i <= MAX ; ++ i)
                    prime[ i ] 
            = true;
                
            for (i = 2 ; i <= MAX ; ++ i)
                    
            if (prime[ i ])
                    
            {
                        ff[
            ++ kf] = i;
                        rank[ i ] 
            = kf;
                        
            for (j = i + i ; j <= MAX ; j += i)
                            prime[ j ] 
            = false;
                    }

                    
                    
            while (scanf ("%d" , &n) , n)
                    
            {
                        g 
            = 2;
                        memset (r , 
            0 , sizeof (r));
                        
            for (i = 2 ; i <= n ; ++ i)
                        
            {
                            
            if (prime[ i ])
                            
            {
                                r[ i ] 
            ++;
                                
            if (i > g)
                                    g 
            = i;
                            }

                            
            else
                                factor (i);
                        }

                        k 
            = 1;
                        ans[ 
            1 ] = 1;
                        
            int tmp = 0;
                        
            for (i = 1 ; i <= rank[ g ] ; ++ i)
                        
            {
                            tt 
            = (r[ff[ i ]] << 1+ 1;
                            
            for (j = 1 ; j <= k ; ++ j)
                            
            {
                                ans[ j ] 
            = ans[ j ] * tt + tmp;
                                tmp 
            = ans[ j ] / 10000 , ans[ j ] %= 10000;
                            }

                            
            while (tmp)
                                ans[
            ++ k] = tmp % 10000 , tmp /= 10000;
                        }

                        printf (
            "%I64d" , ans[ k ]);
                        
            for (i = k - 1 ; i > 0 ; -- i)
                            printf (
            "%04I64d" , ans[ i ]);
                        printf (
            "\n");
                    }

                    
            return 0;
            }
            這題其實是求(n!)^2的約數的個數,注意要大數處理。其余的自己去算吧,orz~~~~~~

            posted on 2009-12-23 12:51 此最相思 閱讀(123) 評論(0)  編輯 收藏 引用

            久久天天躁狠狠躁夜夜av浪潮| 国产精品99久久精品| 青青国产成人久久91网| 亚洲精品乱码久久久久久蜜桃 | 精品无码久久久久久午夜| 久久精品一区二区影院| 波多野结衣AV无码久久一区| 精品久久久久久国产三级| 性欧美大战久久久久久久久| 久久精品国产99久久久香蕉| 18禁黄久久久AAA片| 久久这里有精品| 精品久久久久久久久久中文字幕| 婷婷久久综合九色综合绿巨人| 国产成人精品久久一区二区三区av| 久久久精品人妻一区二区三区蜜桃 | 伊人伊成久久人综合网777| 精品久久久久久亚洲精品| 四虎影视久久久免费观看| 久久夜色精品国产亚洲| 久久精品国产99久久久古代| 色婷婷综合久久久久中文字幕 | 久久国产精品波多野结衣AV| 蜜臀av性久久久久蜜臀aⅴ麻豆| 精品伊人久久大线蕉色首页| 精品久久久久久99人妻| 久久66热人妻偷产精品9| 亚洲色婷婷综合久久| 性高湖久久久久久久久AAAAA| www亚洲欲色成人久久精品| 2022年国产精品久久久久| 国产亚洲精久久久久久无码77777| 久久人搡人人玩人妻精品首页 | 免费精品国产日韩热久久| 久久久久99精品成人片| 久久精品国产一区二区三区不卡| 久久99国产精品久久99| 欧美无乱码久久久免费午夜一区二区三区中文字幕 | 狠狠久久综合| 精品熟女少妇aⅴ免费久久| 91精品婷婷国产综合久久|