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

misschuer

常用鏈接

統(tǒng)計(jì)

積分與排名

百事通

最新評(píng)論

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;
}
這題其實(shí)是求(n!)^2的約數(shù)的個(gè)數(shù),注意要大數(shù)處理。其余的自己去算吧,orz~~~~~~

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


只有注冊(cè)用戶(hù)登錄后才能發(fā)表評(píng)論。
網(wǎng)站導(dǎo)航: 博客園   IT新聞   BlogJava   博問(wèn)   Chat2DB   管理


青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品
  • <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>
            亚洲人线精品午夜| 国产亚洲精品7777| av72成人在线| 99天天综合性| 亚洲一区二区三区精品动漫| 一区二区三区三区在线| 亚洲性感美女99在线| 午夜精品久久久久久久久久久久久 | 亚洲人线精品午夜| 欧美成人a∨高清免费观看| 国产精品视频免费| 国产精品日韩欧美一区二区| 国产一区二区三区精品欧美日韩一区二区三区 | 免费国产一区二区| 欧美日韩成人综合天天影院| 国产精品美女久久| 韩国女主播一区二区三区| 亚洲国产欧美日韩精品| 在线成人性视频| 亚洲欧洲日夜超级视频| 毛片一区二区三区| 欧美精品观看| 国产日韩精品视频一区二区三区| 国内精品嫩模av私拍在线观看| 亚洲福利专区| 亚洲一区二区欧美| 久久伊人免费视频| 中文欧美字幕免费| 欧美www在线| 国产欧美精品xxxx另类| 亚洲激情在线观看| 久久激情视频久久| 亚洲精品欧美一区二区三区| 欧美在线免费播放| 国产精品成人免费| 红桃视频国产一区| 夜夜躁日日躁狠狠久久88av| 日韩亚洲欧美一区二区三区| 亚洲美女免费精品视频在线观看| 欧美国产欧美亚洲国产日韩mv天天看完整 | 欧美激情第4页| 国产一区二区三区久久精品| 一区二区三区不卡视频在线观看 | 午夜精品久久久久| 久久久久久免费| 99视频有精品| 欧美激情小视频| 在线观看成人小视频| 西瓜成人精品人成网站| 亚洲开发第一视频在线播放| 免费在线看成人av| 亚洲第一伊人| 老司机精品久久| 久久精品观看| 激情久久一区| 久久这里有精品视频| 午夜精品剧场| 国产一区二区日韩精品欧美精品| 亚洲伊人伊色伊影伊综合网| 亚洲第一精品夜夜躁人人爽| 午夜欧美不卡精品aaaaa| 久久久久五月天| 免播放器亚洲一区| 亚洲电影第三页| 欧美jjzz| 免费不卡在线观看| 亚洲人成网站精品片在线观看| 裸体素人女欧美日韩| 久久精品二区| 亚洲国产精品精华液网站| 欧美77777| 欧美精品videossex性护士| av成人免费观看| 亚洲天堂av高清| 国产三区二区一区久久| 欧美伊人久久久久久午夜久久久久| 亚洲一区国产视频| 国产伦精品一区二区三区视频孕妇 | 亚洲欧美日韩精品一区二区| 国产精品久久久久天堂| 午夜精品久久久久久久蜜桃app| 亚洲一区二区三区影院| 国产亚洲电影| 欧美成人tv| 欧美日韩视频免费播放| 新67194成人永久网站| 欧美在线网站| 亚洲精品视频在线| 亚洲性av在线| 91久久久一线二线三线品牌| 在线亚洲一区二区| 国内精品久久久久国产盗摄免费观看完整版 | 欧美大片免费久久精品三p| 一本色道**综合亚洲精品蜜桃冫| 一区二区免费在线视频| 国产亚洲人成网站在线观看| 亚洲高清在线视频| 国产精品乱码人人做人人爱| 久久一区二区三区四区五区| 欧美日本高清视频| 久久国产福利国产秒拍| 欧美极品欧美精品欧美视频| 欧美一区国产在线| 欧美韩日亚洲| 久久久久久久精| 欧美三级欧美一级| 欧美aⅴ一区二区三区视频| 国产精品大全| 亚洲电影在线看| 国产一区日韩一区| 一本大道av伊人久久综合| 在线不卡中文字幕播放| 亚洲少妇在线| 99riav国产精品| 久久精品噜噜噜成人av农村| 亚洲综合色网站| 欧美理论电影在线播放| 久久综合九色九九| 国产日韩精品一区| 亚洲永久精品大片| 亚洲小少妇裸体bbw| 欧美福利视频| 亚洲福利av| 亚洲精品永久免费| 欧美gay视频激情| 欧美.日韩.国产.一区.二区| 国产综合香蕉五月婷在线| 亚洲欧美日本日韩| 午夜欧美不卡精品aaaaa| 国产精品成人一区二区三区夜夜夜 | 欧美日韩国产一区精品一区| 欧美成人在线免费观看| 激情欧美丁香| 另类人畜视频在线| 国产一区二区三区免费观看 | 久久国产加勒比精品无码| 久久精品国产久精国产思思| 国产老女人精品毛片久久| 在线视频日韩精品| 亚洲综合电影一区二区三区| 国产精品国产精品| 亚洲中字在线| 久久久蜜桃一区二区人| 激情综合自拍| 免费观看日韩| 亚洲欧洲精品一区二区三区波多野1战4 | 欧美一级夜夜爽| 久久久在线视频| 亚洲国产免费看| 欧美精品自拍偷拍动漫精品| 99热在这里有精品免费| 亚洲欧美色婷婷| 国户精品久久久久久久久久久不卡| 欧美在线地址| 亚洲大片av| 亚洲一区二区在线免费观看视频| 国产精品伦一区| 久久久久九九九| 亚洲人成在线播放| 午夜宅男久久久| 亚洲激情婷婷| 国产精品一区二区久久久久| 久久精品国产一区二区三| 亚洲福利一区| 午夜精品视频在线| 在线成人免费视频| 欧美午夜不卡| 久久久久综合网| 亚洲精品一区二区三区四区高清 | 亚洲视频一二三| 国产综合色在线| 欧美日韩日本网| 久久免费视频网站| 亚洲视频一区二区免费在线观看| 久久午夜精品| 亚洲——在线| 亚洲精品美女| 国产在线精品一区二区夜色| 欧美日韩国产亚洲一区| 久久不射网站| 在线中文字幕一区| 亚洲国产精品t66y| 久久久蜜桃一区二区人| 亚洲在线一区二区三区| 亚洲激情网站| 樱桃国产成人精品视频| 国产精品毛片va一区二区三区 | 亚洲六月丁香色婷婷综合久久| 国产精品理论片在线观看| 欧美黄色成人网| 久久久久久高潮国产精品视| 亚洲无玛一区| 日韩视频永久免费| 亚洲国产高清高潮精品美女| 久久久久久久久久久久久久一区| 亚洲综合色激情五月| 一本久道久久综合狠狠爱| 亚洲精品久久久久久一区二区| 国产在线播放一区二区三区| 国产精品亚洲产品|