据称有高等数学常识的人都可以看出这是从某一个式子求导k次得到的。
我们设定义f(x)=(1+x)^n,即为ΣC(n,i)*(x^i)(i=0~n),
对这个式子进行一次求导,得到f'(x)=n*(1+x)^(n-1) 即 n*ΣC(n-1,i)*(x^i) (i=0~n-1),
因为c(n-1,i)=(n-1)!/(i!*(n-1-i)!)
我们将外面的n乘进去,再给每个c(n-1,i)除以 (i+1) 可将c(n-1,i)变为n!/((i+1)! * (n-(i+1))!)),
那么就将原式变为了ΣC(n,i)*(x^(i-1))*i (i=1~n),成功引入了i的一次方
再给该式子乘x,就变成了ΣC(n,i)*(x^i)*i(i=1~n)
那么我们我们对(1+x)^n次方求k次导,每求导一步之后再乘x,最后x=1时的解就是答案了。
你觉得变成x*n*(1+x)^(n-1)之后,再求导答案就变了?我们再求一次导。
对它再求导乘x,可以得到x*n*(1+x)^n-1 + x*x*n*(n-1)*(1+x)^(n-2)
我们对加号右边的式子使用类似刚刚的整理方法,可以得到其等于ΣC(n,i)*(x^i)*(i-1)*i(i=1~n)
再对加号左边的式子整理得到,得到了ΣC(n,i)*(x^i)*i(i=1~n)
两个相加刚好得到了k=2时候的答案。
所以可以证明每次求导的时候,左边的式子会把右边式子多余的补上,从而维持像题目的式子即可。
最后我们用dp[i][j]表示执行了i次操作,化出来的式子x的指数为j的答案,因为每次求导x的指数与(1+x)的指数和保持为n,所以j+k==n 而且j<=i,所以由高中数学可以得到转移方程(前不导乘后导,后不导乘前导)dp[a][b]=(b*dfs(a-1,b)%mod+(n-b)*dfs(a-1,b+1)%mod)%mod,记忆化搜索即可。
下附AC代码。
#include<iostream>
#include<stdio.h>
#include<string.h>
#include<algorithm>
#define maxn 5005
using namespace std;
typedef long long ll;
const ll mod=1e9+7;
ll quickpow(ll p,ll k)
{
ll ans=1;
while(k)
{
if(k&1)
ans=(ans*p)%mod;
p=(p*p)%mod;
k>>=1;
}
return ans;
}
ll n,k;
ll dp[maxn][maxn];
ll dfs(ll a,ll b)
{
if(b>n) return 0;
if(dp[a][b]!=-1) return dp[a][b];
if(a==0) return dp[a][b]=quickpow(2,n-b);
return dp[a][b]=(b*dfs(a-1,b)%mod+(n-b)*dfs(a-1,b+1)%mod)%mod;
}
int main()
{
memset(dp,-1,sizeof(dp));
scanf("%I64d%I64d",&n,&k);
printf("%I64d\n",dfs(k,0));
}