首页 > 代码库 > BZOJ4449 : [Neerc2015]Distance on Triangulation

BZOJ4449 : [Neerc2015]Distance on Triangulation

首先拓扑,每次取出度数为$2$的点,这样可以把所有三角形都找到。

那么建出对偶图,会发现是一棵树。

对这棵树进行点分治,每次取出重心,DFS求出所有在里面的点,然后从重心$3$个点分别做一次BFS。

对于每个询问,如果不经过重心这个区域,那么递归求解,否则用BFS的结果回答即可。

时间复杂度$O(n\log n)$。

 

#include<cstdio>#include<algorithm>using namespace std;const int N=50010,M=100010,U=M*20;int n,m,i,x,y,z,d[N],g[N],v[N<<2],nxt[N<<2],ed,vis[N],h,t,q[N],ans[M];struct Q{int x,y;}que[M];inline void read(int&a){char c;while(!(((c=getchar())>=‘0‘)&&(c<=‘9‘)));a=c-‘0‘;while(((c=getchar())>=‘0‘)&&(c<=‘9‘))(a*=10)+=c-‘0‘;}inline void add(int x,int y){d[x]++;v[++ed]=y;nxt[ed]=g[x];g[x]=ed;}namespace Graph{int n,m,i,id[N][3],g[N],v[N<<1],nxt[N<<1],ok[N<<1],ed;int all,f[N],son[N],now,pos[N],last[N],cur,dis[3][N];int G[N],V[U],NXT[U],ED,p[N*3],cnt;struct E{int x,y,p;E(){}E(int _x,int _y,int _p){x=_x,y=_y,p=_p;}}e[N*3];inline bool cmp(const E&a,const E&b){return a.x==b.x?a.y<b.y:a.x<b.x;}inline void add(int x,int y){  v[++ed]=y;ok[ed]=1;nxt[ed]=g[x];g[x]=ed;  v[++ed]=x;ok[ed]=1;nxt[ed]=g[y];g[y]=ed;}inline void ADD(int x,int y){V[++ED]=y;NXT[ED]=G[x];G[x]=ED;}inline void newedge(int x,int y,int z){  n++;  id[n][0]=x,id[n][1]=y,id[n][2]=z;  e[++m]=E(min(x,y),max(x,y),n);  e[++m]=E(min(x,z),max(x,z),n);  e[++m]=E(min(y,z),max(y,z),n);}void findroot(int x,int y){  son[x]=1;f[x]=0;  for(int i=g[x];i;i=nxt[i])if(ok[i]&&v[i]!=y){    findroot(v[i],x);    son[x]+=son[v[i]];    if(son[v[i]]>f[x])f[x]=son[v[i]];  }  if(all-son[x]>f[x])f[x]=all-son[x];  if(f[x]<f[now])now=x;}void dfs(int x,int y,int z){  for(int i=0;i<3;i++){    pos[id[x][i]]=z;    last[id[x][i]]=cur;    p[++cnt]=id[x][i];  }  for(int i=g[x];i;i=nxt[i])if(ok[i]&&v[i]!=y)dfs(v[i],x,z);}inline void bfs(int S,int*d){  int i,x,y,h,t;static int q[N];  for(i=1;i<=cnt;i++)d[p[i]]=U;  d[q[h=t=1]=S]=0;  while(h<=t)for(i=::g[x=q[h++]];i;i=::nxt[i]){    y=::v[i];    if(d[y]<U||last[y]<cur)continue;    d[q[++t]=y]=d[x]+1;  }}void solve(int x){  if(!G[x])return;  f[0]=all=son[x],findroot(x,now=0);  int i,j,A,B;  cur++;  for(cnt=0,i=g[now];i;i=nxt[i])if(ok[i])dfs(v[i],now,v[i]);  for(i=0;i<3;i++){    p[++cnt]=A=id[now][i];    pos[A]=now,last[A]=cur;  }  for(i=0;i<3;i++)bfs(id[now][i],dis[i]);  for(cnt=0,i=G[x];i;i=NXT[i])p[++cnt]=V[i];G[x]=0;  for(i=1;i<=cnt;i++){    A=que[p[i]].x,B=que[p[i]].y;    if(pos[A]==pos[B])      if(pos[A]==now)ans[p[i]]=1;      else ADD(pos[A],p[i]);    else for(j=0;j<3;j++)ans[p[i]]=min(ans[p[i]],dis[j][A]+dis[j][B]);  }  for(i=g[now];i;i=nxt[i])if(ok[i])ok[i^1]=0,solve(v[i]);}void Main(){  sort(e+1,e+m+1,cmp);  for(ed=i=1;i<m;i++)if(e[i].x==e[i+1].x&&e[i].y==e[i+1].y)add(e[i].p,e[i+1].p);  son[1]=n;solve(1);}}int main(){  read(n);  for(i=1;i<n;i++)add(i,i+1),add(i+1,i);  add(1,n),add(n,1);  for(i=1;i<=n-3;i++)read(x),read(y),add(x,y),add(y,x);  for(h=i=1;i<=n;i++)if(d[i]==2)q[++t]=i;  while(h<=t){    x=q[h++];    if(d[x]!=2)continue;    vis[x]=1,y=0;    for(i=g[x];i;i=nxt[i])if(!vis[v[i]]){      if(y)z=v[i];else y=v[i];      if((--d[v[i]])==2)q[++t]=v[i];    }    Graph::newedge(x,y,z);  }  read(m);  for(i=1;i<=m;i++){    read(x),read(y);    if(x!=y)que[i].x=x,que[i].y=y,Graph::ADD(1,i),ans[i]=U;  }  Graph::Main();  for(i=1;i<=m;i++)printf("%d\n",ans[i]);  return 0;}

  

BZOJ4449 : [Neerc2015]Distance on Triangulation