• <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>
            posts - 74,  comments - 33,  trackbacks - 0

            Xiaoming has just come up with a new way for encryption, by calculating the key from a publicly viewable number in the following way:
            Let the public key N = AB, where 1 <= A, B <= 1000000, and a0, a1, a2, …, ak-1 be the factors of N, then the private key M is calculated by summing the cube of number of factors of all ais. For example, if A is 2 and B is 3, then N = AB = 8, a0 = 1, a1 = 2, a2 = 4, a3 = 8, so the value of M is 1 + 8 + 27 + 64 = 100.
            However, contrary to what Xiaoming believes, this encryption scheme is extremely vulnerable. Can you write a program to prove it?

            Input

            There are multiple test cases in the input file. Each test case starts with two integers A, and B. (1 <= A, B <= 1000000). Input ends with End-of-File.
            Note: There are about 50000 test cases in the input file. Please optimize your algorithm to ensure that it can finish within the given time limit.

            Output

            For each test case, output the value of M (mod 10007) in the format as indicated in the sample output.

            Sample Input

            2 2
            1 1
            4 7

            Sample Output

            Case 1: 36
            Case 2: 1
            Case 3: 4393


            Author: 2008 Asia Hangzhou Regional Contest Online
            這道題是杭州賽區的網絡預選賽的賽題,當時打死也沒過,現在也沒有當時的代碼了,不知道哪里錯了,今天一看想起了一個公式就是
            1^3+2^3+……n^3=(n)^2*(n+1)*^2/4;記得這個公式還是高二老師教會的著,具體證明我們可以數學歸納法證明一下,簡單的很自己證明吧!
            而那時自己網絡賽的時候不知道怎么推導出過10007是個循環,可以把b的取模,而現在有了這個公式就可以一步搞定為什么對b取模!
            要解題就必須求出a ^b的因子數,這不是難點
            證明如下:
            N=a1^b1*a2^b2*……*aN^bN;根據排列組合原理我們知道N的因子數可以如下解出:
            即(b1+1)*(b2+1)*.......*(bN+1);而現在是a^b只需要把(b1*b+1)*(b2*b+1)*.......*(bN*b+1);極為因子數:
            而這個時候我們知道1<=a,b<=1000000超出整形范圍:我們進一步優化得到
            因為有:m*n%x=(m%x*n%x);(m+n)%x=(m%x+n%x);
            所以我們必須在過程中優化算法!
            這樣這道題目就出來了!當然求因子時必不可少素數打表!可以優化打表1000,減少時間!當然還可以把1-10007的密碼表打出,更優時間啊
            素數打表代碼如下:
            memset(prim,0,sizeof(prim));
            ????
            for(sign=0,i=2;i<MAXN;i++)
            ????????
            if(!prim[i]){
            ????????????num[sign
            ++]=i;
            ????????????
            for(j=0;i*j<MAXN;j++)
            ????????????????prim[i
            *j]=true;????
            ????????}



            posted on 2009-03-23 19:20 KNIGHT 閱讀(530) 評論(0)  編輯 收藏 引用
            <2009年3月>
            22232425262728
            1234567
            891011121314
            15161718192021
            22232425262728
            2930311234

            常用鏈接

            留言簿(8)

            隨筆檔案

            文章檔案

            Friends

            OJ

            搜索

            •  

            最新評論

            閱讀排行榜

            評論排行榜

            国产精品对白刺激久久久| 色综合色天天久久婷婷基地| 免费久久人人爽人人爽av| AV无码久久久久不卡蜜桃| 精品久久久久久亚洲| 久久综合成人网| 久久精品国产秦先生| 久久精品无码一区二区三区免费 | 国产精品无码久久四虎| 久久精品一区二区三区中文字幕| 久久99国产精品久久99小说| 国产午夜福利精品久久2021| 欧美亚洲日本久久精品| 伊人热人久久中文字幕| 久久人人爽人人爽人人片AV不| 99久久伊人精品综合观看| 2021国内久久精品| 久久久WWW成人免费精品| 精品久久一区二区三区| 国产成人无码精品久久久性色| 久久高潮一级毛片免费| 2021久久国自产拍精品| 18禁黄久久久AAA片| 欧美国产精品久久高清| 亚洲精品高清久久| 国产精品久久久久久久久| 色偷偷88888欧美精品久久久| 亚洲中文字幕伊人久久无码 | 中文字幕无码av激情不卡久久| 国产激情久久久久影院老熟女| 潮喷大喷水系列无码久久精品| 九九精品99久久久香蕉| 99久久精品免费看国产一区二区三区 | 久久久久久曰本AV免费免费| 婷婷国产天堂久久综合五月| 四虎影视久久久免费| 日批日出水久久亚洲精品tv| 久久久久99精品成人片牛牛影视| 国产精品九九久久免费视频 | 久久无码一区二区三区少妇| 国内精品久久久久久久涩爱|