#include #include using namespace std; namespace my{ using ml=atcoder::modint1000000007; auto&operator>>(istream&i,ml&x){int t;i>>t;x=t;return i;} auto&operator<<(ostream&o,const ml&x){return o<<(int)x.val();} #define eb emplace_back #define LL(...) ll __VA_ARGS__;lin(__VA_ARGS__) #define VV(n,...) vec__VA_ARGS__;setsize({n},__VA_ARGS__);vin(__VA_ARGS__) #define FO(n) for(ll ij=n;ij-->0;) #define FOR(i,...) for(auto[i,i##stop,i##step]=range(0,__VA_ARGS__);isync_with_stdio(0);cout<r{0,0,1};ll I=0;((r[I++]=a),...);if(!s&&I==1)swap(r[0],r[1]);r[0]-=s;if(s)r[2]*=-1;return r;} constexpr char newline=10; constexpr char space=32; constexpr ll width2(auto x){x|=1;ll r=0;while(x>0)x>>=1,++r;return r;} bool amax(auto&a,const auto&b){return a>auto&sort(auto&a,F f={}){ranges::sort(a,f);return a;} auto&unique(auto&a){sort(a).erase(ranges::unique(a).begin(),a.end());return a;} templatestruct queue:std::queue{ queue(const initializer_list&a={}){fe(a,e)this->emplace(e);} queue(const vector&a){fe(a,e)this->emplace(e);} ll size()const{return std::queue::size();} T pop(){T r=this->front();std::queue::pop();return r;} friend ostream&operator<<(ostream&o,queue q){while(q.size())o<0,space);return o;} }; templateauto pack_kth(const auto&...a){return get(make_tuple(a...));} templateauto pack_slice(const auto&...a){return[&](index_sequence){return array{get(forward_as_tuple(a...))...};}(make_index_sequence{});} templateconcept vectorial=is_base_of_v,V>; templatestruct vec_attr{using core_type=T;static constexpr int rank=0;}; templatestruct vec_attr{using core_type=typename vec_attr::core_type;static constexpr int rank=vec_attr::rank+1;}; templateusing core_t=vec_attr::core_type; templateistream&operator>>(istream&i,vector&v){fe(v,e)i>>e;return i;} templateostream&operator<<(ostream&o,const vector&v){fe(v,e)o<?newline:space);return o;} templatestruct vec; templatestruct tensor_helper{using type=vec::type>;}; templatestruct tensor_helper<0,T>{using type=T;}; templateusing tensor=typename tensor_helper::type; templatestruct vec:vector{ using vector::vector; vec(const vector&v){vector::operator=(v);} templaterequires(sizeof...(A)>=3)vec(A...a){const ll n=sizeof...(a)-1;auto t=pack_slice(a...);ll s[n];fo(i,n)s[i]=t[i];*this=make_vec(s,pack_kth(a...));} templatestatic auto make_vec(const ll(&s)[n],T x){if constexpr(i==n-1)return vec(s[i],x);else{auto X=make_vec(s,x);return vec(s[i],X);}} vec&operator^=(const vec&u){this->insert(this->end(),u.begin(),u.end());return*this;} vec operator^(const vec&u)const{return vec{*this}^=u;} vec&operator+=(const vec&u){vec&v=*this;assert(v.size()==u.size());fo(i,v.size())v[i]+=u[i];return v;} vec&operator-=(const vec&u){vec&v=*this;assert(v.size()==u.size());fo(i,v.size())v[i]-=u[i];return v;} vec operator+(const vec&u)const{return vec{*this}+=u;} vec operator-(const vec&u)const{return vec{*this}-=u;} vec&operator++(){fe(*this,e)++e;return*this;} vec&operator--(){fe(*this,e)--e;return*this;} vec operator-()const{vec v=*this;fe(v,e)e=-e;return v;} vec&operator%=(auto M){vec&v=*this;fe(v,e)e%=M;return v;} vec operator%(auto M)const{return vec{*this}%=M;} ll size()const{return vector::size();} auto pop_front(){auto r=*this->begin();this->erase(this->begin());return r;} template>auto sort(F f={})const{vec v=*this;ranges::sort(v,f);return v;} auto unique()const{auto r=sort();r.erase(ranges::unique(r).begin(),r.end());return r;} }; templaterequires(sizeof...(A)>=2)vec(A...a)->vec(declval>()))>>>; vec(ll)->vec; templatevoid setsize(const ll(&l)[n],A&...a){((a=vec::make_vec(l,core_t{})),...);} void lin(auto&...a){(cin>>...>>a);} void vin(auto&...a){fo(i,(a.size()&...))(cin>>...>>a[i]);} void dec(auto&...a){((--a),...);} templatestruct pow_linear:vec{pow_linear(ll a,ll n):vec(n+1,1){fo(i,n)(*this)[i+1]=(*this)[i]*a;}}; ll pow2ll(ll n){static const auto v=pow_linear(2,(ll)sizeof(ll)*CHAR_BIT-1);return v[n];} constexpr ll at2ll(auto x,ll i){return x>>i&1;} void sort(auto&...a){vec>v;(v.eb(a),...);sort(v);ll I=0;((a=v[I++]),...);} struct dsu{ ll n; vecpar_or_size; dsu(ll n):n(n),par_or_size(n,-1){} ll merge(ll a,ll b){ ll x=leader(a),y=leader(b); if(x==y)return x; if(-par_or_size[x]<-par_or_size[y])swap(x,y); par_or_size[x]+=par_or_size[y]; par_or_size[y]=x; return x; } ll leader(ll a){ assert(0<=a&&a>res(n); fo(u,n)res[leader(u)].eb(u); fe(res,e)ranges::sort(e); unique(res); if(res[0].empty())res.pop_front(); return res; } }; templatestruct edge{ ll from,to; WT wt; ll id; edge()=default; edge(ll from,ll to,WT wt=1,ll id=-1):from(from),to(to),wt(wt),id(id){} auto operator<=>(const edge&e)const{return wt<=>e.wt;} friend ostream&operator<<(ostream&o,const edge&e){return o<<"(to "<struct graph{ vec>>edges; graph()=default; graph(ll n):edges(n){} decltype(auto)operator[](ll i){return edges[i];} decltype(auto)operator[](ll i)const{return edges[i];} ll size()const{return edges.size();} friend ostream&operator<<(ostream&o,const graph&g){ fo(u,g.size()){ o<<"from "<&p){fo(i,p.size())if(p[i]!=-1)edges[p[i]].eb(p[i],i,1,i);} void add_edges(const vec&a,const vec&b){fo(i,a.size())edges[a[i]].eb(a[i],b[i],1,i);} void add_edges(const vec&a,const vec&b,const vec&w){fo(i,a.size())edges[a[i]].eb(a[i],b[i],w[i],i);} void add_edges(const vec>&es){fe(es,e)edges[e.from].eb(e);} }; templatestruct rooted_tree:graph{ vec>rev_edge; ll r; rooted_tree(ll n,ll r):graph(n),rev_edge(n,{-1,-1,WT(),-1}),r(r){} ll root()const{return r;} void add_edge_and_set_rev_edge(const edge&e){ this->edges[e.from].eb(e); this->rev_edge[e.to]=e; } }; templateauto sum_pair_and_tree_dist(const rooted_tree&g){ ll n=g.size(); WT wmax=0; fo(u,n)fe(g[u],e)amax(wmax,e.wt); T res=0; fo(b,width2(wmax)){ dsu uf(n); fo(u,n)fe(g[u],e)if(at2ll(e.wt,b))uf.merge(u,e.to); T t=0; fe(uf.groups(),v)t+=v.size()*(v.size()-1)/2; res+=t*pow2ll(b); } return res; } templateauto bfs_tree(const graph&g,ll s){ rooted_treetree(g.size(),s); vecused(g.size()); queue>q{{-1,s,WT{},-1}}; while(q.size()){ auto pu=q.pop(); ll u=pu.to; if(used[u])continue; used[u]=1; if(pu.from!=-1)tree.add_edge_and_set_rev_edge(pu); fe(g[u],e)if(!used[e.to])q.emplace(e); } return tree; } templateauto bfs_tree(const rooted_tree&g){return bfs_tree(g,g.root());} void pp(auto&&...a){ll n=sizeof...(a);((cout<0,space)),...);cout<g(N,0); g.add_edges(a,b,c); g.add_edges(b,a,c); g=bfs_tree(g); pp(sum_pair_and_tree_dist(g)); }}