【暴力】Codeforces990G GCD Counting

题意:

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

猜你喜欢

转载自blog.csdn.net/qq_34454069/article/details/80732484