讀半天不解其意,傻算了半天才把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;
}