点分治学习笔记
这里就不讲点分治的原理以及一些注意事项了，只讲几道例题吧（懒
[T1 Luogu3806]【模板】点分治1
给定一棵有n个点的树
询问树上距离为k的点对是否存在。
一看标题就是点分治了嘛...
这是之前在机房打的
对每个路径长度进行标记即可，考虑到重复标记，容斥原理搞一下对于子树的答案贡献权值-1即可（就是更新的时候把新增贡献改为-1）
1 #include<cstdio> 2 #include<queue> 3 #include<iostream> 4 #include<cstring> 5 #define R register 6 using namespace std; 7 inline int read(){ 8 int ans=0,f=1;char chr=getchar(); 9 while(!isdigit(chr)){if(chr=='-') f=-1;chr=getchar();} 10 while(isdigit(chr)){ans=(ans<<3)+(ans<<1)+chr-48;chr=getchar();} 11 return ans*f; 12 }const int M=100005<<1; 13 int head[M],ver[M],nxt[M],val[M],tot,n,m,x,root,sz[M],mx[M],vis[M],ANS[M],p,a[M],dis[M],s; 14 inline void add(int x,int y,int z){ver[++tot]=y;val[tot]=z;nxt[tot]=head[x];head[x]=tot;} 15 void Find_Gf(int x,int fa){//找根 16 sz[x]=1;mx[x]=0; 17 for(R int i=head[x];i;i=nxt[i]) 18 if(ver[i]!=fa&&!vis[ver[i]]) 19 Find_Gf(ver[i],x),sz[x]+=sz[ver[i]],mx[x]=max(mx[x],sz[ver[i]]); 20 mx[x]=max(s-sz[x],mx[x]); 21 if(mx[x]<mx[root]) root=x; 22 } 23 void dfs(int x,int len,int fa){//找路径 24 dis[++p]=a[x]; 25 for(R int i=head[x];i;i=nxt[i]) 26 if(ver[i]!=fa&&!vis[ver[i]]) 27 a[ver[i]]=len+val[i],dfs(ver[i],a[ver[i]],x); 28 } 29 void Calc(int x,int len,int w){//标记路径 30 p=0,a[x]=len,dfs(x,len,0); 31 for(R int i=1;i<=p;++i) 32 for(R int j=1;j<=p;++j) 33 if(i!=j) ANS[dis[i]+dis[j]]+=w; 34 } 35 void Solve(int x){//点分治 36 Calc(x,0,1),vis[x]=1; 37 for(R int i=head[x];i;i=nxt[i]) 38 if(!vis[ver[i]]) 39 Calc(ver[i],val[i],-1),s=sz[x],root=0,mx[0]=n,Find_Gf(ver[i],x),Solve(root); 40 }int main(){ 41 n=read();m=read(); 42 for(R int i=1,a,b,z;i<n;++i) 43 a=read(),b=read(),z=read(),add(a,b,z),add(b,a,z); 44 s=n,mx[0]=n,root=0; 45 Find_Gf(1,0),Solve(root); 46 while(m--) puts(ANS[x=read()]?"AYE":"NAY"); 47 return 0; 48 }
[T2 Luogu2634] [国家集训队]聪聪可可
给定一棵树，求树上路径是3的倍数的条数
看到树上路径联想到点分
显然我们不可能开一个桶存边...观察发现只要记录下边权对3取模的结果就好了
并且对一个点来说，它的子树的贡献是这样的：
$ANS=cnt[1]*cnt[2]*2+{cnt[0]}^2$
考虑排除重复计算，在处理子树的子树时减掉即可，模板题
总数显然就是树上随便两个点并且可以重复所以就是${n}^2$
关于答案...除一个gcd就好了
1 #include<cstdio> 2 #include<cstring> 3 #include<iostream> 4 #include<algorithm> 5 #define int long long 6 using namespace std; 7 const int M=100005<<1; 8 int head[M],ver[M],nxt[M],tot,val[M],n,m,s,vis[M],sz[M],mx[M],root,ans,cnt[M]; 9 inline void add(int x,int y,int z){ver[++tot]=y;nxt[tot]=head[x];head[x]=tot;val[tot]=z;} 10 inline int read(){ 11 int ans=0,f=1;char chr=getchar(); 12 while(!isdigit(chr)){chr=getchar();} 13 while(isdigit(chr)){ans=(ans<<3)+(ans<<1)+chr-48;chr=getchar();} 14 return ans; 15 } 16 inline void cmax(int &x,int y){if(x<y)x=y;} 17 void find_gf(int x,int fa){ 18 mx[x]=sz[x]=1; 19 for(int i=head[x];i;i=nxt[i]){ 20 if(vis[ver[i]]||ver[i]==fa) continue; 21 find_gf(ver[i],x); 22 cmax(mx[x],sz[ver[i]]); 23 sz[x]+=sz[ver[i]]; 24 }cmax(mx[x],s-sz[x]); 25 if(mx[root]>mx[x]) root=x; 26 } 27 void dfs(int x,int fa,int dis){ 28 ++cnt[dis%3]; 29 for(int i=head[x];i;i=nxt[i]){ 30 if(ver[i]==fa||vis[ver[i]]) continue; 31 dfs(ver[i],x,(dis+val[i]%3)); 32 } 33 } 34 int Solve(int x,int w){ 35 cnt[0]=cnt[1]=cnt[2]=0; 36 dfs(x,0,w); 37 return cnt[1]*cnt[2]*2+cnt[0]*cnt[0]; 38 } 39 void Divide(int x){ 40 ans+=Solve(x,0); 41 vis[x]=1; 42 for(int i=head[x];i;i=nxt[i]){ 43 if(vis[ver[i]]) continue; 44 ans-=Solve(ver[i],val[i]); 45 root=0,s=sz[ver[i]];mx[ver[i]]=1; 46 find_gf(ver[i],0);Divide(root); 47 } 48 } 49 signed main(){ 50 n=read(); 51 for(int i=1,x,y,z;i<n;++i){ 52 x=read(),y=read(),z=read();z%=3; 53 add(x,y,z),add(y,x,z); 54 } 55 root=0,s=n;mx[root]=0x3f3f3f3f; 56 find_gf(1,0); 57 Divide(root); 58 printf("%lld/%lld",ans/__gcd(ans,n*n),n*n/__gcd(ans,n*n)); 59 return 0; 60 }
T3 CF161D Distance in Tree
求一棵树上边权值和为k的条数
显然淀粉质无疑了...
先考虑权值小于等于K的条数，减去小于等于K-1的条数即可
在计算贡献的时候我们显然不可能双重循环寻找答案，时间复杂度会退化成O(${n}^2$$logn$),这样不如暴力+LCA求解方便快捷
其实这里面是含有一个单调性的
在点分的dfs过程中记录点到root的距离（用一个数组存储）,然后sort这个数组，O（n）扫描即可
时间复杂度O(n${log}^2$n),在可以承受的范围内（还有一些人直接把其中一个logn当常数来看了非严格来讲O($nlogn$)也不算错，毕竟logn最大也只有20不到）
这题还有一个双倍经验（Luogu P4178 Tree）
1 #include<cstdio> 2 #include<cstring> 3 #include<iostream> 4 #include<algorithm> 5 #define int long long 6 using namespace std; 7 const int M=100005<<1; 8 int head[M],ver[M],nxt[M],tot,val[M],n,m,s,vis[M],sz[M],mx[M],root,ans,a[M]; 9 inline void add(int x,int y,int z){ver[++tot]=y;nxt[tot]=head[x];head[x]=tot;val[tot]=z;} 10 inline int read(){ 11 int ans=0,f=1;char chr=getchar(); 12 while(!isdigit(chr)){chr=getchar();} 13 while(isdigit(chr)){ans=(ans<<3)+(ans<<1)+chr-48;chr=getchar();} 14 return ans; 15 } 16 inline void cmax(int &x,int y){if(x<y)x=y;} 17 void find_gf(int x,int fa){ 18 mx[x]=sz[x]=1; 19 for(int i=head[x];i;i=nxt[i]){ 20 if(vis[ver[i]]||ver[i]==fa) continue; 21 find_gf(ver[i],x); 22 cmax(mx[x],sz[ver[i]]); 23 sz[x]+=sz[ver[i]]; 24 }cmax(mx[x],s-sz[x]); 25 if(mx[root]>mx[x]) root=x; 26 } 27 int p; 28 void dfs(int x,int fa,int dis){ 29 a[++p]=dis; 30 for(int i=head[x];i;i=nxt[i]){ 31 if(ver[i]==fa||vis[ver[i]]) continue; 32 dfs(ver[i],x,dis+1); 33 } 34 } 35 int Solve(int x,int w){ 36 p=0;dfs(x,0,w); 37 sort(a+1,a+p+1); 38 int l=1,r=p,ans=0; 39 while(l<r) if(a[l]+a[r]<=m)ans+=r-l,l++; else r--;//计算<=k的 40 l=1,r=p; 41 while(l<r) if(a[l]+a[r]<m)ans-=r-l,l++; else r--;//计算<=k-1的 42 return ans; 43 } 44 void Divide(int x){ 45 ans+=Solve(x,0); 46 vis[x]=1; 47 for(int i=head[x];i;i=nxt[i]){ 48 if(vis[ver[i]]) continue; 49 ans-=Solve(ver[i],val[i]); 50 root=0,s=sz[ver[i]];mx[ver[i]]=1; 51 find_gf(ver[i],0);Divide(root); 52 } 53 } 54 signed main(){ 55 n=read();m=read(); 56 for(int i=1,x,y,z;i<n;++i){ 57 x=read(),y=read(),z=1; 58 add(x,y,z),add(y,x,z); 59 } 60 root=0,s=n;mx[root]=0x3f3f3f3f; 61 find_gf(1,0); 62 Divide(root); 63 cout<<ans; 64 return 0; 65 }
[T4 LuoguP4149] [IOI2011]Race
求一棵树上权值为K且覆盖边数最少的路径
就是把上面的问题结合一下，然后观察到K不是很大...
注意一点当前边权大于K的时候及时退出...不然桶会炸
1 #include<cstdio> 2 #include<queue> 3 #include<iostream> 4 #include<cstring> 5 #include<algorithm> 6 using namespace std; 7 inline int read(){ 8 int ans=0,f=1;char chr=getchar(); 9 while(!isdigit(chr)){if(chr=='-') f=-1;chr=getchar();} 10 while(isdigit(chr)){ans=(ans<<3)+(ans<<1)+chr-48;chr=getchar();} 11 return ans*f; 12 }const int M=200005; 13 int n,kk,head[M<<1],ver[M<<1],nxt[M<<1],val[M<<1],tot,mx[M],vis[M],sz[M],root,s,Ans,dep[M],dis[M],p,b[1000100]; 14 inline void add(int x,int y,int z){ver[++tot]=y;val[tot]=z;nxt[tot]=head[x];head[x]=tot;} 15 void Find_GF(int x,int fa){ 16 sz[x]=1,mx[x]=0; 17 for(register int i=head[x];i;i=nxt[i]){ 18 if(vis[ver[i]]||ver[i]==fa) continue; 19 Find_GF(ver[i],x); 20 sz[x]+=sz[ver[i]],mx[x]=max(mx[x],sz[ver[i]]); 21 }mx[x]=max(mx[x],s-sz[x]); 22 if(mx[x]<mx[root]) root=x; 23 } 24 void DFS(int x,int fa,int d1,int d2){ 25 if(d1>kk) return; 26 dis[++p]=d1;dep[p]=d2; 27 for(register int i=head[x];i;i=nxt[i]){ 28 if(ver[i]==fa || vis[ver[i]]) continue; 29 DFS(ver[i],x,d1+val[i],d2+1); 30 } 31 } 32 inline void Calc(int x){ 33 b[0]=0,p=0; 34 for(register int i=head[x];i;i=nxt[i]){ 35 if(vis[ver[i]]) continue; 36 int lst=p; 37 DFS(ver[i],x,val[i],1); 38 for(register int j=lst+1;j<=p;++j) 39 Ans=min(Ans,b[kk-dis[j]]+dep[j]); 40 for(register int j=lst+1;j<=p;++j) 41 b[dis[j]]=min(b[dis[j]],dep[j]); 42 }for(int i=1;i<=p;i++) b[dis[i]]=1e9; 43 } 44 void Divide(int x){ 45 vis[x]=1,Calc(x); 46 for(register int i=head[x];i;i=nxt[i]){ 47 if(vis[ver[i]]) continue; 48 root=0,s=sz[ver[i]];root=0; 49 Find_GF(ver[i],x); 50 Divide(root); 51 } 52 } 53 int main(){Ans=0x3f3f3f3f; 54 n=read(),kk=read(); 55 for(register int x,y,z,i=1;i<n;++i){ 56 x=read(),y=read(),z=read();++x,++y; 57 add(x,y,z),add(y,x,z); 58 }mx[0]=n+1,root=0,s=n;memset(b,0x3f,sizeof(b)); 59 Find_GF(1,0); 60 Divide(root); 61 if(Ans>=n)puts("-1"); 62 else printf("%d",Ans); 63 return 0; 64 }
有一些难一点的覆盖问题以后可能会贴上来，但是博主tcl...有可能会咕
[UPD4.09]
然而并没有咕...不过是最后一次更新这篇啦
[T5 LuoguP4886]快递员
题目大意：
给定一棵n个节点的树和m个点对，在树上求一点使这个点到一个点对中每个点的距离之和的最大值的最小值...
见guai了...一眼看到根本没想到是点分治（然鹅打起来也不太像点分治...)看了一些大佬的题解才知道
我们随便找一个点比如重心，以其为根节点的儿子节点向外扩展（DFS），记录到每个节点的距离dis[]、它的祖先（根节点的某一个儿子）rt[]...
然后每个点对更新一下maxn值...并且记录那些点对可以取到maxn值...显然由于我们在模拟点分的过程...如果两个满足条件的点对却不连在同一个祖先（rt[]）上也是要舍掉的...
注意这里的DFS过程中不像其他点分治过程一样不能访问已经vis[]过的点，因为在遍历的过程中是要记录每一个节点在当前情况下的距离的！！
如果在某一个情况下访问到了已经访问过的节点只能说明所有情况都访问过了...可以直接输出退出了
1 #include<cstdio> 2 #include<queue> 3 #include<iostream> 4 #include<cstring> 5 #define R register 6 #define GO(i) for(R int i=head[x];i;i=nxt[i]) 7 using namespace std; 8 inline int read(){ 9 int ans=0,f=1;char chr=getchar(); 10 while(!isdigit(chr)){if(chr=='-') f=-1;chr=getchar();} 11 while(isdigit(chr)){ans=(ans<<3)+(ans<<1)+chr-48;chr=getchar();} 12 return ans*f; 13 }const int M=1e5+5; 14 int head[M<<1],ver[M<<1],nxt[M<<1],val[M<<1],tot,n,m,p[M],sz[M],mx[M],vis[M],s,root,px[M],py[M],rt[M],dis[M],Ans; 15 inline void add(int x,int y,int z){ver[++tot]=y;val[tot]=z;nxt[tot]=head[x];head[x]=tot;} 16 inline void cmax(int &x,int y){if(x<y) x=y;} 17 inline void cmin(int &x,int y){if(x>y) x=y;} 18 void Find_GF(int x,int fa){ 19 sz[x]=1,mx[x]=0; 20 GO(i){ 21 if(ver[i]==fa||vis[ver[i]]) continue; 22 Find_GF(ver[i],x); 23 sz[x]+=sz[ver[i]],cmax(mx[x],sz[ver[i]]); 24 }cmax(mx[x],s-sz[x]); 25 if(mx[root]>mx[x]) root=x; 26 } 27 void DFS(int x,int fa,int root){ 28 rt[x]=root; 29 GO(i){ 30 if(ver[i]==fa) continue; 31 dis[ver[i]]=dis[x]+val[i]; 32 DFS(ver[i],x,root); 33 } 34 } 35 void Divide(int x){ 36 if(vis[x]){cout<<Ans;exit(0);} 37 vis[x]=1,dis[x]=0; 38 GO(i)dis[ver[i]]=val[i],DFS(ver[i],x,ver[i]); 39 int maxn=0,siz=0,lst=0; 40 for(R int i=1;i<=m;++i) 41 if(dis[px[i]]+dis[py[i]]>maxn) 42 maxn=dis[px[i]]+dis[py[i]],p[siz=1]=i; 43 else if(dis[px[i]]+dis[py[i]]==maxn) p[++siz]=i; 44 cmin(Ans,maxn); 45 for(R int i=1;i<=siz;++i){ 46 if(rt[px[p[i]]]!=rt[py[p[i]]]){cout<<Ans;exit(0);} 47 else if(!lst) lst=rt[px[p[i]]];else if(lst!=rt[px[p[i]]]){cout<<Ans;exit(0);} 48 }root=0;s=sz[lst]; 49 Find_GF(lst,x);Divide(root); 50 } 51 int main(){Ans=0x7fffffff; 52 n=read(),m=read(); 53 for(R int i=1,x,y,z;i<n;++i){ 54 x=read(),y=read(),z=read(); 55 add(x,y,z),add(y,x,z); 56 }root=0,s=n,mx[0]=0x3f3f3f3f; 57 for(R int i=1;i<=m;++i) px[i]=read(),py[i]=read(); 58 Find_GF(1,0);Divide(root); 59 return 0; 60 }
