bzoj3331 BJOI2013压力（圆方树+树上差分）

题目链接

QWQ我真的好菜啊。。

首先，根据题目的意思，对于每条路径，就是要让他的路径上的割点的$ans+1$。

那我们不妨把圆方树建出来，然后直接用类似树上差分标记路径的方式来标记路径就ok啦

因为在圆方树上，两个点之间的路径上的所有圆点，都是割点

然后直接统计一下就好。

QWQ

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<queue>
#include<map>
#include<set>
#define mk makr_pair
#define ll long long
using namespace std;
inline int read()
{
  int x=0,f=1;char ch=getchar();
  while (!isdigit(ch)) {if (ch=='-') f=-1;ch=getchar();}
  while (isdigit(ch)) {x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
  return x*f;
}
const int maxn = 4e5+1e2;
const int maxm = 2*maxn;
int point[maxn],nxt[maxm],to[maxm];
int cnt,n,m;
int point1[maxn],nxt1[maxm],to1[maxm];
int cnt1;
int low[maxn],f[maxn][21];
int dfn[maxn],s[maxn];
int dis[maxn],ans[maxn];
int deep[maxn];
int top;
int num;
int tot;
void addedge1(int x,int y)
{
 nxt1[++cnt1]=point1[x];
 to1[cnt1]=y;
 point1[x]=cnt1;
} 
void addedge(int x,int y)
{
//cout<<x<<" "<<y<<endl;
 nxt[++cnt]=point[x];
 to[cnt]=y;
 point[x]=cnt;
}
void tarjan(int x,int fa)
{
 dfn[x]=low[x]=++tot;
 s[++top]=x;
 for (int i=point1[x];i;i=nxt1[i])
 {
  int p = to1[i];
  if (p==fa) continue;
  if (!dfn[p])
  {
   tarjan(p,x);
   low[x]=min(low[x],low[p]);
   if(low[p]>=dfn[x])
   {
    ++num;
    addedge(num,x);
    addedge(x,num);
    do{
     addedge(s[top],num);
     addedge(num,s[top]);
     top--;
    }while (s[top+1]!=p);
   }
  }
  else
    low[x]=min(low[x],dfn[p]);
 }
}
void dfs(int x,int fa,int dep)
{
 deep[x]=dep;
 for (int i=point[x];i;i=nxt[i])
 {
  int p = to[i];
  if (p==fa) continue;
  f[p][0]=x;
  dfs(p,x,dep+1);
 }
}
void init()
{
 for (int j=1;j<=20;j++)
   for (int i=1;i<=num;i++)
      f[i][j]=f[f[i][j-1]][j-1];
}
int go_up(int x,int d)
{
 for (int i=0;i<=20;i++)
 {
  if ((1<<i)&d) x=f[x][i];
 }
 return x;
}
int lca(int x,int y)
{
 if (deep[x]>deep[y]) x=go_up(x,deep[x]-deep[y]);
 else y=go_up(y,deep[y]-deep[x]);
 if (x==y) return x;
 for (int i=20;i>=0;i--)
 {
  if (f[x][i]!=f[y][i])
  {
   x=f[x][i];
   y=f[y][i];
  }
 }
 return f[x][0];
}
int q;
void dp(int x,int fa)
{
 for (int i=point[x];i;i=nxt[i])
 {
  int p = to[i];
     if (p==fa) continue;
     dp(p,x);
     ans[x]+=ans[p];
 }
}
int main()
{
  n=read(),m=read();q=read();
  num=n;
  for (int i=1;i<=m;i++)
  {
    int x=read(),y=read();
    addedge1(x,y);
    addedge1(y,x); 
  }
  for (int i=1;i<=n;i++)
  {
   if(!dfn[i]) tarjan(i,0);
  }
  dfs(1,0,1);
  init();
  for (int i=1;i<=q;i++)
  {
    int x=read(),y=read();
    ans[x]++;
    ans[y]++;
    ans[f[lca(x,y)][0]]--;
  ans[lca(x,y)]--; 
  }
  dp(1,0);
  for (int i=1;i<=n;i++)
    cout<<ans[i]<<"\n";
  return 0;
}