[SDOI2018]战略游戏【圆方树+链并】

题目描述:

n个点m条边的无向连通图,每次询问给出SS个点,问有多少不在SS中的点使得在图中删去它后在SS中存在两点u,vu,v不连通。
2n105,n1m2105,2S,S21052\le n\le10^5,n-1\le m\le2*10^5,2\le|S|,\sum|S|\le2*10^5

题目分析:

圆方树上S个点之间的圆点个数即为所求。
求链并当然可以建一波虚树
只需要将S个点按照dfs排序后,设dis[x]表示x到根的圆点个数,依次加入dis[a[i]]+dis[a[i+1]]-2*dis[LCA(a[i],a[i+1])],相当于每次加入a[i]到a[i+1]的路径(没有算lca,相当于加入边),每次加完后除了最外围的一条链,其他的边都被算了两次。这样只需要最后再加上dis[a[k]]+dis[a[1]]-2*dis[LCA(a[k],a[1])],就将路径上的每条边都算了两次,除以2后加上[LCA(a[k],a[1])为圆点],最后减去|S|即可。

Code:

#include <bits/stdc++.h>
#define maxn 400005
using namespace std;
int n,m,Q,sz;
vector<int>G[maxn],E[maxn];
inline void lineG(int x,int y){G[x].push_back(y),G[y].push_back(x);}
inline void lineE(int x,int y){E[x].push_back(y),E[y].push_back(x);}
namespace Graph{
    int dfn[maxn],low[maxn],tim,S[maxn];
    void tarjan(int u){
        dfn[u]=low[u]=++tim;
        S[++S[0]]=u;
        for(int i=G[u].size()-1,v;i>=0;i--)
            if(!dfn[v=G[u][i]]){
                tarjan(v),low[u]=min(low[u],low[v]);
                if(low[v]>=dfn[u]){
                    lineE(++sz,u);
                    do lineE(sz,S[S[0]]); while(S[S[0]--]!=v);
                }
            }
            else low[u]=min(low[u],dfn[v]);
    }
}
const int LOG = 18;
int dis[maxn],dep[maxn],dfn[maxn],tim,f[maxn][LOG+2],seq[maxn];
void dfs(int u,int ff){
    dis[u]=dis[f[u][0]=ff]+(u<=n);
    dep[u]=dep[ff]+1,dfn[u]=++tim;
    for(int i=E[u].size()-1,v;i>=0;i--) if((v=E[u][i])!=ff) dfs(v,u);
}
int LCA(int u,int v){
    if(dep[u]<dep[v]) swap(u,v);
    for(int i=0,d=dep[u]-dep[v];i<=LOG;i++) if(d&1<<i) u=f[u][i];
    if(u==v) return u;
    for(int i=LOG;i>=0;i--) if(f[u][i]!=f[v][i]) u=f[u][i],v=f[v][i];
    return f[u][0];
}
inline bool cmp(int x,int y){return dfn[x]<dfn[y];}
int main()
{
    int T,k,x,y,ans;
    scanf("%d",&T);
    while(T--){
        tim=Graph::tim=Graph::S[0]=0;
        scanf("%d%d",&n,&m),sz=n;
        for(int i=1;i<=n<<1;i++) Graph::dfn[i]=0,G[i].clear(),E[i].clear();
        for(int i=1;i<=m;i++) scanf("%d%d",&x,&y),lineG(x,y);
        Graph::tarjan(1);
        dfs(1,0);
        for(int j=1;j<=LOG;j++) for(int i=1;i<=sz;i++) f[i][j]=f[f[i][j-1]][j-1];
        scanf("%d",&Q);
        while(Q--){
            scanf("%d",&k),ans=0;
            for(int i=1;i<=k;i++) scanf("%d",&seq[i]);
            sort(seq+1,seq+1+k,cmp),seq[k+1]=seq[1];
            for(int i=1;i<=k;i++) ans+=dis[seq[i]]+dis[seq[i+1]]-2*dis[x=LCA(seq[i],seq[i+1])];
            printf("%d\n",ans/2+(x<=n)-k);
        }
    }
}
发布了363 篇原创文章 · 获赞 129 · 访问量 5万+
展开阅读全文

没有更多推荐了,返回首页

©️2019 CSDN 皮肤主题: 精致技术 设计师: CSDN官方博客

分享到微信朋友圈

×

扫一扫,手机浏览