2019CCPC湖南全国邀请赛 I Neko and tree 树形dp

dp[ i ][ j ] 表示在 i 这棵子树中, 距离 i 最远点的距离为 j 的方案数。

转移应该挺显然的把, 记一下mxd[ u ] 表示 u 这颗子树的最大深度, 然后两个for套在一起复杂度就对了。

#include<bits/stdc++.h>
#define LL long long
#define LD long double
#define ull unsigned long long
#define fi first
#define se second
#define mk make_pair
#define PLL pair<LL, LL>
#define PLI pair<LL, int>
#define PII pair<int, int>
#define SZ(x) ((int)x.size())
#define ALL(x) (x).begin(), (x).end()
#define fio ios::sync_with_stdio(false); cin.tie(0);

using namespace std;

const int N = 5000 + 7;
const int inf = 0x3f3f3f3f;
const LL INF = 0x3f3f3f3f3f3f3f3f;
const int mod = 1e9 + 7;
const double eps = 1e-8;
const double PI = acos(-1);

template<class T, class S> inline void add(T& a, S b) {a += b; if(a >= mod) a -= mod;}
template<class T, class S> inline void sub(T& a, S b) {a -= b; if(a < 0) a += mod;}
template<class T, class S> inline bool chkmax(T& a, S b) {return a < b ? a = b, true : false;}
template<class T, class S> inline bool chkmin(T& a, S b) {return a > b ? a = b, true : false;}

int f[N][N];
int g[N];
int c[N];
int n, m, k;
bool is[N];
int mxd[N];

vector<int> G[N];

void dfs(int u, int fa) {
    for(auto& v : G[u]) {
        if(v == fa) continue;
        dfs(v, u);
        chkmax(mxd[u], mxd[v] + 1);
        memcpy(g, f[u], sizeof(g));
        for(int i = 0; i <= mxd[u] && i <= k; i++) {
            for(int j = 0; j <= mxd[v] && j <= k; j++) {
                if(i + j + 1 <= k) {
                    add(g[max(i, j + 1)], 1LL * f[u][i] * f[v][j] % mod);
                }
            }
        }
        for(int j = 0; j <= mxd[v]; j++) add(g[j + 1], f[v][j]);
        for(int i = 0; i <= mxd[u]; i++) f[u][i] = g[i];
    }
    if(is[u]) {
        f[u][0] = 1;
        for(int i = 1; i <= mxd[u] && i <= k; i++) add(f[u][i], f[u][i]);
    }
}

int main() {
    scanf("%d%d%d", &n, &m, &k);
    for(int i = 2; i <= n; i++) {
        int u, v;
        scanf("%d%d", &u, &v);
        G[u].push_back(v);
        G[v].push_back(u);
    }
    for(int i = 1; i <= m; i++) {
        int x; scanf("%d", &x);
        is[x] = true;
    }
    dfs(1, 0);
    int ans = 0;
    for(int i = 0; i <= mxd[1]; i++)
        add(ans, f[1][i]);
    printf("%d\n", ans);
    return 0;
}

/*
*/

 

转载于:https://www.cnblogs.com/CJLHY/p/10890357.html

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值