傳送門
文章目錄
題意：
思路：
這個題可以用各種姿勢a掉，樹啟和線段樹合并都可以，比較無腦。下面給一個解決這種問題比較巧妙的做法。
考慮暴力怎么寫，我們先將每個點的權值離散化一下，每到一個點就把樹狀數組清空，讓后在樹狀數組中加上這個子樹的值，用樹狀數組查詢pj>pip_j>p_ipj?>pi?的數量。
我們優化以上算法，還是先離散化，也是用到了樹狀數組。
考慮對于一個點，當dfsdfsdfs遍歷這個點會遍歷兩次 ，第一次是從父親節點到這個點，第二次是從他的子節點回溯回來。我們在回溯的時候加上這個節點的權值，那么當回溯到這個點的時候，樹狀數組中存的值包含了以這個點為根的子樹的全部權值，當然里面可能還有別的點的值(即"垃圾"值)，如果我們按照之前一遍一清空的話復雜度是O(n2)O(n^2)O(n2)的，我們就是要優化這里。現在想想，如果能把樹狀數組中那些”垃圾“值減去不就行了嗎？我們如何才能知道那些值呢？我們驚奇的發現第一次到這個點的時候樹狀數組中存的不就是這個值嘛？所以我們定義一個ansansans，第一次到這個點的時候減去sum(se)?sum(a[u])sum(se)-sum(a[u])sum(se)?sum(a[u])，回溯的時候加上sum(se)?sum(a[u])sum(se)-sum(a[u])sum(se)?sum(a[u])即可，這個時候即為答案啦。
#include<cstdio>
#include<iostream>
#include<string>
#include<cstring>
#include<map>
#include<cmath>
#include<cctype>
#include<vector>
#include<set>
#include<queue>
#include<algorithm>
#include<sstream>
#include<ctime>
#include<cstdlib>
#define X first
#define Y second
#define L (u<<1)
#define R (u<<1|1)
#define pb push_back
#define mk make_pair
#define Mid (tr[u].l+tr[u].r>>1)
#define Len(u) (tr[u].r-tr[u].l+1)
#define random(a,b) ((a)+rand()%((b)-(a)+1))
#define db puts("---")
#define lowbit(x) ((x)&(-x))
using namespace std
;
typedef long long LL
;
typedef unsigned long long ULL
;
typedef pair
<int,int> PII
;const int N
=1000010,M
=N
*2,mod
=1e9+7,INF
=0x3f3f3f3f;
const double eps
=1e-6;int n
,se
;
int a
[N
];
vector
<int>v
;
int tr
[N
];
int e
[M
],ne
[M
],h
[N
],idx
;
int ans
[N
];void add(int a
,int b
)
{e
[idx
]=b
,ne
[idx
]=h
[a
],h
[a
]=idx
++;
}void addd(int x
,int c
)
{for(int i
=x
;i
<=se
;i
+=lowbit(i
)) tr
[i
]+=c
;
}int sum(int x
)
{int ans
=0;for(int i
=x
;i
;i
-=lowbit(i
)) ans
+=tr
[i
];return ans
;
}int find(int x
)
{return lower_bound(v
.begin(),v
.end(),x
)-v
.begin();
}void dfs(int u
)
{ans
[u
]-=sum(se
)-sum(a
[u
]);for(int i
=h
[u
];~i
;i
=ne
[i
]) { int ver
=e
[i
]; dfs(ver
); }ans
[u
]+=sum(se
)-sum(a
[u
]);addd(a
[u
],1);
}int main()
{
scanf("%d",&n
);memset(h
,-1,sizeof(h
)); idx
=0;for(int i
=1;i
<=n
;i
++) scanf("%d",&a
[i
]),v
.pb(a
[i
]);sort(v
.begin(),v
.end()); v
.erase(unique(v
.begin(),v
.end()),v
.end());se
=v
.size();for(int i
=1;i
<=n
;i
++) a
[i
]=find(a
[i
])+1;for(int i
=2;i
<=n
;i
++){int x
; scanf("%d",&x
);add(x
,i
);}dfs(1);for(int i
=1;i
<=n
;i
++) printf("%d\n",ans
[i
]);return 0;
}
總結
以上是生活随笔為你收集整理的P3605 [USACO17JAN]Promotion Counting P dfs序的全部內容，希望文章能夠幫你解決所遇到的問題。
如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。
