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

            Push Botton Lock poj 3088斯特靈數

            讀半天不解其意,傻算了半天才把sample湊出來.其實很簡單,連高精度都沒有.第一類斯特靈數S(n,m)就是把n元集合分成m部的個數,有遞推關系S(n,m)=S(n-1,m-1)+mS(n-1,m).所求還要全排列一下.再乘以m!就可以了累加1~B個數全部用上,就是結果.F(B)=sum(C(B,i)*(Sum(Stir( i,j )* j ! ) ) )就是結果.其中下標i從1變到B,j從1變到i.
            /*Source CodeProblem: 3088        User: y09shendazhi
            Memory: 164K        Time: 0MS
            Language: C++        Result: Accepte/

            Source Code
            */

            #include 
            <iostream>
            using namespace std;
            // /n\
            // | |
            // \m/
            __int64 comb(__int64 m,__int64 n)
            {
                
            if(n<m)
                    
            return 0;
                __int64 result
            =1;
                m
            =m<n-m?m:n-m;
                
            for(__int64 i=1;i<=m;i++)
                    result
            =(result*(n-m+i))/i;    
                
            return result;
            }

            __int64 fun(__int64 n)
            {
                __int64 ans
            =1;
                
            for(__int64 i=1;i<=n;i++)
                    ans
            *=i;
                
            return ans;
            }

            int main(__int64 argc, char *argv[])
            {
                __int64 i,j,k;
                __int64 stir[
            15][15]={0};
                
            for(i=1;i<12;i++)
                
            {
                    stir[i][
            1]=1;
                    
            for(j=2;j<i;j++)
                    
            {
                        stir[i][j]
            =stir[i-1][j-1]+j*stir[i-1][j];
                    }

                    stir[i][i]
            =1;
                }

                __int64 ans[
            12]={0,1};
                
            for(i=2;i<12;i++)
                
            {
                    
                    
            for(j=1;j<=i;j++)
                    
            {
                        __int64 temp
            =0;
                        
            for(k=1;k<=j;k++)
                        
            {
                            temp
            +=stir[j][k]*fun(k);
                        }

                        ans[i]
            +=temp*comb(j,i);
                    }

                }

                __int64 
            in=0;
                __int64 t
            =0;
                scanf(
            "%d",&t);
                
            for(i=0;i<t;i++)
                
            {
                    scanf(
            "%d",&in);
                    printf(
            "%I64d %I64d %I64d\n",i+1,in,ans[in]);
                }



                
                
            return 0;
            }

            posted on 2010-08-21 18:20 若余 閱讀(390) 評論(0)  編輯 收藏 引用

            導航

            <2010年8月>
            25262728293031
            1234567
            891011121314
            15161718192021
            22232425262728
            2930311234

            統計

            常用鏈接

            留言簿

            隨筆檔案(16)

            搜索

            最新隨筆

            最新評論

            評論排行榜

            久久九九兔免费精品6| 99久久精品免费看国产一区二区三区 | 久久久久人妻一区精品性色av| 日韩乱码人妻无码中文字幕久久| 国内精品久久久久| 欧美与黑人午夜性猛交久久久| 久久久久久久97| 久久人人爽人人精品视频| 久久夜色精品国产噜噜亚洲AV| 久久99精品久久久久久9蜜桃| 色诱久久久久综合网ywww | 久久久久久久久波多野高潮| 国产一久久香蕉国产线看观看| 亚洲人AV永久一区二区三区久久| 久久综合久久综合九色| 久久久亚洲裙底偷窥综合| 无码8090精品久久一区| 免费精品99久久国产综合精品| 亚洲精品无码专区久久久| 三级韩国一区久久二区综合 | 新狼窝色AV性久久久久久| 国产成人综合久久精品尤物| 国产Av激情久久无码天堂| 久久婷婷五月综合国产尤物app| 精品久久久久久无码人妻蜜桃| 久久亚洲AV成人出白浆无码国产| 日产精品久久久久久久| 亚洲国产成人精品无码久久久久久综合 | 中文字幕热久久久久久久| 欧美粉嫩小泬久久久久久久| 99久久综合国产精品二区| 91精品国产高清91久久久久久 | 精品久久久久久无码中文字幕 | 香蕉久久夜色精品升级完成| 精产国品久久一二三产区区别| 午夜精品久久久久成人| 一级做a爰片久久毛片免费陪| 亚洲国产成人精品久久久国产成人一区二区三区综 | 国产精品无码久久久久| 国产精品久久久久乳精品爆 | 久久人搡人人玩人妻精品首页|