2016弱校联盟十一专场10.2---Around the World(深搜+组合数、逆元)
题目链接
https://acm.bnu.edu.cn/v3/problem_show.php?pid=52305
problem description
In ICPCCamp, there are n cities and (n−1) (bidirectional) roads between cities. The i-th road is between the ai-th and bi-th cities. It is guaranteed that cities are connected. Recently, there are 2×ci −1 new roads built between the ai-th and bi-th cities. Bobo soon comes up with an idea to travel around the world! He plans to start in city 1 and returns to city 1 after traveling along every road exactly once. It is clear that Bobo has many plans to choose from. He would like to find out the number of different plans, modulo (109 + 7). Note that two plans A and B are considered different only if there exists an i where the i-th traveled road in plan A is different from the i-th road in plan B.
Input
The first line contains an integer n (2 ≤ n ≤ 105 ). The i-th of the following (n − 1) lines contains 3 integers ai , bi , ci (1 ≤ ai , bi ≤ n, ci ≥ 1, c1 + c2 + · · · + cn−1 ≤ 106 ).
Output
An integer denotes the number of plans modulo (109 + 7).
Examples
standard input
3
1 2 1
2 3 1
3
1 2 1
1 3 2
standard output
4
144
题意:输入n 表示由n个节点构成的树,然后输入n-1条边,n-1行每行输入ai bi c 表示ai点和bi点之间有2*c条边相连,现在求从1号点开始走完所有边且每条边只走一次的方案数?
思路:深搜,然后排列组合,举两个例子:

第一个图片中:定义sum=1,从1号点开始向下深搜,到达孩子2号点的时候sum*=2! 表示从1到2号点走遍边的方法数,然后继续向下深搜到4号点,sum*=2! 然后返回再到5号点,sum*=4! 再返回走到6号点,sum*=2! ,然后返回2号点,这时sum*=4!/2! 为什么呢? 因为从2号点开始走完它的直系(儿子)孩子节点时,它的儿子路径之间存在先后顺序,所以乘上路径数的阶乘,但要除去重复的部分比如从2号点到5号点时这两条路径的先后顺序在之前的4!中已经考虑了,接下来就是相同的过程了......注意的是下面的孩子节点算完了,再计算上层节点时就不需要考虑了,但是相邻两层的节点之间还是有影响的,这个图看不出来,我用下面一个图作解释。
第二个图片中:按照之前的算法为sum=4!*2!=48 其实正确的解是96 ,少乘了一个2,为什么呢? 因为上面一层的路径对下面一层的路径产生了影响,上一层有两条路径,下面一层只有一条路径,那么可以分析存在两种情况:1、 1->2->3->2->1->2->1 2、 1->2->1->2->3->2->1 怎么产生的? 其实是在走上一层的路径时,在哪一次走的这一层的路径,用挡板法,设上一层有t条路径,这一层有s条路径,那么得在t次把下面的s条路径走完,也就是C(t+s-1,t-1) 。
代码如下:
#include <iostream>
#include <stdio.h>
#include <algorithm>
#include <string.h>
typedef long long LL;
const LL maxn=1e5+;
const LL mod=1e9+;
using namespace std;
struct Tree
{
LL to,next,c;
} edge[maxn*];
LL Inv[*maxn];
LL tot,head[maxn];
LL jc[maxn*];
LL sum;
void Init()
{
Inv[]=;
Inv[] = ;
for(LL i=; i<*maxn; i++)
Inv[i] = (mod-mod/i)*Inv[mod%i]%mod;
for(LL i=; i<*maxn; i++)
Inv[i] = (Inv[i-]*Inv[i])%mod;
}
void init()
{
tot=;
memset(head,-,sizeof(head));
}
void add_edge(LL u,LL v,LL c)
{
edge[tot].to=v;
edge[tot].c=c;
edge[tot].next=head[u];
head[u]=tot++;
} void init2()
{
jc[]=;
for(LL i=; i<maxn*; i++)
jc[i]=((jc[i-]*i)%mod);
}
LL road[maxn];
void dfs(LL u,LL pre)
{
road[u]=;
for(LL i=head[u]; i!=-; i=edge[i].next)
{
LL v=edge[i].to;
if(v==pre) continue;
dfs(v,u);
road[u]+=edge[i].c;
sum=(((sum*jc[edge[i].c*])%mod)*Inv[edge[i].c])%mod;
sum=((sum* jc[(road[v]+edge[i].c-)]%mod)*Inv[edge[i].c-] )%mod;
sum=(sum*Inv[road[v]])%mod;
}
sum=(sum*jc[road[u]])%mod;
} int main()
{
Init();
init2();
LL n,u,v,c;
while(scanf("%lld",&n)!=-)
{
sum=;
init();
for(LL i=; i<n; i++)
{
scanf("%lld%lld%lld",&u,&v,&c);
add_edge(u,v,c);
add_edge(v,u,c);
}
dfs(,);
printf("%lld\n",sum);
}
return ;
}
2016弱校联盟十一专场10.2---Around the World(深搜+组合数、逆元)的更多相关文章
- 2016弱校联盟十一专场10.3---Similarity of Subtrees(深搜+hash、映射)
题目链接 https://acm.bnu.edu.cn/v3/problem_show.php?pid=52310 problem description Define the depth of a ...
- 2016弱校联盟十一专场10.5---As Easy As Possible(倍增)
题目链接 https://acm.bnu.edu.cn/v3/contest_show.php?cid=8506#problem/A problem description As we know, t ...
- (2016弱校联盟十一专场10.3) D Parentheses
题目链接 把左括号看成A右括号看成B,推一下就行了.好久之前写的,推到最后发现是一个有规律的序列. #include <bits/stdc++.h> using namespace std ...
- (2016弱校联盟十一专场10.3) B.Help the Princess!
题目链接 宽搜一下就行. #include <iostream> #include<cstdio> #include<cstring> #include<qu ...
- (2016弱校联盟十一专场10.3) A.Best Matched Pair
题目链接 #include<cstdio> #include<cstring> #include<algorithm> #include<stack> ...
- 2016弱校联盟十一专场10.3---We don't wanna work!(STL--set的使用)
题目链接 https://acm.bnu.edu.cn/v3/contest_show.php?cid=8504#problem/C 代码如下: #include <iostream> # ...
- (2016弱校联盟十一专场10.2) A.Nearest Neighbor Search
题目链接 水题,算一下就行. #include <bits/stdc++.h> using namespace std; typedef long long ll; ll x[],y[], ...
- (2016弱校联盟十一专场10.2) E.Coins
题目链接 很久之前写的了,好像是对拍打表过的,推一下就行了. #include <bits/stdc++.h> using namespace std; typedef long long ...
- (2016弱校联盟十一专场10.5) F. Fibonacci of Fibonacci
题目链接 题目大意就是这个,先找出下标的循环节,再快速幂对20160519取余就行了. 找出下标循环节: #include <cstdio> #include <iostream&g ...
随机推荐
- 每天一个linux命令(16):which命令
我们经常在linux要查找某个文件,但不知道放在哪里了,可以使用下面的一些命令来搜索: which 查看可执行文件的位置. whereis 查看文件的位置. ...
- eclipse打包apk
eclipse打包apk
- 经典String str = new String("abc")内存分配问题
出自:http://blog.csdn.net/ycwload/article/details/2650059 今天要找和存储管理相关的一些知识,网上搜了半天也没有找到完善的(30%的程度都不到),没 ...
- Liferay7 BPM门户开发之39: Form表单提交的ProcessAction处理
在v6.2开始后,需要设置<requires-namespaced-parameters>false</requires-namespaced-parameters> 来避免 ...
- react-native —— 在Windows下搭建React Native Android开发环境
在Windows下搭建React Native Android开发环境 前段时间在开发者头条收藏了 @天地之灵_邓鋆 分享的<在Windows下搭建React Native Android开发环 ...
- 让我们一起写出更有效的CSharp代码吧,少年们!
周末空闲,选读了一下一本很不错的C#语言使用的书,特此记载下便于对项目代码进行重构和优化时查看. Standing On Shoulders of Giants,附上思维导图,其中标记的颜色越深表示在 ...
- Oracle 11g系列:数据库
1.创建Oracle数据库 创建Oracle数据库的最常用工具为Database Configuration Assistant(数据库配置助手),依次选择[开始]|[所有程序]|[Oracle-Or ...
- js连续指定两次或者多次的click事件(解决办法)
setTimeout (表达式,延时时间)setTimeout(表达式,交互时间)延时时间/交互时间是以豪秒为单位的(1000ms=1s) setTimeout 在执行时,是在载入后延迟指定时间后, ...
- java.util.Scanner简单应用
import java.util.Scanner; import java.io.*; public class FileScannerTest{ public static void main(St ...
- Bounce.js – 快速创建漂亮的 CSS3 动画效果
Bounce.js 是一个用于制作漂亮的 CSS3 关键帧动画的 JavaScript 库,使用其特有的方式生成的动画效果.只需添加一个组件,选择预设,然后你就可以得到一个短网址或者导出为 CSS 代 ...