#include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; typedef long long int LL; typedef pair P; typedef pair LP; const int INF=1<<30; const LL MAX=1e9+7; void array_show(int *array,int array_n,char middle=' '){ for(int i=0;i &vec_s,int vec_n=-1,char middle=' '){ if(vec_n==-1)vec_n=vec_s.size(); for(int i=0;i &vec_s,int vec_n=-1,char middle=' '){ if(vec_n==-1)vec_n=vec_s.size(); for(int i=0;i uft_q1; vector uft_parent; vector uft_num; public: void init(){ uft_parent.assign(uft_n,-1); uft_num.assign(uft_n,1); } union_find_tree(int uft_n_init){ assert(uft_n_init>=0); uft_n=uft_n_init; init(); } union_find_tree(){ uft_n=uft_N; init(); } int check_parent(int uft_x){ assert(uft_x>=0 && uft_x=0 && uft_x=0 && uft_yuft_num[uft_y])swap(uft_x,uft_y); uft_parent[uft_x]=uft_y; uft_num[uft_y]+=uft_num[uft_x]; return false; } int size(int pos){ pos=check_parent(pos); return uft_num[pos]; } }; long long int pow_mod(long long int p_a,long long int p_n,long long int p_p=1e9+7){ //p_a^p_n mod p_p long long int p_b=1,p_t=1; for(;p_b<=p_n;p_b<<=1); for(p_b>>=1;p_b>0;p_b>>=1){ p_t*=p_t; if(p_t>=p_p)p_t%=p_p; if(p_n&p_b)p_t*=p_a; if(p_t>=p_p)p_t%=p_p; } return p_t; } LL n,m,p; vector

path[220000]; LL sz[220000]; LL check(int pos,int bef){ LL s=0; LL a=0; for(auto node:path[pos]){ if(node.first==bef)continue; s+=check(node.first,pos); a+=sz[node.first]; s+=sz[node.first]*(n-sz[node.first])%MAX*pow_mod(p,node.second)%MAX; } sz[pos]=a+1; s%=MAX; return s; } int main(){ int i,j,k; int a,b,c; cin>>n>>m>>p; union_find_tree ua(n); for(i=0;i>a>>b>>c; a--,b--; if(ua.connect(a,b))continue; path[a].push_back(make_pair(b,c)); path[b].push_back(make_pair(a,c)); } cout<