【BZOJ 4229】 4229: 选择 (线段树+树链剖分)

4229: 选择

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 67  Solved: 41

Description

现在，我想知道自己是否还有选择。

给定n个点m条边的无向图以及顺序发生的q个事件。

每个事件都属于下面两种之一：

1、删除某一条图上仍存在的边

2、询问是否存在两条边不相交的路径可以从点u出发到点v

Input

第一行三个整数n,m,q

接下来m行，每行两个整数u,v，表示u和v之间有一条边

接下来q行，每行一个大写字母o和2个整数u、v，依次表示按顺序发生的q个事件：

当o为’Z’时，表示删除一条u和v之间的边

当o为’P’时，表示询问是否存在两条边不相交的路径可以从点u出发到点v

Output

对于每组询问，如果存在，输出Yes,否则输出No

Sample Input

7 8 7
1 2
1 3
1 4
2 3
3 4
3 7
7 4
5 6
Z 1 4
P 1 3
P 2 4
Z 1 3
P 1 3
Z 6 5
P 5 6

Sample Output

Yes
Yes
No
No

HINT

对于全部数据，max(n,m,q)<=100000

Source

 

 

【分析】

　　我的做法跟BZOJ 1969大致一样【我改了改代码交的。。

　　这题没说删完的图联通，所以最后还要处理一下。。

　　然后每次询问x到y的关键边数量是否为0就好了。

　　【详细做法见BZOJ 1969

 

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
#define Maxn 100010
#define Maxm 100010

struct node
{
    int x,y,next;
    int bj,id;
}t[Maxn*2],tt[Maxm*4];
int first[Maxn],len;

void ins(int x,int y)
{
    // printf("%d -> %d\n",x,y);
    t[++len].x=x;t[len].y=y;
    t[len].next=first[x];first[x]=len;
}

struct nnode
{
    int l,r,lc,rc,sm;
}tr[Maxn*2];

void upd(int x)
{
    if(tr[x].l==tr[x].r||tr[x].sm!=0) return;
    int lc=tr[x].lc,rc=tr[x].rc;
    tr[lc].sm=tr[rc].sm=0;
}

int tot;
int build(int l,int r)
{
    int x=++tot;
    tr[x].l=l;tr[x].r=r;
    if(l!=r)
    {
        int mid=(l+r)>>1;
        tr[x].lc=build(l,mid);
        tr[x].rc=build(mid+1,r);
    }
    else tr[x].lc=tr[x].rc=0;
    tr[x].sm=r-l+1;
    return x;
}

void change(int x,int l,int r)
{
    if(tr[x].l==l&&tr[x].r==r)
    {
        tr[x].sm=0;
        return;
    }
    upd(x);
    int mid=(tr[x].l+tr[x].r)>>1,lc=tr[x].lc,rc=tr[x].rc;
    if(r<=mid) change(lc,l,r);
    else if(l>mid) change(rc,l,r);
    else
    {
        change(lc,l,mid);
        change(rc,mid+1,r);
    }
    tr[x].sm=tr[lc].sm+tr[rc].sm;
}

int query(int x,int l,int r)
{
    if(tr[x].sm==0) return 0;
    if(tr[x].l==l&&tr[x].r==r) return tr[x].sm;
    upd(x);
    int mid=(tr[x].l+tr[x].r)>>1,lc=tr[x].lc,rc=tr[x].rc;
    if(r<=mid) return query(lc,l,r);
    else if(l>mid) return query(rc,l,r);
    else return query(lc,l,mid)+query(rc,mid+1,r);
}

int tp[Maxn],sum[Maxn],son[Maxn],dfn[Maxn],dep[Maxn];
int ff[Maxn];
int cnt;
void dfs(int x,int f)
{
    son[x]=0;sum[x]=1;dep[x]=dep[f]+1;
    ff[x]=f;
    for(int i=first[x];i;i=t[i].next) if(t[i].y!=f)
    {
        int y=t[i].y;
        dfs(y,x);
        sum[x]+=sum[y];
        if(son[x]==0||sum[son[x]]<sum[y]) son[x]=y;
    }
}

void dfs2(int x,int tpp)
{
    dfn[x]=++cnt;tp[x]=tpp;
    if(son[x]) dfs2(son[x],tpp);
    for(int i=first[x];i;i=t[i].next) if(t[i].y!=ff[x]&&t[i].y!=son[x])
    {
        int y=t[i].y;
        dfs2(y,y);
    }
}

void fchange(int x,int y)
{
    while(tp[x]!=tp[y])
    {
        if(dep[tp[x]]<dep[tp[y]]) swap(x,y);
        change(1,dfn[tp[x]],dfn[x]);
        x=ff[tp[x]];
    }
    if(dep[x]<dep[y]) swap(x,y);
    if(x!=y) change(1,dfn[y]+1,dfn[x]);
}

int fquery(int x,int y)
{
    int ans=0;
    while(tp[x]!=tp[y])
    {
        if(dep[tp[x]]<dep[tp[y]]) swap(x,y);
        ans+=query(1,dfn[tp[x]],dfn[x]);
        x=ff[tp[x]];
    }
    if(dep[x]<dep[y]) swap(x,y);
    if(x!=y)
    {
        // if(dfn[y]+1>dfn[x]) while(1);
        ans+=query(1,dfn[y]+1,dfn[x]);
    }
    return ans;
}

int fa[Maxn];
int ffa(int x)
{
    if(fa[x]!=x) fa[x]=ffa(fa[x]);
    return fa[x];
}

bool cmp(node x,node y)
{
    if(x.x==y.x&&x.y==y.y) return x.bj<y.bj;
    return (x.x==y.x)?(x.y<y.y):(x.x<y.x);
}
bool cmp2(node x,node y) {return x.id<y.id;}

int ans[Maxn];

char s[10];

void solve()
{
    int n,m,q;
    scanf("%d%d%d",&n,&m,&q);
    int ll=0;
    for(int i=1;i<=m;i++)
    {
        int x,y;
        ll++;
        scanf("%d%d",&tt[ll].x,&tt[ll].y);
        if(tt[ll].x>tt[ll].y) swap(tt[ll].x,tt[ll].y);
        tt[ll].bj=1;//cha ru
    }
    int ct=0;
    while(q--)
    {
        scanf("%s",s);
        ll++;
        scanf("%d%d",&tt[ll].x,&tt[ll].y);
        if(tt[ll].x>tt[ll].y) swap(tt[ll].x,tt[ll].y);
        if(s[0]=='Z') tt[ll].bj=-1;//shan chu
        else tt[ll].bj=0;//xun wen
        tt[ll].id=++ct;
    }
    sort(tt+1,tt+1+ll,cmp);
    for(int i=1;i<=n;i++) fa[i]=i;
    ct++;
    int pp=0;
    for(int i=1;i<=ll;i++)
    {
        if(tt[i].bj==0) continue;
        if(tt[i].bj==1)
        {
            if(pp==0||tt[i].x!=tt[pp].x||tt[i].y!=tt[pp].y)
            {
                if(ffa(tt[i].x)==ffa(tt[i].y))
                {
                    tt[i].bj=-1;
                    tt[i].id=ct;
                }
                else
                {
                    ins(tt[i].x,tt[i].y);
                    ins(tt[i].y,tt[i].x);
                    fa[ffa(tt[i].x)]=tt[i].y;
                }
            }
        }
        pp=i;
    }
    sort(tt+1,tt+1+ll,cmp2);
    for(int i=ll;i>=1;i--) if(tt[i].bj==-1)
    {
        if(ffa(tt[i].x)!=ffa(tt[i].y))
        {
            ins(tt[i].x,tt[i].y);
            ins(tt[i].y,tt[i].x);
            fa[ffa(tt[i].x)]=tt[i].y;
            tt[i].bj=1;
        }
    }
    dep[0]=0;cnt=0;
    for(int i=1;i<=n;i++) if(ffa(i)==i)
    {
        dfs(i,0);
        dfs2(i,i);
    }
    build(1,n);
    
    ans[0]=0;
    for(int i=ll;i>=1;i--)
    {
        if(tt[i].bj==1) continue;
        if(tt[i].bj==-1)
        {
            fchange(tt[i].x,tt[i].y);
        }
        else
        {
            if(ffa(tt[i].x)!=ffa(tt[i].y)) ans[++ans[0]]=-1;
            else ans[++ans[0]]=fquery(tt[i].x,tt[i].y);
        }
    }
    for(int i=ans[0];i>=1;i--)
    {
        if(ans[i]==0) printf("Yes\n");
        else printf("No\n");
    }
}

int main()
{
    solve();
    return 0;
}

【是不是打的有点丑

 

2017-03-27 10:33:13