当前位置: 代码迷 >> 综合 >> 【暴力】Codeforces990G GCD Counting
  详细解决方案

【暴力】Codeforces990G GCD Counting

热度:129   发布时间:2023-09-27 07:29:37.0

题意:

在一颗树上,每个点有一个权值,求任意两点间路径上所有点的权值的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]);
}