C. Ilya And The Tree(树形搜索+dp)

C. Ilya And The Tree
time limit per test
2 seconds
memory limit per test
256 megabytes
input
standard input
output
standard output

Ilya is very fond of graphs, especially trees. During his last trip to the forest Ilya found a very interesting tree rooted at vertex 1. There is an integer number written on each vertex of the tree; the number written on vertex i is equal to ai.

Ilya believes that the beauty of the vertex x is the greatest common divisor of all numbers written on the vertices on the path from the root to x, including this vertex itself. In addition, Ilya can change the number in one arbitrary vertex to 0 or leave all vertices unchanged. Now for each vertex Ilya wants to know the maximum possible beauty it can have.

For each vertex the answer must be considered independently.

The beauty of the root equals to number written on it.

Input

First line contains one integer number n — the number of vertices in tree (1 ≤ n ≤ 2·105).

Next line contains n integer numbers ai (1 ≤ i ≤ n1 ≤ ai ≤ 2·105).

Each of next n - 1 lines contains two integer numbers x and y (1 ≤ x, y ≤ nx ≠ y), which means that there is an edge (x, y) in the tree.

Output

Output n numbers separated by spaces, where i-th number equals to maximum possible beauty of vertex i.

Examples
input
2
6 2
1 2
output
6 6 
input
3
6 2 3
1 2
1 3
output
6 6 6 
input
1
10
output
10 
 
        
思路 : 对于树上的每一个节点,它的最大值或者等于把该节点变为0得到的结果,即从根节点到父亲节点没有变0得到的公分母;或者父亲节点里面所有可能的gcd和该点的最大公分母的最大值。
 
        

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cmath>
#include <queue>
#include <stack>
#include <cstring>
#include <map>
#include <climits>
#include <vector>
using namespace std;

int gcd(int a,int b)
{
    return b == 0 ? a : gcd(b,a%b);
}

vector <int > vc[1000000];
vector <int> dp_[1000006];//记录每个点所有可能的gcd
int a[1000006];
int dp[1000005];// 记录不变0的gcd
int mk[1000000];
void dfs(int cur)
{
    for(int i = 0; i < vc[cur].size(); i++)
    {
        int tmp = vc[cur][i];
        if(mk[tmp]) continue;
        dp[tmp] = gcd(a[tmp],dp[cur]);
        dp_[tmp].push_back(dp[cur]);
        for(int j = 0; j < dp_[cur].size(); j++)
        {
            int x = gcd(a[tmp],dp_[cur][j]);
            dp_[tmp].push_back(x);
        }
        sort(dp_[tmp].begin(),dp_[tmp].end());//去重,减少内存消耗。
        dp_[tmp].erase(unique(dp_[tmp].begin(),dp_[tmp].end()),dp_[tmp].end());
        mk[tmp] = 1;
        dfs(tmp);
        mk[tmp] = 0;
    }
}

int main()
{
   int n;
   cin >> n;
   for(int i = 1; i <= n; i++)
       scanf("%d",&a[i]);
   for(int i = 1; i < n; i++)
   {
       int x,y;
       scanf("%d%d",&x,&y);
       vc[x].push_back(y);
       vc[y].push_back(x);
   }
   dp[1] = DP[1] = a[1];
   dp_[1].push_back(0);
   dp_[1].push_back(a[1]);
   mk[1] = 1;
   dfs(1);
   for(int i = 1; i < n; i++)
   {
       printf("%d ",dp_[i].back());
   }
   printf("%d\n",dp_[n].back());
   return 0;
}


  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值