zoukankan      html  css  js  c++  java
  • bzoj5329: [Sdoi2018]战略游戏

    大家都说是裸题。。。

    圆方树把图变成树,看虚树上有多少圆点即可

    #include<cstdio>
    #include<iostream>
    #include<cstring>
    #include<cstdlib>
    #include<algorithm>
    #include<cmath>
    #include<vector>
    using namespace std;
    const int _=1e2;
    const int maxn=1e5+_;
    const int maxm=2*1e5+_;
    int read()
    {
        int x=0,f=1;char ch=getchar();
        while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
        while(isdigit(ch)){x=x*10+ch-'0';ch=getchar();}
        return x*f;
    }
    void write(int x)
    {
        if(x>=10)write(x/10);
        putchar(x%10+'0');
    }
    struct node
    {
        int x,y,next;
    }a[2*2*maxm];int len,last[2*maxn];
    void ins(int x,int y)
    {
    //    printf("%d %d
    ",x,y);
        len++;
        a[len].x=x;a[len].y=y;
        a[len].next=last[x];last[x]=len;
    }
    
    //-----------------------------------------def----------------------------------------------------
    
    int vz,dfn[maxn],low[maxn];
    int top,sta[2*2*maxn];
    int cnt;vector<int>vec[maxn];
    void V_DCC(int x)
    {
        dfn[x]=low[x]=++vz;
        sta[++top]=x;
        for(int k=last[x];k;k=a[k].next)
        {
            int y=a[k].y;
            if(dfn[y]==0)
            {
                V_DCC(y);
                low[x]=min(low[x],low[y]);
                if(low[y]>=dfn[x])
                {
                    cnt++;vec[cnt].clear();
                    int i;
                    do
                    {
                        i=sta[top];top--;
                        vec[cnt].push_back(i);
                    }while(i!=y);
                    vec[cnt].push_back(x);
                }
            }
            else low[x]=min(low[x],dfn[y]);
        }
    }
    void maketree(int id)
    {
        len=1;memset(last,0,sizeof(last));
        for(int i=1;i<=cnt;i++)
        {
            id++;
            for(int j=0;j<vec[i].size();j++)
                ins(id,vec[i][j]),ins(vec[i][j],id);
        }
    }
    
    //-----------------------------------------------composition------------------------------------------
    
    int lim,num[2*maxn];
    int z,ys[2*maxn];
    int    fa[30][2*maxn],dep[2*maxn];
    void dfs(int x)
    {
        for(int i=1;(1<<i)<=dep[x];i++)fa[i][x]=fa[i-1][fa[i-1][x]];
        ys[x]=++z;
        for(int k=last[x];k;k=a[k].next)
        {
            int y=a[k].y;
            if(y!=fa[0][x])
            {
                fa[0][y]=x;
                dep[y]=dep[x]+1;
                num[y]=num[x]+(y<=lim);
                dfs(y);
            }
        }
    }
    int LCA(int x,int y)
    {
        if(dep[x]<dep[y])swap(x,y);
        for(int i=25;i>=0;i--)
            if(dep[x]-dep[y]>=(1<<i))x=fa[i][x];
        if(x==y)return x;
        for(int i=25;i>=0;i--)
            if(dep[x]>=(1<<i)&&fa[i][x]!=fa[i][y])x=fa[i][x],y=fa[i][y];
        return fa[0][x];
    }
    
    //-----------------------------------------------prepare--------------------------------------------
    
    int p[maxn],ans;
    bool cmp(int p1,int p2){return ys[p1]<ys[p2];}
    void VirtualTree(int m)
    {
        sort(p+1,p+m+1,cmp);
        top=1;sta[top]=p[1];
        for(int i=2;i<=m;i++)
        {
            if(LCA(p[i],sta[top])==sta[top])
                {sta[++top]=p[i];continue;}
            
            while(top>1&&LCA(p[i],sta[top-1])!=sta[top-1])
                ans+=num[sta[top]]-num[sta[top-1]],top--;
                
            if(top==1||LCA(p[i],sta[top])!=sta[top-1])
            {
                int lca=LCA(p[i],sta[top]);
                ans+=num[sta[top]]-num[lca],top--;
                sta[++top]=lca;
                sta[++top]=p[i];
            }
            else
            {
                ans+=num[sta[top]]-num[sta[top-1]],top--;
                sta[++top]=p[i];
            }
        }
        while(top>1)
            ans+=num[sta[top]]-num[sta[top-1]],top--;
        ans+=num[sta[top]]-num[fa[0][sta[top]]];
    }
    
    int main()
    {
        freopen("a.in","r",stdin);
        freopen("a.out","w",stdout);
        int T;
        T=read();
        while(T--)
        {
            int n,m,x,y;
            n=read(),m=read();
            len=1;memset(last,0,sizeof(last));
            for(int i=1;i<=m;i++)
            {
                x=read(),y=read();
                ins(x,y),ins(y,x);
            }
    //        puts("");
            vz=top=cnt=0;
            memset(dfn,0,sizeof(dfn));
            V_DCC(1);
            maketree(n);
            num[1]=1;lim=n;dfs(1);
    //        puts("");
            
            //....pre.....
            
            int Q;
            Q=read();
            while(Q--)
            {
                m=read();
                for(int i=1;i<=m;i++)p[i]=read();
                ans=0;
                VirtualTree(m);
                write(ans-m),putchar('
    ');
            }
    //        break;
        }
        return 0;
    }
  • 相关阅读:
    Android.mk高级写法
    Android Bitmap和Canvas学习笔记
    cocos2dx 3.x 开发环境搭建
    quick cocos2d-x 2.2.4 window环境调试
    Android应用开发相关下载资源(2014/12/14更新)
    quick-cocos2dx-2.2.4环境搭建
    U盘重装Windows系统
    Visual Studio 32位64位的问题和如何编译32位64位工程的问题
    使用UE4.16开发Oculus需要Oculus软件版本高于1.11
    UE4打包程序没有声音-需要安装UE4PrereqSetup_x64.exe
  • 原文地址:https://www.cnblogs.com/AKCqhzdy/p/10563122.html
Copyright © 2011-2022 走看看