結果

問題 No.2949 Product on Tree
ユーザー Tatsu_mrTatsu_mr
提出日時 2024-10-25 22:20:30
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 606 ms / 2,000 ms
コード長 3,786 bytes
コンパイル時間 5,933 ms
コンパイル使用メモリ 321,312 KB
実行使用メモリ 92,172 KB
最終ジャッジ日時 2024-10-25 22:21:03
合計ジャッジ時間 32,764 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 46
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#include <atcoder/all>
#define For(i, a, b) for(long long i = a; i < b; i++)
#define rep(i, n) For(i, 0, n)
#define rFor(i, a, b) for(long long i = a; i >= b; i--)
#define ALL(v) (v).begin(), (v).end()
#define rALL(v) (v).rbegin(), (v).rend()
using namespace std;

using lint = long long;
using ld = long double;
using mint = atcoder::modint998244353;

int INF = 2000000000;
lint LINF = 1000000000000000000;

template <class T>
struct Edge {
    int from, to;
    T cost;
    int idx;
    
    Edge() {}
    Edge(int to_) : to(to_) {}
    Edge(int to_, T cost_) : to(to_), cost(cost_) {}
    Edge(int from_, int to_, int idx_) : from(from_), to(to_), idx(idx_) {}
    Edge(int from_, int to_, T cost_, int idx_) : from(from_), to(to_), cost(cost_), idx(idx_) {}
};

template <class T> using Graph = vector<vector<Edge<T>>>;
using graph = Graph<long long>;
using edge = Edge<long long>;

#define add emplace_back

template <class T, class E, T (*merge)(T, T), T (*add_edge)(Edge<E>, T), T (*add_node)(int, T), T (*id)()>
struct RerootingDP {
    private:
    int n;
    vector<vector<Edge<E>>> g;
    vector<vector<T>> dp, l, r;
    
    T dfs1(int v, int pv) {
        T res = id();
        for (int i = 0; i < (int)g[v].size(); i++) {
            int nv = g[v][i].to;
            if (nv == pv) {
                continue;
            }
            dp[v][i] = dfs1(nv, v);
            res = merge(add_edge(g[v][i], dp[v][i]), res);
        }
        return add_node(v, res);
    }
    
    void dfs2(int v, int pv, T num) {
        int sz = g[v].size();
        for (int i = 0; i < sz; i++) {
            int nv = g[v][i].to;
            if (nv == pv) {
                dp[v][i] = add_node(pv, num);
            }
        }
        l[v][0] = r[v][sz] = id();
        for (int i = 0; i < sz; i++) {
            l[v][i + 1] = merge(add_edge(g[v][i], dp[v][i]), l[v][i]);
            r[v][sz - i - 1] = merge(add_edge(g[v][sz - i - 1], dp[v][sz - i - 1]), r[v][sz - i]);
        }
        for (int i = 0; i < sz; i++) {
            int nv = g[v][i].to;
            if (nv == pv) {
                continue;
            }
            T nnum = merge(l[v][i], r[v][i + 1]);
            dfs2(nv, v, nnum);
        }
    }
    
    public:
    RerootingDP(vector<vector<Edge<E>>> g_) : n(g_.size()), g(g_), dp(n), l(n), r(n) {
        for (int i = 0; i < n; i++) {
            dp[i].resize(g[i].size());
            l[i].resize(g[i].size() + 1);
            r[i].resize(g[i].size() + 1);
        }
        dfs1(0, -1);
        dfs2(0, -1, id());
    }
    
    T get(int v) {
        T res = id();
        for (int i = 0; i < (int)g[v].size(); i++) {
            res = merge(add_edge(g[v][i], dp[v][i]), res);
        }
        return add_node(v, res);
    }
    
    void dump() {
        for (int v = 0; v < n; v++) {
            for (int i = 0; i < (int)g[v].size(); i++) {
                cout << v + 1 << "->" << g[v][i].to + 1 << " : ";
                cout << dp[v][i];
                cout << endl;
            }
        }
    }
};

vector<mint> A(200010);

mint merge(mint a, mint b) {
    return a + b;
}

mint add_edge(Edge<int> e, mint a) {
    return a;
}

mint add_node(int v, mint a) {
    return A[v] * a + A[v];
}

mint id() {
    return mint(0);
}

int main() {
    int n;
    cin >> n;
    rep(i, n) {
        int x;
        cin >> x;
        A[i] = mint(x);
    }
    Graph<int> g(n);
    rep(i, n - 1) {
        int u, v;
        cin >> u >> v;
        u--;
        v--;
        g[u].add(v);
        g[v].add(u);
    }
    RerootingDP<mint, int, merge, add_edge, add_node, id> rd(g);
    mint ans = mint(0);
    rep(i, n) {
        ans += rd.get(i);
    }
    ans -= accumulate(ALL(A), mint(0));
    ans /= mint(2);
    cout << ans.val() << endl;
}
0