标签：Yoi int sum 乐队 ret choose 350 brace ll

题面

求：\(\sum_{i=0}^{n} f(i)2^{n-i}{n \choose i}\),其中\(f(x)\)是一个\(K\)次多项式,\(n\leq 10^9,K\leq 5000\)

分析

关键在于: \(i^j=\sum_{h=0}^{i}{n \choose i}{j \brace i}\)
证明：考虑意义：
LHS表示将i个不同的球放入j个不同的盒子中的方案数
RHS表示枚举有多少空盒子，用第二类斯特拉林数计算对应的方案数
\( \begin {aligned}\sum_{i=0}^{n} f(i)*2^{n-i}*{n \choose i} &= \sum_{i=0}^{n}\sum_{j=0}^{k}a_ji^j2^{n-i}{n\choose i} \\&=\sum_{j=0}^{k}a_j\sum_{i=0}^{n} i^j{n \choose i}2^{n-i} \\&=\sum_{j=0}^{k} a_j\sum_{i=0}^{n}\sum_{h=0}^{j}{j \brace h}{i \choose h} h! {n\choose i}2^{n-i} \\&=\sum_{j=0}^{k} a_j\sum_{h=0}^{j}{j \brace h}h!\sum_{i=0}^{n}{i \choose h}{n \choose i}2^{n-i} \\ &=\sum_{j=0}^{k} a_j\sum_{h=0}^{j}{j \brace h}h!{n \choose h}\sum_{i=0}^{n}{n-h \choose i-h}2^{n-i} \\&=\sum_{j=0}^{k} a_j\sum_{h=0}^{j}{j \brace h}h!{n \choose h}3^{n-h} \end {aligned} \)

#include<bits/stdc++.h>
#define ll long long
const int p=998244853;
using namespace std;

inline int ksm(ll a,int b) {
	ll ret=1;
	while(b) {
		if(b&1) ret=ret*a%p;
		a=a*a%p,b>>=1;
	}
	return ret;
}
const int N=5005;
int d[N][N],n,K,a[N],s[N];

int main() {
	freopen("b.in","r",stdin);
	freopen("b.out","w",stdout);
	scanf("%d%d",&n,&K);
	for(int i=0;i<=K;i++) scanf("%d",&a[i]);
	d[0][0]=1;
	for(int i=1;i<=K;i++) {
		for(int j=1;j<=i;j++) {
			d[j][i]=((ll)d[j][i-1]*j+d[j-1][i-1])%p;
		}
	}
	int jc=1,mi=ksm(3,n),inv3=ksm(3,p-2)%p,c=1;
	s[0]=mi;
	for(int i=1;i<=K;i++) {
		c=(ll)c*(n-i+1)%p*ksm(i,p-2)%p;
		mi=(ll)mi*inv3%p;
		jc=(ll)jc*i%p;
		s[i]=(ll)c*jc%p*mi%p;
	}
	int ans=0;
	for(int i=0;i<=K;i++) {
		ll sum=0;
		for(int j=0;j<=i;j++) {
			sum=(sum+(ll)d[j][i]*s[j])%p;
		}
		ans=((ll)ans+sum*a[i])%p;
	}
	printf("%d\n",ans);
	return 0;
}

标签：Yoi,int,sum,乐队,ret,choose,350,brace,ll
来源： https://www.cnblogs.com/wsfwsf/p/13903223.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。