1. 程式人生 > >bzoj 4196 [Noi2015]軟體包管理器 (樹鏈剖分+線段樹)

bzoj 4196 [Noi2015]軟體包管理器 (樹鏈剖分+線段樹)

4196: [Noi2015]軟體包管理器

Time Limit: 10 Sec  Memory Limit: 512 MB
Submit: 2852  Solved: 1668
[Submit][Status][Discuss]

Description

 Linux使用者和OSX使用者一定對軟體包管理器不會陌生。通過軟體包管理器,你可以通過一行命令安裝某一個軟體包,然後軟體包管理器會幫助你從軟體源下載軟體包,同時自動解決所有的依賴(即下載安裝這個軟體包的安裝所依賴的其它軟體包),完成所有的配置。Debian/Ubuntu使用的apt-get,Fedora/CentOS使用的yum,以及OSX下可用的homebrew都是優秀的軟體包管理器。

你決定設計你自己的軟體包管理器。不可避免地,你要解決軟體包之間的依賴問題。如果軟體包A依賴軟體包B,那麼安裝軟體包A以前,必須先安裝軟體包B。同時,如果想要解除安裝軟體包B,則必須解除安裝軟體包A。現在你已經獲得了所有的軟體包之間的依賴關係。而且,由於你之前的工作,除0號軟體包以外,在你的管理器當中的軟體包都會依賴一個且僅一個軟體包,而0號軟體包不依賴任何一個軟體包。依賴關係不存在環(若有m(m≥2)個軟體包A1,A2,A3,…,Am,其中A1依賴A2,A2依賴A3,A3依賴A4,……,Am−1依賴Am,而Am依賴A1,則稱這m個軟體包的依賴關係構成環),當然也不會有一個軟體包依賴自己。 現在你要為你的軟體包管理器寫一個依賴解決程式。根據反饋,使用者希望在安裝和解除安裝某個軟體包時,快速地知道這個操作實際上會改變多少個軟體包的安裝狀態(即安裝操作會安裝多少個未安裝的軟體包,或解除安裝操作會解除安裝多少個已安裝的軟體包),你的任務就是實現這個部分。注意,安裝一個已安裝的軟體包,或解除安裝一個未安裝的軟體包,都不會改變任何軟體包的安裝狀態,即在此情況下,改變安裝狀態的軟體包數為0。  

Input

輸入檔案的第1行包含1個正整數n,表示軟體包的總數。軟體包從0開始編號。

隨後一行包含n−1個整數,相鄰整數之間用單個空格隔開,分別表示1,2,3,…,n−2,n−1號軟體包依賴的軟體包的編號。 接下來一行包含1個正整數q,表示詢問的總數。 之後q行,每行1個詢問。詢問分為兩種: installx:表示安裝軟體包x uninstallx:表示解除安裝軟體包x 你需要維護每個軟體包的安裝狀態,一開始所有的軟體包都處於未安裝狀態。對於每個操作,你需要輸出這步操作會改變多少個軟體包的安裝狀態,隨後應用這個操作(即改變你維護的安裝狀態)。  

Output

輸出檔案包括q行。

輸出檔案的第i行輸出1個整數,為第i步操作中改變安裝狀態的軟體包數。  

Sample Input

7
0 0 0 1 1 5
5
install 5
install 6
uninstall 1
install 4
uninstall 0

Sample Output

3
1
3
2
3

HINT

 

 一開始所有的軟體包都處於未安裝狀態。


安裝 5 號軟體包,需要安裝 0,1,5 三個軟體包。
之後安裝 6 號軟體包,只需要安裝 6 號軟體包。此時安裝了 0,1,5,6 四個軟體包。
解除安裝 1 號軟體包需要解除安裝 1,5,6 三個軟體包。此時只有 0 號軟體包還處於安裝狀態。
之後安裝 4 號軟體包,需要安裝 1,4 兩個軟體包。此時 0,1,4 處在安裝狀態。
最後,解除安裝 0 號軟體包會解除安裝所有的軟體包。
 
 
n=100000
q=100000   題意: 其實就是給你兩個操作: 1.。一條鏈權值全部變成1 2.一個子樹權值全部變成0 問每次操作都有多少個點被改變了   思路: 比較裸的樹鏈剖分+線段樹,因為是修改後直接詢問,我們可以直接把修改和詢問操作合併,詢問鏈的時候因為要在不同鏈上跳來跳去我們需要加個solve函式,子樹的話因為樹鏈剖分的特性一個子樹在剖下來的區間中是連續的區間,我們直接區間詢問就好了。 忘了自己是開的雙向邊。。陣列開1e5,RE了一發,自閉。 實現程式碼:
#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define mid int m = (l + r) >> 1
const int mod = 201314;
const int M = 2e5 + 10;  //建了雙向邊 *個2
struct node{
    int to,next;
}e[M];
int cnt,cnt1,n;
int son[M],siz[M],head[M],fa[M],top[M],dep[M],tid[M],mx[M],rk[M];
int sum[M<<2],lazy[M<<2];
void add(int u,int v){
    e[++cnt].to = v;e[cnt].next = head[u];head[u] = cnt;
}

void dfs1(int u,int faz,int deep){
    dep[u] = deep;
    fa[u] = faz;
    siz[u] = 1;
    for(int i = head[u];i;i = e[i].next){
        int v = e[i].to;
        if(v == faz) continue;
        dfs1(v,u,deep+1);
        siz[u] += siz[v];
        if(siz[v] > siz[son[u]]||son[u] == -1)
            son[u] = v;
    }
}

void dfs2(int u,int t){
    top[u] = t;
    mx[u] = cnt1;
    tid[u] = cnt1;
    rk[cnt1] = u;
    cnt1++;
    if(son[u] == -1) return ;
    dfs2(son[u],t),mx[u] = max(mx[u],mx[son[u]]);
    for(int i = head[u];i;i=e[i].next){
        int v = e[i].to;
        if(v != fa[u]&&v != son[u])
            dfs2(v,v),mx[u] = max(mx[u],mx[v]);
    }
}

void pushup(int rt){
    sum[rt] = sum[rt<<1] + sum[rt<<1|1];
}

void pushdown(int l,int r,int rt){
    if(lazy[rt]==0){
        sum[rt<<1] = sum[rt<<1|1] = 0;
        lazy[rt<<1] = lazy[rt<<1|1] = 0;
        lazy[rt] = -1;
    }
    else if(lazy[rt] == 1){
        mid;
        sum[rt<<1] = m-l+1  ;
        sum[rt<<1|1] = r-m;
        lazy[rt<<1] = lazy[rt<<1|1] = 1;
        lazy[rt] = -1;
    }
}

void build(int l,int r,int rt){
    lazy[rt] = -1;
    if(l == r){
        sum[rt] = 0;
        lazy[rt] = -1;
        return ;
    }
    mid;
    build(lson); build(rson);
}

int update(int L,int R,int c,int l,int r,int rt){
   if(L <= l&&R >= r){
       if(c==0){
           int cnt = sum[rt];
           sum[rt] = 0; lazy[rt] = 0;
           return cnt;
       }
       else{
            int cnt = r-l+1-sum[rt];
            sum[rt] = r-l+1;
            lazy[rt] = 1;
            return cnt;
       }
   }
   pushdown(l,r,rt);
   mid; int ret = 0;
   if(L <= m) ret += update(L,R,c,lson);
   if(R > m) ret += update(L,R,c,rson);
   pushup(rt);
   return ret;
}

int solve(int x,int y){
    int fx = top[x],fy = top[y];
    int ans = 0;
    while(fx != fy){
        if(dep[fx] < dep[fy]) swap(fx,fy),swap(x,y);
        ans += update(tid[fx],tid[x],1,1,n,1);
        x = fa[fx]; fx = top[x];
    }
    if(dep[x] > dep[y]) swap(x,y);
    ans += update(tid[x],tid[y],1,1,n,1);
    return ans;
}

char s[100];
int main()
{
    int x,m;
    cnt = 1,cnt1 = 1;
    scanf("%d",&n);
    memset(son,-1,sizeof(son));
    for(int i = 2;i <= n;i ++){
        scanf("%d",&x);
        add(x+1,i);add(i,x+1);
    }
    dfs1(1,0,1); dfs2(1,0);
    build(1,n,1);
    scanf("%d",&m);
    for(int i = 1;i <= m;i ++){
        scanf("%s",s);
        if(s[0] == 'i'){
            scanf("%d",&x);
            x++;
            printf("%d\n",solve(x,1));
        }
        else{
            scanf("%d",&x);
            x++;
            printf("%d\n",update(tid[x],mx[x],0,1,n,1));
        }
    }
    return 0;
}