2015年辽宁省赛Interesting Tree

题目描述

Recently, Miss Huang want to receive a Tree as her birthday gift! (What a interesting person!) As a interesting person, Miss Huang likes interesting tree.

The interesting tree is a tree with the Max interesting degree. The interesting degree equals LCM(H(1),H(2), .. , H(N)), H(x) means the height of node x in the tree and LCM means least common multiple.

Now Mr Chen has some non-rooted tree, he wants you to choose a root of the non-rooted tree and make the strange degree of tree be maximum;

The answer may be very large please print the answer mod 10^9+7.

输入

Your program will be tested on one or more test cases. In each test case:

First line a number N(1<=N<=10^5) represent the number of node.

The next N-1 line has two number u, v indicates that there is a edge between u, v.

输出

For every test case, print the following line:

answer

where answer is the maximum LCM.(mod 10^9+7)

样例输入

41 21 31 421 2

样例输出

题解:问题分两步:

1.求树的直径:即树种距离最长的两个点之间的距离

2.求从1到n这n个数的最小公倍数.

先说第一个问题,有两种方法:深搜和广搜.

先说深搜:对于树的某个节点,它的最深的两个儿子m1,m2.

那么很显然m1+m2+1有可能是树的直径.用这个数更新树的直径

再说广搜:1.随机选取一个节点A进行广搜,从而找到距离它最远的结点B.当然B可能不唯一.

2.B必然是此树某个直径的一个端点(树的直径可能不唯一).

从B开始广搜找到离B最远的C结点.BC极为直径.

再说第二个问题:有两种方法,其中一种是错的.

先说正确的:比如从1到9.既然有8,那么4,2就没法说话了,山上有老虎,猴子怎能称大王.

既然有9,那么3就没法说话了.

规律出来了:全体质数打表求出来,2,3,5,7.

再说错误的: 对于任意一个序列a[],求其全部元素的最小公倍数.怎么求?

求其全部元素的最大公约数怎么求?这个你会.两两求,从头扫描到尾.

但是这道题里不行,必须进行质因数分解,因为有取模运算!!!!!!!!!! 取模之后就没法再求最大公倍数了!!!!!!!!!!!!!!!!

比赛的时候,我就错在了这里,我真傻真的

#include<iostream>
#include<list>
#include<math.h>
#include<stdio.h>
#include<string.h>
using namespace std;
const int maxn = 1e5 + 7;
const int big = 1e9 + 7;
int N;
struct Edge{int to, next;
}e[maxn*2];
int g[maxn],ei;
bool vis[maxn];
int maxLen;
int prime[maxn / 5], psize;
void init(){bool is[maxn];psize = 0;memset(is, 1, sizeof(is));for (int i = 2; i < maxn; i++){if (is[i]){prime[psize++] = i;}for (int j = 0; i*prime[j] < maxn; j++){is[i*prime[j]] = 0;if (i%prime[j] == 0)break;}}
}
void push_back(int from, int to){ e[ei].next = g[from];e[ei].to = to;g[from] = ei++;
}
int deep(int x){int m1, m2;vis[x] = 1, m1 = m2 = 0;for (int i = g[x]; i; i = e[i].next){int t = e[i].to;if (vis[t] == false){int d = deep(t);if (m2<d){if (m1<d){ m2 = m1, m1 = d; }else m2 = d;}}}int len = m1 + m2 + 1;if (len > maxLen)maxLen = len;return m1 + 1;
}
int main(){//freopen("in.txt", "r", stdin);init();while (~scanf("%d", &N)){ ei = 1, memset(g, 0, sizeof(int)*(N + 1));int x, y; for (int i = 1; i < N; i++)scanf("%d%d", &x, &y), push_back(x, y), push_back(y, x); memset(vis, 0, sizeof(bool)*(N + 1));  maxLen = 0;deep(1);long long int ans = 1;for (int i = 0; i < psize; i++){int mi = floor(log(maxLen) / log(prime[i]));if (mi == 0)break;ans *= pow(prime[i], mi);ans %= big;}cout << ans << endl; }return 0;
}
/**************************************************************Problem: 1578User: 20124003Language: C++Result: 正确Time:181 msMemory:9568 kb
****************************************************************/

转载于:https://www.cnblogs.com/weiyinfu/p/5013884.html