F.DFS dfs序+线段树

F.DFS
题意:给你一颗初始树,然后两种操作,一个是原来两个点没有边然后加边,另外一个是两个点有边,删除边(不包括初始树边)。问你每次操作后,有多少个合法的点。合法点的定义:以这点为dfs的起点建出来的树与初始树是否不一样,如果每次都是不一样的,那么这就是不合格的点,否则就是合格的点。
题解:
在这里插入图片描述
在这里插入图片描述
我只要求得x,y这个环,不包括x,y内的所有点,和与环内点相连的点,都是不合格的点,我用总数减去不合格的。建议去手推一下具体是哪些点是不合格的,这是比较容易能看出来的。

#include<bits/stdc++.h>
#define ls o*2
#define rs o*2+1
#define m (l+r)/2
using namespace std;
const int maxn = 2e5+10;
int tr[maxn*4],lz[maxn*4];
int f[maxn][20],d[maxn],l[maxn],r[maxn],cnt;
vector<int>G[maxn];
map<int,int>mp[maxn];
void push(int o,int l,int r)
{
    if(lz[o])tr[o] = r-l+1;
    else if(l==r) tr[o] = 0;
    else tr[o] = tr[ls]+tr[rs];
}
void up(int o,int l,int r,int ql,int qr,int p)
{
    if(ql>qr)return ;
    if(ql<=l&&qr>=r) {
        lz[o] += p ;
       push(o,l,r);
        return ;
    }
    if(ql<=m) up(ls,l,m,ql,qr,p);
    if(qr>m) up(rs,m+1,r,ql,qr,p);
   push(o,l,r);
}
void dfs(int u,int fa,int dep)
{
   f[u][0] = fa;
   d[u] = dep;
   l[u] = ++cnt;
   for(auto v:G[u])
   if(v!=fa) dfs(v,u,dep+1);
   r[u] = cnt;
}
int main()
{
    int u,v;
    int n,q;scanf("%d%d",&n,&q);
    for(int i=1;i<n;i++){
      scanf("%d%d",&u,&v);
      G[u].push_back(v);
      G[v].push_back(u);
    }
    dfs(1,0,1);
    for(int j=1;j<20;j++)
    for(int i=1;i<=n;i++)
    f[i][j] = f[f[i][j-1]][j-1];
    int p;
    for(int i=1;i<=q;i++)
    {
      scanf("%d%d",&u,&v);
      if(l[u]>l[v]) swap(u,v);
      if(!mp[u][v]) mp[u][v] = 1 ,p = 1;
      else mp[u][v] = 0 , p = -1;
      if(r[u]>=r[v]){
        int x = v;
        for(int j=19;j>=0;j--)
        if(d[f[x][j]]>d[u]) x = f[x][j];
        up(1,1,n,l[x],l[v]-1,p) , up(1,1,n,r[v]+1,r[x],p);
      }
      else up(1,1,n,1,l[u]-1,p) , up(1,1,n,r[v]+1,n,p) , up(1,1,n,r[u]+1,l[v]-1,p);
      printf("%d\n",n-tr[1]);
     }
}

注意:就算那个节点lz[o]的值为0,它还需要tr[o] = tr[ls] + tr[rs];

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
#include<bits/stdc++.h> using namespace std; const int mx=1e5+1; int n,Q,x,y,d[mx],fa[mx],siz[mx],ev[mx],a[mx],son[mx],dfn[mx],cnt,id[mx],top[mx],ans[mx]; struct edge{int c,w,id,u,v;}e[mx*2]; struct que{int u,v,x,y;}q[mx*2]; struct tree{int l,r,lzy1,lzy2;}t[mx*4]; vector<edge> v[mx]; vector<int> es[mx]; vector<int> qs[mx]; //以下树剖 void dfs1(int f,int u) { d[u]=d[f]+1,fa[u]=f,siz[u]=1; int len=v[u].size(); for(int i=0;i<len;i++) { edge next=v[u][i]; int nv=next.v; if(nv==f) continue; ev[next.id]=nv,a[nv]=next.w; dfs1(u,nv); siz[u]+=siz[nv]; if(siz[nv]>siz[son[u]]) son[u]=nv; } } void dfs2(int f,int u) { dfn[u]=++cnt,id[cnt]=u,top[u]=f; if(son[u]) dfs2(f,son[u]); int len=v[u].size(); for(int i=0;i<len;i++) { int nv=v[u][i].v; if(nv==fa[u] || nv==son[u]) continue; dfs2(nv,nv); } } //以上树剖 //以下线段树 void pushup1(int x){t[x].lzy1=t[x<<1].lzy1+t[x<<1|1].lzy1;} void pushup2(int x){t[x].lzy2=t[x<<1].lzy2+t[x<<1|1].lzy2;} void build(int x,int l,int r) { t[x].l=l,t[x].r=r; if(l==r) { t[x].lzy1=a[id[l]],t[x].lzy2=0; return; } int mid=(l+r)/2; build(x<<1,l,mid);build(x<<1|1,mid+1,r); pushup1(x); } void chang1(int x,int obx,int w) { if(t[x].l==t[x].r){t[x].lzy1=w;return;} int mid=(t[x].l+t[x].r)>>1; if(obx<=mid) chang1(x<<1,obx,w); else chang1(x<<1|1,obx,w); pushup1(x); } void chang2(int x,int obx,int w) { if(t[x].l==t[x].r){t[x].lzy2=w;return;} int mid=(t[x].l+t[x].r)>>1; if(obx<=mid) chang2(x<<1,obx,w); else chang2(x<<1|1,obx,w); pushup2(x); } int find1(int x,int l,int r) { if(l<=t[x].l && r>=t[x].r) return t[x].lzy1; int mid=(l+r)>>1,s=0; if(l<=mid) s+=find1(x<<1,l,r); if(r>mid) s+=find1(x<<1|1,l,r); return s; } int find2(int x,int l,int r) { if(l<=t[x].l && r>=t[x].r) return t[x].lzy2; int mid=(l+r)>>1,s=0; if(l<=mid) s+=find2(x<<1,l,r); if(r>mid) s+=find2(x<<1|1,l,r); return s; } //以上线段树 int fans(int x,int y,int k) { int ans=0; while(top[x]!=top[y]) { if(d[top[x]]<d[top[y]]) swap(x,y); ans+=find1(1,dfn[top[x]],dfn[x]); ans+=find2(1,dfn[top[x]],dfn[x]); x=fa[top[x]]; } if(d[x]>d[y]) swap(x,y); if(x!=y) { ans+=find1(1,dfn[x]+1,dfn[y]); ans+=k*find2(1,dfn[x]+1,dfn[y]); } return ans; } int main() { cin >> n >> Q; for(int i=1;i<n;i++) { cin >> e[i].u >> e[i].v >> e[i].c >> e[i].w; e[i].id=i; v[e[i].u].push_back({e[i].u,e[i].v,e[i].c,e[i].w,e[i].id}); v[e[i].v].push_back({e[i].v,e[i].u,e[i].c,e[i].w,e[i].id}); es[e[i].c].push_back(i); } for(int i=1;i<=Q;i++) { cin >> q[i].x >> q[i].y >> q[i].u >> q[i].v; qs[q[i].x].push_back(i); } dfs1(1,1);dfs2(1,1);build(1,1,n); for(int i=1;i<n;i++) { int len=es[i].size(); for(int j=0;j<len;j++) { int k=ev[es[i][j]]; find1(1,dfn[k],0); find2(1,dfn[k],1); } for(int j=0;j<len;j++) { int k=qs[i][j]; ans[k]=fans(q[k].u,q[k].v,q[k].y); } for(int j=0;j<len;j++) { int k=ev[es[i][j]]; find1(1,dfn[k],e[es[i][j]].w); find2(1,dfn[k],0); } } for(int i=1;i<=Q;i++) cout<<ans[i]<<"\n"; return 0; }
最新发布
07-14

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值