HDU4812 D Tree（树的点分治）
题目大概说给一棵有点权的树，输出字典序最小的点对，使这两点间路径上点权的乘积模1000003的结果为k。
树的点分治搞了。因为是点权过根的两条路径的LCA会被重复统计，而注意到1000003是质数，所以这个用乘法逆元搞一下就OK了。还有要注意“治”的各个实现，把时间复杂度“控制”在O(nlogn)。
WA了几次，WA在漏了点到子树根的路径，还有每次分治忘了清空数组。
1 #include<cstdio> 2 #include<cstring> 3 #include<algorithm> 4 using namespace std; 5 #define INF (1<<30) 6 #define MAXN 111111 7 struct Edge{ 8 int v,next; 9 }edge[MAXN<<1]; 10 int NE,head[MAXN]; 11 void addEdge(int u,int v){ 12 edge[NE].v=v; edge[NE].next=head[u]; head[u]=NE++; 13 } 14 bool vis[MAXN]; 15 int mini,cen,size[MAXN]; 16 void getSize(int u,int fa){ 17 size[u]=1; 18 for(int i=head[u]; i!=-1; i=edge[i].next){ 19 int v=edge[i].v; 20 if(v==fa || vis[v]) continue; 21 getSize(v,u); 22 size[u]+=size[v]; 23 } 24 } 25 void getCen(int u,int fa,int &tot){ 26 int res=tot-size[u]; 27 for(int i=head[u]; i!=-1; i=edge[i].next){ 28 int v=edge[i].v; 29 if(v==fa || vis[v]) continue; 30 getCen(v,u,tot); 31 res=max(res,size[v]); 32 } 33 if(res<mini) mini=res,cen=u; 34 } 35 int getCen(int u){ 36 getSize(u,u); 37 mini=INF; 38 getCen(u,u,size[u]); 39 return cen; 40 } 41 long long ine(long long a){ 42 long long res=1,n=1000001; 43 while(n){ 44 if(n&1) res*=a,res%=1000003; 45 a*=a; a%=1000003; 46 n>>=1; 47 } 48 return res; 49 } 50 int n,k,val[MAXN]; 51 int ansx,ansy; 52 int record[1000003],tn,tmpx[MAXN],tmpy[MAXN],all[MAXN],an; 53 void dfs(int u,int fa,long long dist,int &top){ 54 int v=record[ine(dist)*k%1000003*top%1000003]; 55 if(v){ 56 if(u<v){ 57 if(u<ansx) ansx=u,ansy=v; 58 else if(u==ansx && v<ansy) ansy=u,ansy=v; 59 }else{ 60 if(v<ansx) ansx=v,ansy=u; 61 else if(v==ansx && u<ansy) ansy=v,ansy=u; 62 } 63 } 64 tmpx[tn]=u; tmpy[tn]=dist; ++tn; 65 all[an++]=dist; 66 for(int i=head[u]; i!=-1; i=edge[i].next){ 67 int v=edge[i].v; 68 if(v==fa || vis[v]) continue; 69 dfs(v,u,dist*val[v]%1000003,top); 70 } 71 } 72 void conquer(int u){ 73 an=0; 74 all[an++]=val[u]; 75 record[val[u]]=u; 76 for(int i=head[u]; i!=-1; i=edge[i].next){ 77 int v=edge[i].v; 78 if(vis[v]) continue; 79 tn=0; 80 dfs(v,v,(long long)val[u]*val[v]%1000003,val[u]); 81 for(int j=0; j<tn; ++j){ 82 if(record[tmpy[j]]==0 || record[tmpy[j]]>tmpx[j]) record[tmpy[j]]=tmpx[j]; 83 } 84 } 85 for(int i=0; i<an; ++i) record[all[i]]=0; 86 } 87 void divide(int u){ 88 u=getCen(u); 89 vis[u]=1; 90 conquer(u); 91 for(int i=head[u]; i!=-1; i=edge[i].next){ 92 int v=edge[i].v; 93 if(vis[v]) continue; 94 divide(v); 95 } 96 } 97 int main(){ 98 int a,b; 99 while(~scanf("%d%d",&n,&k)){ 100 for(int i=1; i<=n; ++i){ 101 scanf("%d",val+i); 102 } 103 NE=0; 104 memset(head,-1,sizeof(head)); 105 for(int i=1; i<n; ++i){ 106 scanf("%d%d",&a,&b); 107 addEdge(a,b); 108 addEdge(b,a); 109 } 110 memset(vis,0,sizeof(vis)); 111 ansx=ansy=INF; 112 divide(1); 113 if(ansx==INF) puts("No solution"); 114 else printf("%d %d\n",ansx,ansy); 115 } 116 return 0; 117 }
