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

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>
            在线国产日韩| 欧美日韩成人网| 国产视频在线观看一区二区三区 | 性欧美xxxx视频在线观看| 亚洲日本成人网| 欧美精品黄色| 亚洲综合精品自拍| 亚洲免费视频网站| 国产亚洲欧洲997久久综合| 久久欧美中文字幕| 免费的成人av| 亚洲视频精选| 亚洲欧美区自拍先锋| 精品91视频| 99精品国产99久久久久久福利| 国产精品久久一卡二卡| 久久精品国产清高在天天线| 老鸭窝毛片一区二区三区| 一区二区免费在线播放| 欧美伊人久久| 亚洲精品免费观看| 亚洲欧美视频一区二区三区| 在线观看视频欧美| 日韩视频在线观看一区二区| 日韩一级精品视频在线观看| 久久综合色播五月| 亚洲一二三区在线观看| 欧美在线短视频| 一本一本久久a久久精品综合麻豆| 亚洲伊人网站| 亚洲国产精品va在线看黑人| 这里只有视频精品| 亚洲韩国精品一区| 亚洲午夜久久久久久久久电影院| 亚洲电影第三页| 亚洲欧美国产77777| 亚洲精品一区二区在线| 欧美一区二区三区四区在线观看地址 | 久久精品一区二区| 欧美日韩午夜精品| 久久夜色精品国产噜噜av| 欧美日一区二区在线观看| 美女黄网久久| 国产一区二区三区奇米久涩| 一区二区三欧美| 日韩视频永久免费观看| 久久人人97超碰人人澡爱香蕉| 午夜久久一区| 欧美色精品在线视频| 亚洲国产精品日韩| 伊大人香蕉综合8在线视| 亚洲制服av| 午夜精品偷拍| 国产精品激情| 亚洲五月六月| 亚洲一区二区三区精品在线观看 | 国产综合婷婷| 亚洲制服少妇| 久久av资源网| 国产麻豆综合| 午夜欧美不卡精品aaaaa| 香蕉成人伊视频在线观看| 国产精品mm| 亚洲视频在线看| 亚洲永久免费| 国产精品亚洲综合色区韩国| 亚洲一品av免费观看| 香蕉av777xxx色综合一区| 国产日韩欧美二区| 欧美一区二区三区婷婷月色 | 午夜伦理片一区| 国产麻豆日韩| 久久美女艺术照精彩视频福利播放| 久久免费黄色| 最近中文字幕日韩精品 | 午夜精品福利在线| 久久久国产视频91| 在线成人欧美| 欧美成年人网| 日韩亚洲欧美一区二区三区| 亚洲欧美99| 久久亚洲精品一区二区| 欧美日韩高清不卡| 在线中文字幕一区| 久久精品视频导航| 在线观看欧美精品| 欧美精品国产精品| 亚洲欧美成人网| 欧美成人影音| 亚洲一区在线免费观看| 国产亚洲毛片在线| 榴莲视频成人在线观看| 日韩一区二区精品| 久久精品99国产精品日本| 亚洲国产精品一区| 国产精品v片在线观看不卡| 欧美中文字幕视频| 亚洲人成小说网站色在线| 亚欧成人在线| 最新日韩在线视频| 国产精品天天看| 免费在线观看一区二区| 亚洲午夜一区二区三区| 欧美大尺度在线观看| 亚洲伊人伊色伊影伊综合网| 伊人久久大香线蕉综合热线 | 亚洲视频在线观看网站| 欧美wwwwww| 性一交一乱一区二区洋洋av| 亚洲精品欧美日韩| 国产啪精品视频| 欧美精品亚洲精品| 久久精品国产v日韩v亚洲| 亚洲精品免费在线| 欧美ed2k| 久久久精品一区| 亚洲一区二区三区在线看| 亚洲第一福利视频| 国产女主播在线一区二区| 欧美日韩1234| 老司机aⅴ在线精品导航| 亚洲免费一区二区| 99成人在线| 亚洲精品美女在线观看播放| 欧美电影免费观看高清完整版| 欧美一区观看| 香蕉免费一区二区三区在线观看 | 国产一区二区三区免费在线观看| 国产精品国产三级欧美二区| 欧美日韩国产一区精品一区| 蜜臀99久久精品久久久久久软件| 欧美在线地址| 欧美在线视屏 | 亚洲福利视频三区| 欧美波霸影院| 欧美成人精品一区二区三区| 久久夜色撩人精品| 久久精品视频导航| 久久噜噜亚洲综合| 久久女同精品一区二区| 久久噜噜噜精品国产亚洲综合| 久久久www| 久久一区精品| 欧美91视频| 亚洲高清成人| 亚洲人成久久| 一本色道久久88亚洲综合88| 夜夜精品视频一区二区| 一区二区欧美在线| 亚洲一区二区三区精品视频| 99精品国产在热久久婷婷| 一区二区三区四区国产精品| 中文成人激情娱乐网| 亚洲欧美日韩久久精品| 欧美日韩另类视频| 久久青草福利网站| 久久综合九色九九| 欧美精品电影在线| 欧美日韩在线不卡| 国产精品乱码久久久久久| 国产精品一区二区视频| 国内自拍亚洲| 最新亚洲激情| 午夜久久久久久| 麻豆精品视频在线观看| 亚洲国产午夜| 亚洲视频二区| 久久久亚洲午夜电影| 欧美乱人伦中文字幕在线| 欧美视频一区二区三区| 国产一区视频观看| 亚洲精品一区在线| 性欧美暴力猛交69hd| 欧美成人国产va精品日本一级| 亚洲另类一区二区| 亚洲欧美高清| 欧美福利视频一区| 国产免费成人| 亚洲日韩视频| 欧美专区亚洲专区| 欧美激情一区二区久久久| 亚洲色图自拍| 蜜桃av噜噜一区二区三区| 国产精品区一区二区三区| 亚洲国产成人久久| 欧美一区二区三区视频| 亚洲电影天堂av| 欧美一区二区精品在线| 欧美激情亚洲自拍| 国产综合久久久久久| 亚洲一区二区三区中文字幕| 欧美激情精品久久久久久蜜臀| 亚洲影视九九影院在线观看| 欧美成人按摩| 黄网动漫久久久| 久久av红桃一区二区小说| 中文日韩在线| 欧美日韩成人在线播放| 亚洲国产精品一区二区第一页 | 亚洲国产欧美在线|