问题:设
T(n,k)=nk
T
(
n
,
k
)
=
n
k
,
S(n,k)=∑ni=1T(i)
S
(
n
,
k
)
=
∑
i
=
1
n
T
(
i
)
。给出n和k,求S(n)。
例如
k=2,n=5,S(n,k)=12+22+32+42+52=55
k
=
2
,
n
=
5
,
S
(
n
,
k
)
=
1
2
+
2
2
+
3
2
+
4
2
+
5
2
=
55
。
由于结果很大,输出
S(n)MOD1000000007
S
(
n
)
M
O
D
1000000007
的结果即可。
伯努利数是18世纪瑞士数学家雅各布·伯努利引入的一个数。设伯努利数为
Bn
B
n
,它的定义为:
tet−1=∑ni=1Bnn!∗tn
t
e
t
−
1
=
∑
i
=
1
n
B
n
n
!
∗
t
n
这里|t|<2 。由计算知:
B0=1,B1=−12
B
0
=
1
,
B
1
=
−
1
2
一般地,当
n≥2
n
≥
2
时,我们可以通过
∑ni=0C(n+1,i)∗Bk=0
∑
i
=
0
n
C
(
n
+
1
,
i
)
∗
B
k
=
0
通过这个公式我们就得到了
Bn
B
n
的公式:
**
Bn=−1n+1∗(C0n+1∗B0+C1n+1∗B1+...+Cn−1n+1∗Bn−1)) B n = − 1 n + 1 ∗ ( C n + 1 0 ∗ B 0 + C n + 1 1 ∗ B 1 + . . . + C n + 1 n − 1 ∗ B n − 1 ) )
**
那么 我现在给出 求
1k+2k+...+nk
1
k
+
2
k
+
.
.
.
+
n
k
的关于伯努利的公式:
∑ni=1ik=1k+1∑k+1i=1Cik+1∗B[k+1−i]∗(n+1)i ∑ i = 1 n i k = 1 k + 1 ∑ i = 1 k + 1 C k + 1 i ∗ B [ k + 1 − i ] ∗ ( n + 1 ) i
通过这个式子就可以在 O(K) O ( K ) 的复杂度求出 S(n,k) S ( n , k )
#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int MOD=1e9+7;
const int N=1010;
int c[N][N],inv[N];
LL B[N];
void init(int n)
{
c[0][0]=1;
for(int i=1;i<=n;i++)
{ c[i][0]=c[i][i]=1;
for(int j=1;j<i;j++)
c[i][j]=(c[i-1][j-1]+c[i-1][j])%MOD;
}
inv[1]=1;
for(int i=2;i<=n;i++)
inv[i]=(LL)(MOD-MOD/i)*inv[MOD%i]%MOD;
B[0]=1;
for(int i=1;i<=n;i++)
{
for(int j=0;j<i;j++)
{
B[i]=(B[j]*c[i+1][j]%MOD+B[i])%MOD;
}
B[i]=(B[i]*(-inv[i+1])%MOD+MOD)%MOD;
}
}
// B[100]:94103270
LL Bnum(int n,int k)
{
LL ans=0,tem=1;
for(int i=1;i<=k+1;i++)
{
tem=tem*(n+1)%MOD;
ans=(c[k+1][i]*B[k+1-i]%MOD*tem%MOD+ans)%MOD;
}
ans=(ans*inv[k+1]%MOD+MOD)%MOD;
return ans;
}
LL pod(LL x,LL n)
{
LL ret=1;
while(n)
{
if(n&1)ret=ret*x%MOD;
x=x*x%MOD;
n>>=1;
}
return ret;
}
int main()
{
init(1005);
int n,k;
while(scanf("%d%d",&n,&k)!=EOF)cout<<Bnum(n,k)<<endl;
}