BZOJ5250: [2018多省省隊聯測]祕密襲擊-樹形DP

傳送門

題意：

給一棵n個點的樹，每個點的點權在 1到 W之間

求所有連通塊的權值第k大的和模 64123

k≤n≤1666,W≤1666

Solution：

正解貌似是線段樹合併+FFT 但是我並不會寫QAQ

所以說我們考慮暴力碾標算：

我們可以考慮每個點對於答案的貢獻：

我們把大於它的點看成1，小於它的點看成0，最後只要求包含它的和爲k-1的聯通塊個數即可，這個可以用一個樹形DP來實現

複雜度爲O(n2k)$O(n^{2}k)$

加一些剪枝就可以過了而且跑得比正解快

代碼：

#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
const int N=2010;
const int mod=64123;
int n,k,w,S;
int d[N],head[N],size;
int f[2010][2010],ans;
struct edg{
    int to,next;
}e[2*N];
void add(int x,int y)
{
    size++;e[size].to=y;e[size].next=head[x];head[x]=size; 
}
void dfs(int x,int fa)
{
    if ((d[S]<d[x])||((d[x]==d[S])&&x<S))
        for (int i=1;i<k;i++) {f[x][i+1]+=f[fa][i];if (f[x][i+1]>=mod) f[x][i+1]-=mod;}
    else for (int i=1;i<=k;i++) {f[x][i]+=f[fa][i];if (f[x][i]>=mod) f[x][i]-=mod;}
    for (int i=head[x];i;i=e[i].next)
    {
        int y=e[i].to;
        if (fa==y) continue;
        dfs(y,x);
    }
    for (int i=1;i<=k;i++) {f[fa][i]+=f[x][i];if (f[fa][i]>=mod) f[fa][i]-=mod;}
}
int main()
{
    scanf("%d%d%d",&n,&k,&w);
    for (int i=1;i<=n;i++) scanf("%d",&d[i]);
    for (int x,y,i=1;i<n;i++) scanf("%d%d",&x,&y),add(x,y),add(y,x);
    for (int i=1;i<=n;i++)
    {
        S=i;int tot=0;
        for (int j=1;j<=n;j++) if ((d[j]>d[i])||((d[i]==d[j])&&i>j)) tot++;
        if (tot<k-1) continue;
        memset(f,0,sizeof(f));
        f[i][1]=1;
        for (int j=head[i];j;j=e[j].next)
        {
            int y=e[j].to;
            dfs(y,i);
        }
        ans+=(1ll*f[i][k]*d[i])%mod;
        if (ans>=mod) ans-=mod;
    }
    printf("%d",ans);
}