#include #include #include namespace nachia{ template class CsrArray{ public: struct ListRange{ using iterator = typename std::vector::iterator; iterator begi, endi; iterator begin() const { return begi; } iterator end() const { return endi; } int size() const { return (int)std::distance(begi, endi); } Elem& operator[](int i) const { return begi[i]; } }; struct ConstListRange{ using iterator = typename std::vector::const_iterator; iterator begi, endi; iterator begin() const { return begi; } iterator end() const { return endi; } int size() const { return (int)std::distance(begi, endi); } const Elem& operator[](int i) const { return begi[i]; } }; private: int m_n; std::vector m_list; std::vector m_pos; public: CsrArray() : m_n(0), m_list(), m_pos() {} static CsrArray Construct(int n, const std::vector>& items){ CsrArray res; res.m_n = n; std::vector buf(n+1, 0); for(auto& [u,v] : items){ ++buf[u]; } for(int i=1; i<=n; i++) buf[i] += buf[i-1]; res.m_list.resize(buf[n]); for(int i=(int)items.size()-1; i>=0; i--){ res.m_list[--buf[items[i].first]] = items[i].second; } res.m_pos = std::move(buf); return res; } static CsrArray FromRaw(std::vector list, std::vector pos){ CsrArray res; res.m_n = pos.size() - 1; res.m_list = std::move(list); res.m_pos = std::move(pos); return res; } ListRange operator[](int u) { return ListRange{ m_list.begin() + m_pos[u], m_list.begin() + m_pos[u+1] }; } ConstListRange operator[](int u) const { return ConstListRange{ m_list.begin() + m_pos[u], m_list.begin() + m_pos[u+1] }; } int size() const { return m_n; } int fullSize() const { return (int)m_list.size(); } }; } // namespace nachia #include #include namespace nachia{ class StaticTopTree{ public: int m_n; nachia::CsrArray m_children; std::vector m_parent; std::vector m_parentEdge; int m_root; struct Node{ int p = -1, l = -1, r = -1; enum Type{ TyCompress, TyRake1, TyRake2, TyEdge } ty = TyEdge; }; std::vector m_node; StaticTopTree(int n, std::vector> tree, int root = 0){ m_n = n; assert((int)tree.size() == n-1); assert(0 <= root && root < n); m_root = root; nachia::CsrArray adj; { std::vector> dedges(n*2-2); for(int i=0; i::Construct(n, dedges); } m_parent.assign(n, -2); m_parentEdge.assign(n, -1); std::vector bfs = {m_root}; bfs.reserve(n); m_parent[m_root] = -1; for(int i=0; i<(int)bfs.size(); i++){ int p = bfs[i]; for(int nx : adj[p]) if(m_parent[nx] == -2){ m_parent[nx] = p; bfs.push_back(nx); } } for(int i=0; i::Construct(n, tree); std::vector nd(n, 1); for(int i=n-1; i>=1; i--) nd[m_parent[bfs[i]]] += nd[bfs[i]]; for(int p=0; p troot(n, -1); int trp = n*2-3; troot[bfs[0]] = --trp; for(int s : bfs) if(m_parent[s] < 0 || adj[m_parent[s]][0] != s){ struct SzNode { int sz, vid, nx; }; std::vector sznode; std::vector Hid = {0}; std::vector boundarySize; if(m_parent[s] >= 0){ sznode.push_back({ 1, m_parentEdge[s], adj[m_parent[s]][0] }); Hid.push_back(Hid.back()+1); } for(int p=s; ; p=adj[p][0]){ if(adj[p].size() == 0) break; for(int e=1; e Que = { { troot[s], 0, (int)sznode.size() } }; Que.reserve(sznode.size() * 2); for(int i=0; i<(int)Que.size(); i++){ int tp = Que[i].p, l = Que[i].l, r = Que[i].r; if(r-l == 1){ troot[sznode[l].nx] = tp; continue; } int m = Que[i].l; while(boundarySize[m] - boundarySize[l] < boundarySize[r] - boundarySize[m+1]) m++; if(Hid[l] == Hid[m]) m_node[tp].ty = Node::TyRake2; else if(Hid[m] == Hid[r]) m_node[tp].ty = Node::TyRake1; else m_node[tp].ty = Node::TyCompress; int pl = (m-l == 1) ? sznode[l].vid : -1; if(pl == -1) pl = --trp; int pr = (r-m == 1) ? sznode[m].vid : -1; if(pr == -1) pr = --trp; m_node[tp].l = pl; m_node[tp].r = pr; m_node[pl].p = tp; m_node[pr].p = tp; Que.push_back({ pl, l, m }); Que.push_back({ pr, m, r }); } } std::swap(m_children, adj); } StaticTopTree() : StaticTopTree(2, {{0,1}}) {} }; } // namespace nachia #include #include #include #include using Modint = atcoder::static_modint<998244353>; struct FPS{ std::vector A; int size() const { return A.size(); } }; FPS operator+(const FPS& l, const FPS& r){ FPS res; res.A.resize(std::max(l.size(), r.size())); for(int i=0; i,2>; FPS getSolid(const Val& v){ FPS res; res.A.resize(std::max({ v[0][0].size(), v[0][1].size()+1, v[1][0].size()+1, v[1][1].size()+2 })); for(int s : {0,1}) for(int t : {0,1}) for(int i=0; i> E(N-1); for(auto& e : E){ int u,v; scanf("%d%d", &u, &v); e = {u-1,v-1}; } auto tree = nachia::StaticTopTree(N, E); std::vector res(tree.m_node.size()); using NodeTy = nachia::StaticTopTree::Node::Type; for(int i=0; i<(int)tree.m_node.size(); i++){ auto& node = tree.m_node[i]; if(node.ty == NodeTy::TyEdge){ res[i][0][0] = FPS{ {1} }; res[i][0][1] = FPS{ {1} }; res[i][1][0] = FPS{ {1} }; res[i][1][1] = FPS{ {} }; } if(node.ty == NodeTy::TyCompress){ for(int s : {0,1}) for(int t : {0,1}){ res[i][s][t] = (res[node.l][s][0] * res[node.r][0][t]) + ((res[node.l][s][1] * res[node.r][1][t]) << 1); } } if(node.ty == NodeTy::TyRake1){ for(int s : {0,1}) for(int t : {0,1}){ res[i][s][t] = (res[node.l][s][t] * res[node.r][t][0]) + ((res[node.l][s][t] * res[node.r][t][1]) << 1); } } if(node.ty == NodeTy::TyRake2){ for(int s : {0,1}) for(int t : {0,1}){ res[i][s][t] = (res[node.l][s][0] * res[node.r][s][t]) + ((res[node.l][s][1] * res[node.r][s][t]) << 1); } } } FPS ans = getSolid(res.back()); ans.A.resize(N+1); for(int i=0; i<=N; i++){ if(i) printf(" "); printf("%u", ans.A[i].val()); } printf("\n"); return 0; }