#include <bits/stdc++.h>
using namespace std;
using ll = long long;

const int MOD = 998244353;

void solve()
{
    int n, k;
    cin >> n >> k;
    vector<vector<int>> node(n + 1);
    for (int i = 1; i <= n - 1; ++i) {
        int u, v;
        cin >> u >> v;
        node[u].push_back(v);
        node[v].push_back(u);
    }
    int sq = 500;
    if (k <= sq) {
        vector<vector<ll>> dp(n + 1, vector<ll>(k + 2));
        vector<int> sz(n + 1);
        function<void(int, int)> dfs = [&](int x, int f) {
            dp[x][1] = 1;
            sz[x] = 1;
            for (auto e : node[x]) {
                if (e == f)
                    continue;
                vector<ll> ndp(k + 2);
                dfs(e, x);
                for (int i = 0; i <= min(k + 1, sz[e]); ++i) {
                    for (int j = 0; j <= min(k + 1, sz[x]); ++j) {
                        if (i + j > k + 1)
                            continue;
                        ndp[i + j] += dp[e][i] * dp[x][j] % MOD;
                    }
                }

                sz[x] += sz[e];
                for (int i = 0; i <= min(k + 1, sz[x]); ++i)
                    ndp[i] %= MOD;
                swap(dp[x], ndp);
            }
            dp[x][0] = (dp[x][k] + dp[x][k + 1]) % MOD;
            return;
        };
        dfs(1, 0);
        cout << dp[1][0] << '\n';
    }
    else {
        vector<vector<vector<ll>>> dp(n + 1, vector<vector<ll>>(n / k + 1, vector<ll>(2)));
        vector<int> sz(n + 1);
        function<void(int, int)> dfs = [&](int x, int f) {
            sz[x] = 1;
            dp[x][0][0] = 1;
            for (auto e : node[x]) {
                if (e == f)
                    continue;
                vector<vector<ll>> ndp(n / k + 1, vector<ll>(2));
                dfs(e, x);
                for (int i = 0; i <= sz[e] / k; ++i) {
                    for (int j = 0; j <= sz[x] / k; ++j) {
                        for (int b1 = 0; b1 < 2; ++b1) {
                            for (int b2 = 0; b2 < 2; ++b2) {
                                int rem = (sz[e] - i * k) % (k + 1) + (sz[x] - j * k) % (k + 1) + (b1 + b2) * (k + 1);
                                if (rem > k + 1)
                                    continue;
                                ndp[i + j][rem == k + 1] += dp[e][i][b1] * dp[x][j][b2] % MOD;
                            }
                        }
                    }
                }
                sz[x] += sz[e];
                for (int i = 0; i <= sz[e] / k + sz[x] / k; ++i) {
                    ndp[i][0] %= MOD;
                    ndp[i][1] %= MOD;
                }
                swap(dp[x], ndp);
            }
            for (int i = 0; i <= n / k; ++i) {
                (dp[x][i][0] += dp[x][i][1]) %= MOD;
                if ((sz[x] - k * i) % (k + 1) == k)
                    (dp[x][i + 1][0] += dp[x][i][0]) %= MOD;
            }
            return;
        };
        dfs(1, 0);
        ll ans = 0;
        for (int i = 0; i <= n / k; ++i) {
            if ((n - i * k) % (k + 1) == 0)
                (ans += dp[1][i][0]) %= MOD;
        }
        cout << ans << '\n';
    }
    return;
}

int main()
{
    ios::sync_with_stdio(0);
    cin.tie(0);
    cout.tie(0);
    int T = 1;
    cin >> T;
    while (T--)
        solve();
    return 0;
}