Ostap & Tree
La solución que describe el editorial es pésima, pero podemos hacerlo mejor. Esto corre en . (¿demostración?)
#include <bits/stdc++.h>
using namespace std;
using ll = long long;
using ld = long double;
using db = double;
using str = string; // estilo python
using pi = pair<int, int>;
using pl = pair<ll, ll>;
using pd = pair<db, db>;
using vi = vector<int>;
using vb = vector<bool>;
using vl = vector<ll>;
using vd = vector<db>;
using vs = vector<str>;
using vpi = vector<pi>;
using vpl = vector<pl>;
using vpd = vector<pd>;
#define tcT template <class T
#define tcTU tcT, class U
// ^ esto deja todo con pinta rara, pero lo pruebo igual
tcT > using V = vector<T>;
tcT, size_t SZ > using AR = array<T, SZ>;
tcT > using PR = pair<T, T>;
// pares
#define mp make_pair
#define f first
#define s second
// vectores
// ojo: size(x), rbegin(x), rend(x) necesitan C++17
#define sz(x) int((x).size())
#define bg(x) begin(x)
#define all(x) bg(x), end(x)
#define rall(x) x.rbegin(), x.rend()
#define sor(x) sort(all(x))
#define rsz resize
#define ins insert
#define ft front()
#define bk back()
#define pb push_back
#define eb emplace_back
#define pf push_front
#define lb lower_bound
#define ub upper_bound
tcT > int lwb(V<T> &a, const T &b) { return int(lb(all(a), b) - bg(a)); }
// bucles
#define FOR(i, a, b) for (int i = (a); i < (b); ++i)
#define F0R(i, a) FOR(i, 0, a)
#define ROF(i, a, b) for (int i = (b) - 1; i >= (a); --i)
#define R0F(i, a) ROF(i, 0, a)
#define trav(a, x) for (auto &a : x)
const int MOD = 1e9 + 7; // 998244353;
const int MX = 2e5 + 5;
const ll INF = 1e18; // no demasiado cerca de LLONG_MAX
const ld PI = acos((ld)-1);
const int dx[4] = {1, 0, -1, 0}, dy[4] = {0, 1, 0, -1}; // para todo problema de grilla
mt19937 rng((uint32_t)chrono::steady_clock::now().time_since_epoch().count());
template <class T> using pqg = priority_queue<T, vector<T>, greater<T>>;
// operaciones de bits
// también ver https://gcc.gnu.org/onlinedocs/gcc/Other-Builtins.html
constexpr int pct(int x) { return __builtin_popcount(x); } // cantidad de bits en 1
constexpr int bits(int x) { // assert(x >= 0); // compatible con C++11 hasta que
// USACO actualice ...
return x == 0 ? 0 : 31 - __builtin_clz(x);
} // floor(log2(x))
constexpr int p2(int x) { return 1 << x; }
constexpr int msk2(int x) { return p2(x) - 1; }
ll cdiv(ll a, ll b) {
return a / b + ((a ^ b) > 0 && a % b);
} // divide a entre b redondeando hacia arriba
ll fdiv(ll a, ll b) {
return a / b - ((a ^ b) < 0 && a % b);
} // divide a entre b redondeando hacia abajo
tcT > bool ckmin(T &a, const T &b) { return b < a ? a = b, 1 : 0; } // a = min(a,b)
tcT > bool ckmax(T &a, const T &b) { return a < b ? a = b, 1 : 0; }
tcTU > T fstTrue(T lo, T hi, U f) {
hi++;
assert(lo <= hi); // asumiendo que f es creciente
while (lo < hi) { // primer índice tal que f es verdadero
T mid = lo + (hi - lo) / 2;
f(mid) ? hi = mid : lo = mid + 1;
}
return lo;
}
tcTU > T lstTrue(T lo, T hi, U f) {
lo--;
assert(lo <= hi); // asumiendo que f es decreciente
while (lo < hi) { // primer índice tal que f es verdadero
T mid = lo + (hi - lo + 1) / 2;
f(mid) ? lo = mid : hi = mid - 1;
}
return lo;
}
tcT > void remDup(vector<T> &v) { // ordenar y quitar duplicados
sort(all(v));
v.erase(unique(all(v)), end(v));
}
tcTU > void erase(T &t, const U &u) { // no borrar
auto it = t.find(u);
assert(it != end(t));
t.erase(it);
} // un elemento que no existe del (multi)conjunto
// ENTRADA
#define tcTUU tcT, class... U
tcT > void re(complex<T> &c);
tcTU > void re(pair<T, U> &p);
tcT > void re(V<T> &v);
tcT, size_t SZ > void re(AR<T, SZ> &a);
tcT > void re(T &x) { cin >> x; }
void re(db &d) {
str t;
re(t);
d = stod(t);
}
void re(ld &d) {
str t;
re(t);
d = stold(t);
}
tcTUU > void re(T &t, U &...u) {
re(t);
re(u...);
}
tcT > void re(complex<T> &c) {
T a, b;
re(a, b);
c = {a, b};
}
tcTU > void re(pair<T, U> &p) { re(p.f, p.s); }
tcT > void re(V<T> &x) { trav(a, x) re(a); }
tcT, size_t SZ > void re(AR<T, SZ> &x) { trav(a, x) re(a); }
tcT > void rv(int n, V<T> &x) {
x.rsz(n);
re(x);
}
// CONVERSIÓN A STRING
#define ts to_string
str ts(char c) { return str(1, c); }
str ts(const char *s) { return (str)s; }
str ts(str s) { return s; }
str ts(bool b) {
#ifdef LOCAL
return b ? "true" : "false";
#else
return ts((int)b);
#endif
}
tcT > str ts(complex<T> c) {
stringstream ss;
ss << c;
return ss.str();
}
str ts(V<bool> v) {
str res = "{";
F0R(i, sz(v)) res += char('0' + v[i]);
res += "}";
return res;
}
template <size_t SZ> str ts(bitset<SZ> b) {
str res = "";
F0R(i, SZ) res += char('0' + b[i]);
return res;
}
tcTU > str ts(pair<T, U> p);
tcT > str ts(T v) { // contenedores con begin(), end()
#ifdef LOCAL
bool fst = 1;
str res = "{";
for (const auto &x : v) {
if (!fst) res += ", ";
fst = 0;
res += ts(x);
}
res += "}";
return res;
#else
bool fst = 1;
str res = "";
for (const auto &x : v) {
if (!fst) res += " ";
fst = 0;
res += ts(x);
}
return res;
#endif
}
tcTU > str ts(pair<T, U> p) {
#ifdef LOCAL
return "(" + ts(p.f) + ", " + ts(p.s) + ")";
#else
return ts(p.f) + " " + ts(p.s);
#endif
}
// SALIDA
tcT > void pr(T x) { cout << ts(x); }
tcTUU > void pr(const T &t, const U &...u) {
pr(t);
pr(u...);
}
void ps() { pr("\n"); } // imprimir con espacios
tcTUU > void ps(const T &t, const U &...u) {
pr(t);
if (sizeof...(u)) pr(" ");
ps(u...);
}
// DEPURACIÓN
void DBG() { cerr << "]" << endl; }
tcTUU > void DBG(const T &t, const U &...u) {
cerr << ts(t);
if (sizeof...(u)) cerr << ", ";
DBG(u...);
}
#ifdef LOCAL // compilar con -DLOCAL, chk -> assert falso
#define dbg(...) \
cerr << "Line(" << __LINE__ << ") -> [" << #__VA_ARGS__ << "]: [", DBG(__VA_ARGS__)
#define chk(...) \
if (!(__VA_ARGS__)) \
cerr << "Line(" << __LINE__ << ") -> function(" << __FUNCTION__ \
<< ") -> CHK FAILED: (" << #__VA_ARGS__ << ")" << "\n", \
exit(0);
#else
#define dbg(...) 0
#define chk(...) 0
#endif
void setPrec() { cout << fixed << setprecision(15); }
void unsyncIO() { cin.tie(0)->sync_with_stdio(0); }
// E/S POR ARCHIVOS
void setIn(str s) { freopen(s.c_str(), "r", stdin); }
void setOut(str s) { freopen(s.c_str(), "w", stdout); }
void setIO(str s = "") {
unsyncIO();
setPrec();
// cin.exceptions(cin.failbit);
// lanza excepción si se hace algo ilegal
// ej. intentar leer una letra en un int
if (sz(s)) setIn(s + ".in"), setOut(s + ".out"); // para USACO
}
/**
* Descripción: operaciones de aritmética modular
* Fuente:
* KACTL
* https://codeforces.com/blog/entry/63903
* https://codeforces.com/contest/1261/submission/65632855 (tourist)
* https://codeforces.com/contest/1264/submission/66344993 (ksun)
* también ver https://github.com/ecnerwala/cp-book/blob/master/src/modnum.hpp
* (ecnerwal) Verificación: https://open.kattis.com/problems/modulararithmetic
*/
template <int MOD, int RT> struct mint {
static const int mod = MOD;
static constexpr mint rt() { return RT; } // raíz primitiva para FFT
int v;
explicit operator int() const {
return v;
} // explicit -> no convertir a int en silencio
mint() { v = 0; }
mint(ll _v) {
v = int((-MOD < _v && _v < MOD) ? _v : _v % MOD);
if (v < 0) v += MOD;
}
friend bool operator==(const mint &a, const mint &b) { return a.v == b.v; }
friend bool operator!=(const mint &a, const mint &b) { return !(a == b); }
friend bool operator<(const mint &a, const mint &b) { return a.v < b.v; }
friend void re(mint &a) {
ll x;
re(x);
a = mint(x);
}
friend str ts(mint a) { return ts(a.v); }
mint &operator+=(const mint &m) {
if ((v += m.v) >= MOD) v -= MOD;
return *this;
}
mint &operator-=(const mint &m) {
if ((v -= m.v) < 0) v += MOD;
return *this;
}
mint &operator*=(const mint &m) {
v = (ll)v * m.v % MOD;
return *this;
}
mint &operator/=(const mint &m) { return (*this) *= inv(m); }
friend mint pow(mint a, ll p) {
mint ans = 1;
assert(p >= 0);
for (; p; p /= 2, a *= a)
if (p & 1) ans *= a;
return ans;
}
friend mint inv(const mint &a) {
assert(a.v != 0);
return pow(a, MOD - 2);
}
mint operator-() const { return mint(-v); }
mint &operator++() { return *this += 1; }
mint &operator--() { return *this -= 1; }
friend mint operator+(mint a, const mint &b) { return a += b; }
friend mint operator-(mint a, const mint &b) { return a -= b; }
friend mint operator*(mint a, const mint &b) { return a *= b; }
friend mint operator/(mint a, const mint &b) { return a /= b; }
};
typedef mint<MOD, 5> mi; // 5 es raíz primitiva de ambos módulos habituales
typedef vector<mi> vmi;
typedef pair<mi, mi> pmi;
typedef vector<pmi> vpmi;
vector<vmi> scmb; // combinaciones chicas
void genComb(int SZ) {
scmb.assign(SZ, vmi(SZ));
scmb[0][0] = 1;
FOR(i, 1, SZ)
F0R(j, i + 1) scmb[i][j] = scmb[i - 1][j] + (j ? scmb[i - 1][j - 1] : 0);
}
vmi yes[MX], no[MX];
vi adj[MX];
int n, k;
void dfs(int x, int y) {
yes[x] = no[x] = {1}; // negro, no negro
// dist. del vértice bueno más cercano
// o dist. del vértice malo más lejano
auto ad = [](vmi &a, int b, mi c) {
while (sz(a) <= b) a.pb(0);
a[b] += c;
};
trav(t, adj[x]) if (t != y) {
dfs(t, x);
yes[t].insert(begin(yes[t]), 0);
no[t].insert(begin(no[t]), 0);
if (sz(no[t]) > k + 1) no[t].pop_back();
vmi YES, NO;
F0R(i, sz(yes[x]))
F0R(j, sz(yes[t])) ad(YES, min(i, j), yes[x][i] * yes[t][j]);
F0R(i, sz(no[x]))
F0R(j, sz(no[t])) ad(NO, max(i, j), no[x][i] * no[t][j]);
auto yesNo = [&](vmi good, vmi bad) {
F0R(i, sz(good)) F0R(j, sz(bad)) {
if (i + j <= k) ad(YES, i, good[i] * bad[j]);
else ad(NO, j, good[i] * bad[j]);
}
};
yesNo(yes[x], no[t]);
yesNo(yes[t], no[x]);
swap(yes[x], YES);
swap(no[x], NO);
}
}
int main() {
setIO();
re(n, k);
F0R(i, n - 1) {
int a, b;
re(a, b);
adj[a].pb(b), adj[b].pb(a);
}
dfs(1, 0);
mi ans = 0;
trav(t, yes[1]) ans += t;
ps(ans);
}