#include using namespace std; #define DEBUG(x) cerr<<#x<<": "< #define vl vector #define vii vector< vector > #define vll vector< vector > #define vs vector #define pii pair #define pis pair #define psi pair #define pll pair template pair operator+(const pair &s, const pair &t) { return pair(s.first + t.first, s.second + t.second); } template pair operator-(const pair &s, const pair &t) { return pair(s.first - t.first, s.second - t.second); } template ostream& operator<<(ostream& os, pair p) { os << "(" << p.first << ", " << p.second << ")"; return os; } #define X first #define Y second #define rep(i,n) for(int i=0;i<(int)(n);i++) #define rep1(i,n) for(int i=1;i<=(int)(n);i++) #define rrep(i,n) for(int i=(int)(n)-1;i>=0;i--) #define rrep1(i,n) for(int i=(int)(n);i>0;i--) #define REP(i,a,b) for(int i=a;i bool chmax(T &a, const T &b) { if (a bool chmin(T &a, const T &b) { if (a>b) { a = b; return 1; } return 0; } #define UNIQUE(v) v.erase(std::unique(v.begin(), v.end()), v.end()); const ll inf = 1000000001; const ll INF = (ll)1e18 + 1; const long double pi = 3.1415926535897932384626433832795028841971L; #define Sp(p) cout<= 0; i--) { rfact[i] = (i + 1) * rfact[i + 1] % M; } } //http://drken1215.hatenablog.com/entry/2018/06/08/210000 //n���傫��fact���v�Z�ł��Ȃ��Ƃ��̂ق��̌v�Z���@�ɂ‚��ď����Ă��� ll nCr(ll n, ll r, ll M = MOD) { if (r > n) return 0; assert(fact[2] == 2); ll ret = fact[n]; if (rfact[r] == 0) { rfact[r] = mod_inverse(fact[r], M); } ret = (ret*rfact[r]) % M; if (rfact[n - r] == 0) { rfact[n - r] = mod_inverse(fact[n - r], M); } ret = (ret*rfact[n - r]) % M; return ret; } ll nHr(ll n, ll r) { return nCr(n+r-1, r); } const int N = 312345; // const int N = 3; vii G(N); vll dp(N, vl(8)); vl sub(N); vl dfs(int now, int par) { vl res(8); sub[now] = 1; rep (i, 8) res[i] = 1; for (int v: G[now]) { if (v == par) continue; vl temp = dfs(v, now); // if (now == 0) { // DEBUG_VEC(temp); // DEBUG_VEC(res); // } vl res2(8); rep (i, 8) { rep (j, 8) { if ((i & j) != 0) continue; (res2[i | j] += res[i] * temp[j] % MOD) %= MOD; } } swap(res, res2); sub[now] += sub[v]; } dp[now] = res; (res[0] += mod_pow(2, sub[now] - 1)) %= MOD; return res; } signed main() { int n; cin >> n; rep (i, n - 1) { int u, v; cin >> u >> v; u--; v--; G[u].push_back(v); G[v].push_back(u); } dfs(0, -1); // DEBUG_MAT(dp); ll ans = 0; rep (i, n) { ll c = 1; if (i != 0) c = mod_pow(2, n - sub[i] - 1); ans += dp[i][7] * c % MOD; ans %= MOD; } ans %= MOD; cout << ans << endl; }