题意:
在一颗树上,每个点有一个权值,求任意两点间路径上所有点的权值的gcd值,以及每个值出现的次数。
分析:
由于每个权值的范围在20000以内,所以能成为多个数的gcd的值很少,于是就可以暴力卡一卡就能过了。
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<map>
#include<vector>
#define SF scanf
#define PF printf
#define MAXN 200010
using namespace std;
typedef long long ll;
map<int,ll> c[MAXN];
int val[MAXN];
vector<int> a[MAXN];
ll cnt[MAXN];
int gcd(int x,int y){
if(y==0)
return x;
return gcd(y,x%y);
}
void dfs(int x,int fa){
c[x][val[x]]++;
for(int i=0;i<a[x].size();i++){
if(a[x][i]==fa)
continue;
int u=a[x][i];
dfs(u,x);
for(map<int,ll>::iterator it1=c[x].begin();it1!=c[x].end();it1++)
for(map<int,ll>::iterator it2=c[u].begin();it2!=c[u].end();it2++){
cnt[gcd(it1->first,it2->first)]+=it1->second*it2->second;
}
for(map<int,ll>::iterator it1=c[u].begin();it1!=c[u].end();it1++)
c[x][gcd(val[x],it1->first)]+=it1->second;
c[u].clear();
}
}
int n,u,v;
int main(){
SF("%d",&n);
for(int i=1;i<=n;i++){
SF("%d",&val[i]);
cnt[val[i]]++;
}
for(int i=1;i<n;i++){
SF("%d%d",&u,&v);
a[u].push_back(v);
a[v].push_back(u);
}
dfs(1,0);
for(int i=1;i<=200000;i++)
if(cnt[i])
PF("%d %I64d\n",i,cnt[i]);
}