直徑
vector很喜歡研究樹,他覺得求的直徑方法很妙。
最近他給帶點權的樹定義了一種特殊的直徑,如果一條鏈上的點兩不互質那麼這直徑,鏈就是合法的,特殊直徑最長。
他很快的秒了這道題,但是調皮的gg把他的代碼刪掉了,現在他把這題交給你。
題面描述:
在樹上求最長直徑,滿足直徑上的點的權值兩兩互質。
這道題早上想了另種方法,
第一種:
找到一個根,然後左子樹只能往上走,右子樹只能往下走,每個節點記錄一個vector,記錄每個可以最長可以達到的公因數,以及可以達到該公因數的值,這樣的時間複雜度是O(n2),但是時間複雜度也要O(n2),而且代碼好像很複雜,於是沒有打。
第二種:
先歐拉篩求出2到100000的質數(),因爲每個數最多有logn的質數,所以枚舉這些質數,並找出來這棵樹上 哪些節點可以整除,然後除完了之後找最長的整數鏈,(這個就是樹的最長直徑)。
但是這次題意有錯。。第二種方案雖然是正解(原題意的正解(即這條鏈上所有數的公共gcd不爲1)嗎,但是本題是所有數兩兩不互質,比如說……(6,10,15)這組數據是過不了的。)
代碼
#include<bits/stdc++.h>
using namespace std;
const int maxn=100005;
int n,l,c[maxn],len,value[maxn],visit[maxn],V[maxn];
bool b[maxn];
map< int, vector<int> > factor;
struct node
{
node *nxt;
int where;
} a[maxn<<1],*first[maxn];
inline void makelist(int x,int y)
{
a[++l].where=y;
a[l].nxt=first[x];
first[x]=&a[l];
}
pair<int, int> bfs(int init,int v,int round)
{
c[1]=init;
visit[init]=1;
int pos=0,will=0;
int k=1,l=1;
for (;l<=k; l++)
{
int m=c[l];
if (visit[m]>will)
{
will=visit[m];
pos=m;
}
for (node *x=first[m]; x; x=x->nxt)
if (!(value[x->where]%v)&&!visit[x->where])
{
visit[x->where]=visit[m]+1;
c[++k]=x->where;
}
}
if (round==0)
{
for (int i=1; i<=k; i++)
visit[c[i]]=0;
return make_pair(pos,will);
}
}
int calc(int v)
{
vector<int> idx=factor[v];
int will=0;
for (int i=0; i<idx.size(); i++)
if (!visit[idx[i]])
will=max(will,bfs(bfs(idx[i],v,0).first,v,1).second);
for (int i=0; i<idx.size(); i++)
visit[idx[i]]=0;
return will;
}
int main()
{
freopen("diameter.in","r",stdin);
freopen("diameter.out","w",stdout);
len=0;
memset(b,false,sizeof(b));
for (int i=2; i<=100000; i++)
{
if (!b[i]) c[++len]=i;
for (int j=1; j<=len; j++)
if (c[j]*i>100000) break;
else
{
b[c[j]*i]=true;
if (!(i%c[j])) break;
}
}
scanf("%d",&n);
memset(first,0,sizeof(first));
l=0;
for (int i=1; i<n; i++)
{
int x,y;
scanf("%d%d",&x,&y);
makelist(x, y);
makelist(y, x);
}
factor.clear();
for (int i=1; i<=n; i++)
{
int x;
scanf("%d",&x);
value[i]=x;
for (int j=1; c[j]*c[j]<=x; ++j)
if (!(x%c[j]))
{
if (factor.find(c[j])==factor.end())
factor[c[j]].clear();
factor[c[j]].push_back(i);
for (;!(x%c[j]);)
x/=c[j];
}
if (x!=1)
{
if (factor.find(x)==factor.end())
factor[x].clear();
factor[x].push_back(i);
}
}
memset(visit,0,sizeof(visit));
memset(V,0,sizeof(V));
int ans=0;
for (map< int, vector<int> >::iterator itr=factor.begin(); itr!=factor.end(); itr++)
ans=max(ans,calc(itr->first));
printf("%d",ans);
}