Description
给出一棵 n n n个点的树,边有边权,对于 1 1 1~ n n n的这 n ! n! n!个排列,记排列 p 1 , . . . , p n p_1,...,p_n p1,...,pn这个排列的权值为在树上 p i p_i pi到 p i + 1 p_{i+1} pi+1的最短路长度之和,求所有排列的权值和
Input
第一行一整数 n n n表示点数,之后 n − 1 n-1 n−1行每行输入三个整数 u , v , w u,v,w u,v,w表示 u , v u,v u,v之间有一条权值为 w w w的树边
( 1 ≤ n ≤ 1 0 5 , 1 ≤ w ≤ 1 0 9 ) (1\le n\le 10^5,1\le w\le 10^9) (1≤n≤105,1≤w≤109)
Output
输出所有排列的权值和,结果模 1 0 9 + 7 10^9+7 109+7
Sample Input
3
1 2 1
2 3 1
3
1 2 1
1 3 2
Sample Output
16
24
Solution
考虑树上每条边对答案的贡献,对于边 u , v , w u,v,w u,v,w,假设 u u u的深度小于 v v v的深度,假设 v v v的子树中有 m m m个节点,那么从这 m m m个点中任取一点 x x x,从剩余 n − m n-m n−m个点中任取一点 y y y,若 x , y x,y x,y在排列中相邻,那么 x , y x,y x,y之间的最短路径必然经过 u v uv uv边,故 u v uv uv边对答案的贡献为 2 ⋅ m ⋅ ( n − m ) ⋅ w ⋅ ( n − 1 ) ! 2\cdot m\cdot (n-m)\cdot w\cdot (n-1)! 2⋅m⋅(n−m)⋅w⋅(n−1)!,一遍树形 D P DP DP求出每个点为根的子树中节点个数即可
Code
#include<cstdio>
#include<vector>
#include<algorithm>
using namespace std;
typedef long long ll;
typedef pair<int,int>P;
#define maxn 100005
#define mod 1000000007
int add(int x,int y)
{
x+=y;
if(x>=mod)x-=mod;
return x;
}
int mul(int x,int y)
{
ll z=1ll*x*y;
return z-z/mod*mod;
}
int T,n,Size[maxn],ans,fact[maxn];
vector<P>g[maxn];
void dfs(int u,int fa)
{
Size[u]=1;
for(int i=0;i<g[u].size();i++)
{
int v=g[u][i].first,w=g[u][i].second;
if(v==fa)continue;
dfs(v,u);
ans=add(ans,mul(w,mul(Size[v],n-Size[v])));
Size[u]+=Size[v];
}
}
int main()
{
fact[0]=1;
for(int i=1;i<=100000;i++)fact[i]=mul(fact[i-1],i);
while(~scanf("%d",&n))
{
for(int i=1;i<=n;i++)g[i].clear();
for(int i=1;i<n;i++)
{
int u,v,w;
scanf("%d%d%d",&u,&v,&w);
g[u].push_back(P(v,w));
g[v].push_back(P(u,w));
}
ans=0;
dfs(1,0);
ans=mul(add(ans,ans),fact[n-1]);
printf("%d\n",ans);
}
return 0;
}