[LOJ#6159][美团 CodeM 初赛 Round A][暴力即正解]最长树链

关注

原创最新推荐文章于 2026-05-20 23:06:32 发布 983 阅读 AI 写同款· GEO 优化 ›

FDTD Solutions-入射波长模式经常在文献中看到仿真光源使用TM/TE模式，这个确定是十分必要的，决定了仿真结果。但是也不能局限于选择TM/TE，因为不同领域、不同软件之间的使用习惯和命名可能不同。并且在现代光波导器件中，几乎全部都是混合模式，或者称之为quasi-TE（准TE模式）和quasi-TM（准TM模式），多数情况下都简称为TE模式和TM模式。要知道到底是TE模式还是TM模式，需要查看其电场的主分量（即振幅... 阅读详情

枚举每一个质数，权值是这个质数的倍数的节点设为1，否则为0，那么最长长链的长度就是点权为1的点构成的最长的链。

当然不用枚举所以质数，因为每个数最多有log个不同的质因数，所以只要把所有点权质因数分解，枚举出现的质数就可以了。

枚举质数后dfs，复杂度是nlogn的，因为每个点有log个不同的质因数，也就是说最多枚举到log个质数要dfs这个节点，所以复杂度是对的

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <vector>
#include <map>
#define pb push_back

using namespace std;

const int N=100010;

int n,u,v,cnt,L,ans;
int p[N],vis[N],a[N],pnt;

map<int,vector<int> > M;
vector<int> e[N];

inline char nc(){
  static char buf[100000],*p1=buf,*p2=buf;
  return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}

inline void rea(int &x){
  char c=nc(); x=0;
  for(;c>'9'||c<'0';c=nc());for(;c>='0'&&c<='9';x=x*10+c-'0',c=nc());
}

inline void Pre(int n){
  for(int i=2;i<=n;i++){
    if(!p[i]) p[++pnt]=i;
    for(int j=1;j<=pnt&&1LL*i*p[j]<=n;j++){
      p[p[j]*i]=1;
      if(i%p[j]==0) break;
    }
  }
}

void dfs(int u,int f,int p){
  L++; vis[u]=p;
  for(int v : e[u])
    if(v!=f&&a[v]%p==0) dfs(v,u,p);
}

int main(){
  rea(n); Pre(100000);
  for(int i=1;i<n;i++)
    rea(u),rea(v),e[u].pb(v),e[v].pb(u);
  for(int i=1;i<=n;i++){
    int x; rea(x); a[i]=x;
    for(int j=1;j<=pnt&&1LL*p[j]*p[j]<=x;j++)
      if(x%p[j]==0){
    while(x%p[j]==0) x/=p[j];
    M[p[j]].pb(i);
    if(x==1) break;
      }
    if(x!=1) M[x].pb(i);
  }
  for(auto s : M){
    int p=s.first;
    vector<int> v=s.second;
    for(int u : v)
      if(vis[u]!=p){
    L=0; dfs(u,0,p);
    ans=max(ans,L);
      }
  }
  printf("%d\n",ans); return 0;
}