BZOJ2159: Crash 的文明世界

这篇写差分表和斯特林数介绍的不错

---

这题就是要计算这个东西

$$S(i)=\sum_{j=1}^ndist(i,j)^k$$

这个东西很难维护，我们把$dist(i,j)^k$拆一下

$S(u)=\sum_v\sum_{j=0}^k\{^k_j\}\frac{d(u,v)!}{(d(u,v)-j)!}$

$S(u)=\sum_v\sum_{j=0}^k\{^k_j\}j!(^{d(u,v)}_j)$

$S(u)=\sum_{j=0}^k\{^k_j\}j!\sum_v(^{d(u,v)}_j)$

令f[u][j]表示u得子树内的$\sum_v(^{d(u,v)}_j)$，因为$(^i_j)=(^{i-1}_j)+(^{i-1}_{j-1})$，可以用f[v][j-1]+f[v][j]转移到f[u][j]这样转移

code：

#include<set>
#include<map>
#include<deque>
#include<queue>
#include<stack>
#include<cmath>
#include<ctime>
#include<bitset>
#include<string>
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<climits>
#include<complex>
#include<iostream>
#include<algorithm>
#define ll long long
using namespace std;

const int mod = 10007;
const int maxn = 51000;
const int maxk = 210;
inline void add(int &a,const int &b){a+=b;if(a>=mod)a-=mod;}
inline void dec(int &a,const int &b){a-=b;if(a<0)a+=mod;}

int s[maxk],str[maxk][maxk];
void pre()
{
    s[0]=1; for(int i=1;i<maxk;i++) s[i]=s[i-1]*i%mod;
    str[0][0]=1;
    for(int i=1;i<maxk;i++)
    {
        for(int j=1;j<=i;j++) str[i][j]=(str[i-1][j-1]+str[i-1][j]*j%mod)%mod;
    }
}

int n,K;
struct edge{int y,nex;}a[maxn<<1]; int len,fir[maxn];
inline void ins(const int x,const int y){a[++len]=(edge){y,fir[x]};fir[x]=len;}
void Read()
{
    int now,A,B,Q,L;
    scanf("%d%d%d",&n,&K,&L);
    scanf("%d%d%d%d",&now,&A,&B,&Q);
    for(int i=1;i<n;i++)
    {
        now=(now*A+B)%Q; int tmp=i<L?i:L;
        int x=i-now%tmp,y=i+1;
        ins(x,y); ins(y,x);
    }
    /*scanf("%d%d",&n,&K);
    for(int i=1;i<n;i++)
    {
        int x,y; scanf("%d%d",&x,&y);
        ins(x,y); ins(y,x);
    }*/
}

int f[maxn][maxk],g[maxn][maxk];
int temp[maxk];
int fa[maxn];
void dpf(const int x)
{
    f[x][0]=1;
    for(int k=fir[x],y=a[k].y;k;k=a[k].nex,y=a[k].y) if(y!=fa[x])
    {
        fa[y]=x;
        dpf(y);
        add(f[x][0],f[y][0]);
        for(int j=1;j<=K;j++) add(f[x][j],(f[y][j-1]+f[y][j])%mod);
    }
}
void dpg(const int x)
{
    for(int j=0;j<=K;j++) temp[j]=f[x][j],add(temp[j],g[x][j]);
    for(int k=fir[x],y=a[k].y;k;k=a[k].nex,y=a[k].y) if(y!=fa[x])
    {
        dec(temp[0],f[y][0]);
        for(int j=1;j<=K;j++) dec(temp[j],(f[y][j-1]+f[y][j])%mod);

        g[y][0]=temp[0];
        for(int j=1;j<=K;j++) g[y][j]=(temp[j-1]+temp[j])%mod;

        add(temp[0],f[y][0]);
        for(int j=1;j<=K;j++) add(temp[j],(f[y][j-1]+f[y][j])%mod);
    }
    for(int k=fir[x],y=a[k].y;k;k=a[k].nex,y=a[k].y) if(y!=fa[x])
        dpg(y);
}

void print()
{
    for(int i=1;i<=n;i++)
    {
        int ans=0;
        for(int j=0;j<=K;j++) add(ans,str[K][j]*s[j]%mod*((f[i][j]+g[i][j])%mod)%mod);
        printf("%d\n",ans);
    }
}

int main()
{
    //freopen("tmp.in","r",stdin);
    //freopen("tmp.out","w",stdout);

    pre();

    Read();

    dpf(1);
    dpg(1);

    print();

    return 0;
}