20210524 刷题赛 乐队_题解【另一种思路】

20210524 刷题赛 乐队_题解【另一种思路】

题解给的做法是第二类斯特林数通常幂转下降幂，其实更简单的也可以直接转成牛顿级数。

把 $F(x)$ 写成牛顿级数
$$F(x)=\sum _{j=0}^k{c}_j\left(\genfrac{}{}{0ex}{}{x}{j}\right)$$
目标式子为
$$\begin{array}{rl}& \sum _{i=0}^n\left(\genfrac{}{}{0ex}{}{n}{i}\right)2^{n-i}F(i)\\ =& \sum _{i=0}^n\sum _{j=0}^k\left(\genfrac{}{}{0ex}{}{n}{i}\right)2^{n-i}{c}_j\left(\genfrac{}{}{0ex}{}{i}{j}\right)\\ =& \sum _{j=0}^k\left(\genfrac{}{}{0ex}{}{n}{j}\right)c_j\sum _{i=0}^n{2}^{n-i}\left(\genfrac{}{}{0ex}{}{n-j}{i-j}\right)\\ =& \sum _{j=0}^k\left(\genfrac{}{}{0ex}{}{n}{j}\right)c_j{3}^{n-j}\end{array}$$
由于我太菜了，在最后一步想了好久，最后发现可以强行给 ${\sum }_i\left(\genfrac{}{}{0ex}{}{a}{i}\right)2^i$ 一个组合意义：

现有 $a$ 个白球，你选了 $i$ 个出来，并给这个 $i$ 个球任意涂上红色或蓝色，问有几种涂色方案。这个就等价于现有 $a$ 个球，可以给每个球任意涂上白红蓝三种颜色的一种，问方案数。

所以就有了 ${\sum }_i\left(\genfrac{}{}{0ex}{}{a}{i}\right)2^i=3^a$ ，这个可以扩展到 ${\sum }_i\left(\genfrac{}{}{0ex}{}{a}{i}\right)X^i=(X+1{)}^a$ 。

#include <bits/stdc++.h>
#define K 5003 
using namespace std;
typedef long long ll;
const int mod=998244853;
ll n,k,fac[K],inv[K];
ll a[K],A[K][K];
ll ksm(ll x,ll y){
    ll res=1;
    while(y){
        if(y&1) res=res*x%mod; 
        x=x*x%mod; y>>=1;
    }
    return res;
}
void init(){
    fac[0]=1;
    for(int i=1;i<=k;i++) fac[i]=fac[i-1]*i%mod;
    inv[k]=ksm(fac[k],mod-2);
    for(int i=k-1;i>=0;i--) 
        inv[i]=inv[i+1]*(i+1)%mod;
}
ll get_A(ll x){
    ll res=0,X=1;
    for(int i=0;i<=k;i++,X=X*x%mod) res=(res+X*a[i]%mod)%mod;
    return res;
}
ll C[K];
int main(){
//    freopen("band2.in","r",stdin);
    cin>>n>>k;
    init();
    for(int i=0;i<=k;i++) cin>>a[i],a[i]%=mod;
    for(int i=0;i<=k;i++) A[0][i]=get_A(i);
    for(int i=k,op=1;i>=1;i--,op++)
        for(int j=0;j<=i-1;j++)
            A[op][j]=(A[op-1][j+1]-A[op-1][j]+mod)%mod;
    for(int i=0;i<=k;i++) a[i]=A[i][0];
    C[0]=1;
    for(int i=1;i<=k;i++) C[i]=C[i-1]*(n-i+1)%mod;
    ll res=0;
    for(int i=0;i<=k;i++) 
		(res+=a[i]*C[i]%mod*inv[i]%mod*ksm(3,n-i)%mod)%=mod;
	cout<<res;
}