通過這道題了解到了神奇的樹上dsu演算法，起因是最近想再學習一波線段樹，然後就找了一套線段樹總結的題目。在那個總結裡面看到了這道題，感覺很有意思。然後晚上回了宿舍就和室友在討論，然後室友告訴我樹上啟發式合併隨便做？？？

當時一臉懵逼，樹上啟發式合併是什麼鬼，其實個人對啟發式合併還是有一點了解的。之前做了一道可持久化並查集，那裡面並查集的合併沒有用正常的路徑壓縮，是用判斷兩個並查集大小的方式進行合併，最後保證複雜度維持在nlogn。

但是樹上啟發式合併之前完全沒有聽說過。室友信誓旦旦的說這做法可以解決大部分的子樹查詢問題。表示懷疑，因為我最開始想到的解法其實是dfs序之後直接莫隊。線段樹我們都不懂怎麼寫= = 然後一起研究了一下線段樹的寫法，感覺確實太麻煩了，室友說他第二天要用dsu解決這道題。我第二天用了莫隊A了，他用dsuA了。。。於是我就把他的資料拿過來研究了一下，發現樹上dsu是真的神奇，搞懂以後就感覺是個純暴力但是稍加修改後它的複雜度是nlogn！！！

但是這個演算法的侷限性也比較大，首先好像不能插入修改操作，其次的話只能應對子樹的查詢問題。而且怎麼說呢，如果要查詢的東西比較複雜的話感覺也比較難搞。

不過就算如此，也是一個非常優秀的演算法了，在應對子樹查詢類問題的時候可以優先考慮一下。而且，技多不壓身啊。。。

貼一下樹上dsu的學習資料：

原文連結：CF原文

題目大意：

給你一顆樹，每個結點有自己的權值，求權值恰好出現次數為k的個數。

解題思路：

dfs序後直接莫隊或者線段樹都可，這裡用樹上dsu來寫 = = 

Ac程式碼：

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=1e5+5;
int n,k,m,res,a[maxn],sum[maxn],sz[maxn],son[maxn],q[maxn],ans[maxn];
bool vis[maxn];
vector<int> v[maxn],rv;
int getid(int x) { return lower_bound(rv.begin(),rv.end(),x)-rv.begin()+1; }
void init()
{
    rv.clear();
    for(int i=0;i<=n;i++) v[i].clear();
    for(int i=0;i<=n;i++) sum[i]=son[i]=vis[i]=0;
}
void dfs1(int x,int fa) //求出每個結點的重兒子
{
    sz[x]=1;
    for(int i=0;i<v[x].size();i++)
    {
        int u=v[x][i];
        if(u==fa) continue;
        dfs1(u,x); sz[x]+=sz[u];
        if(sz[u]>sz[son[x]]) son[x]=u;
    }
}
void edit(int x,int fa,int sp)  //更新每個點的貢獻
{
    if(sum[a[x]]==k) res--;
    sum[a[x]]+=sp;
    if(sum[a[x]]==k) res++;
    for(int i=0;i<v[x].size();i++)
    {
        int u=v[x][i];
        if(u==fa||vis[u]) continue;
        edit(u,x,sp);
    }
}
void dfs(int x,int fa,int sp)   //sp表示貢獻是否被清空
{
    for(int i=0;i<v[x].size();i++)  //遍歷所有非重兒子的結點
    {
        int u=v[x][i];
        if(u==fa||u==son[x]) continue;
        dfs(u,x,0); //sp傳0
    }
    if(son[x]) dfs(son[x],x,1),vis[son[x]]=1;   //sp傳1
    edit(x,fa,1);   //遍歷x的所有結點 新增重兒子以外的貢獻
    ans[x]=res; //記錄答案
    if(son[x]) vis[son[x]]=0;   //重兒子的標記取消
    if(!sp) edit(x,fa,-1);  //貢獻取消
}
int main()
{
    int QAQ,kase=0,flag=0;
    scanf("%d",&QAQ);
    while(QAQ--)
    {
        if(flag) printf("\n");
        scanf("%d%d",&n,&k);
        init();
        for(int i=1;i<=n;i++)   //這裡把權值離散一下
            scanf("%d",&a[i]),rv.push_back(a[i]);
        sort(rv.begin(),rv.end()),rv.erase(unique(rv.begin(),rv.end()),rv.end());
        for(int i=1;i<=n;i++) a[i]=getid(a[i]);
        for(int i=1;i<n;i++)
        {
            int x,y;
            scanf("%d%d",&x,&y);
            v[x].push_back(y);
            v[y].push_back(x);
        }
        scanf("%d",&m);
        for(int i=1;i<=m;i++) scanf("%d",&q[i]);    //儲存詢問
        dfs1(1,0);dfs(1,0,0);   //遍歷跑出結果
        printf("Case #%d:\n",++kase);
        for(int i=1;i<=m;i++)
            printf("%d\n",ans[q[i]]);
        flag=1;
    }
}