codeforces 932E teamwork

         题意:求ΣC(n,i)*i^k。(i=1~n)

         据称有高等数学常识的人都可以看出这是从某一个式子求导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));
} 


  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值