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

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 此最相思 閱讀(127) 評論(0)  編輯 收藏 引用

青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品
  • <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>
            欧美成人亚洲成人日韩成人| 亚洲综合第一页| 免费中文日韩| 一本久久精品一区二区| 欧美福利视频一区| 欧美大片免费观看| 欧美激情一级片一区二区| 久久综合一区二区| 国产一区二区av| 国产精品一区二区在线观看| 一区二区三区国产精华| 国产视频亚洲精品| 欧美一级在线亚洲天堂| 先锋亚洲精品| 日韩一级精品| 欧美在线视频一区| 久久成人综合网| 国产精品99免视看9| 久久国产精品99国产| 国产精品美女久久久免费| 国产农村妇女毛片精品久久麻豆| 亚洲无线视频| 久久米奇亚洲| 久久色在线观看| 欧美精品二区| 亚洲小少妇裸体bbw| 久久国产99| 久久av一区二区三区亚洲| 欧美精品日日鲁夜夜添| 亚洲精品国偷自产在线99热| 午夜精品在线视频| 亚洲综合国产精品| 欧美激情视频一区二区三区在线播放 | 欧美日韩午夜| 国产一区二区三区四区在线观看 | 欧美高清视频一区二区| 亚洲一区日本| 午夜精品区一区二区三| 欧美精品大片| 亚洲欧美日韩视频一区| 欧美激情网友自拍| 欧美三级电影一区| 亚洲欧洲免费视频| 亚洲美女黄网| 美女视频黄免费的久久| 亚洲免费观看高清完整版在线观看熊| 亚洲欧美日韩综合aⅴ视频| 欧美在线高清视频| 国产精品大片wwwwww| 久久久久久久久久久久久久一区| 亚洲欧洲精品一区| 尤物在线精品| 欧美国产免费| 香蕉久久夜色精品| 亚洲激情网址| 久久色中文字幕| 欧美韩国一区| 亚洲欧洲精品一区二区| 亚洲一区视频在线观看视频| 欧美久色视频| 久久久久网站| 欧美在线网址| 日韩亚洲欧美中文三级| 亚洲国产成人久久| 国产精品手机视频| 欧美在线观看天堂一区二区三区| 欧美插天视频在线播放| 亚洲高清激情| 欧美在线999| 一区二区亚洲| 亚洲一区二区三区涩| 国产精品久久久久aaaa樱花| 欧美成人第一页| 蜜桃伊人久久| 久久久夜精品| 久久久久久久91| 午夜激情综合网| 亚洲欧美视频在线观看视频| 国产精品sss| 久久精品国产精品亚洲精品| 亚洲激情综合| 欧美午夜免费影院| 国产午夜亚洲精品理论片色戒| 亚洲欧洲一级| 国产精品扒开腿爽爽爽视频| 欧美成人午夜剧场免费观看| 欧美福利电影网| 欧美r片在线| 欧美经典一区二区| 亚洲黄网站在线观看| 国产精品jizz在线观看美国| 亚洲精品在线观看免费| 国产精品久久福利| 亚洲视频网站在线观看| 国产视频亚洲精品| 久久国产日本精品| 亚洲免费av网站| 欧美日韩1区2区| 先锋影音久久| 国产香蕉97碰碰久久人人| 女人香蕉久久**毛片精品| 樱桃成人精品视频在线播放| 亚洲电影第三页| 亚洲欧洲在线看| 亚洲午夜视频在线观看| 欧美亚洲尤物久久| 欧美韩国日本一区| 一区二区三区精密机械公司 | 亚洲综合好骚| 午夜精品久久久久久99热软件| 久久精品免费观看| 欧美激情a∨在线视频播放| 亚洲欧洲一区二区三区| 国产欧美一级| 久久国产免费看| 国产精品亚洲产品| 欧美一级夜夜爽| 一区二区免费在线播放| 国产精品乱看| 亚洲国产精品久久久久| 日韩亚洲在线| 久久激情视频久久| 欧美激情第10页| 国产日韩欧美一区二区三区在线观看 | 欧美日韩中文字幕在线| 亚洲精品美女在线观看| 依依成人综合视频| 另类天堂视频在线观看| 欧美在线精品免播放器视频| 国产香蕉97碰碰久久人人| 欧美一级网站| 亚洲国产裸拍裸体视频在线观看乱了中文 | 欧美风情在线观看| 亚洲网站在线| 亚洲精品永久免费| 国产欧美日韩中文字幕在线| 亚洲国产精品日韩| 欧美一区二区免费视频| 噜噜噜久久亚洲精品国产品小说| 99精品99| 亚洲美女中文字幕| 国产精品一区2区| 亚洲免费观看高清在线观看 | 久久久女女女女999久久| 欧美午夜精品久久久久免费视 | 久久大逼视频| 国产欧美视频一区二区| 久久久久久久久蜜桃| 亚洲国产精品成人精品| 国产伦精品一区二区三区高清版| 久久米奇亚洲| 亚洲欧美偷拍卡通变态| 久久久久欧美精品| 亚洲欧美大片| 日韩写真在线| 亚洲第一精品电影| 夜久久久久久| 欧美激情视频在线播放| 国产亚洲一区在线播放| 国产精品v日韩精品| 久久精品91| 亚洲欧美日韩区| 久久久人成影片一区二区三区| 亚洲资源在线观看| 一区二区欧美亚洲| 亚洲精品乱码久久久久| 亚洲第一区中文99精品| 久久精品五月婷婷| 国产精品国产三级欧美二区| 欧美国产日韩一区二区三区| 女女同性精品视频| 国产精品99一区| 欧美成人激情视频| 这里只有精品在线播放| 亚洲人成艺术| 亚洲一区在线播放| 国产日韩欧美电影在线观看| 久久久久国产精品人| 亚洲国产欧美日韩| 午夜精品久久久久久久男人的天堂| 国产欧美日韩综合一区在线观看 | 韩国福利一区| 午夜精品美女久久久久av福利| 久久久天天操| 久久精品国产欧美激情| 在线精品视频一区二区三四| 狠狠色丁香婷婷综合| 免费91麻豆精品国产自产在线观看| 久久久久一区二区三区| 亚洲国产欧美一区| 亚洲国产精品久久久久秋霞蜜臀 | 9l国产精品久久久久麻豆| 性色av一区二区怡红| 欧美在线一二三区| 亚洲国产日韩一区二区| 日韩一区二区免费看| 国产精品一区二区三区久久久| 国产女人18毛片水18精品| 久久久久久久尹人综合网亚洲| 小辣椒精品导航|