傳送門

文章目錄

• 題意：
• 思路：

題意：

思路：

由于$a,b$都比$c$厲害，那么$a,b$一定是某個是某個的祖先。那么就分為兩種情況了：
$(1)$ $b$在$a$上面，約定$depth[1]=1$，此時答案顯然為$min(depth[a]?1,k)?(se[a]?1)$。
$(2)$$b$在$a$的下下面，這個時候就不是那么容易搞了，問題轉(zhuǎn)化成我們要求$a$這個子樹中$depth$范圍在$[depth[a]+1,depth[a]+k]$內(nèi)的所有點的$se[i]?1$，先考慮暴力怎么寫呢？顯然我們可以暴力對這顆子樹的深度建線段樹，這樣就變成了區(qū)間查詢的問題了。而這樣復雜度是肯定不行的，所以我們考慮建可持久化線段樹，根據(jù)樹的$dfs$序建主席樹，以節(jié)點深度為下標，那么查詢就變成了$query(root[dfn[p]],root[dfn[p]+se[p]?1],1,n,depth[p]+1,min(depth[p]+k,n))$，再加上$(1)$的答案即可。

我們還可以將其轉(zhuǎn)換成二維數(shù)點的問題。
通過以上分析不難發(fā)現(xiàn)我們要找的點就是
深度范圍是$[depth[p]+1,depth[p]+k]$，$dfs$序范圍是$[dfn[p],dfn[p]+se[p]?1]$，那么我們以深度建立$x$軸，以$dfs$序建立$y$軸，讓后統(tǒng)計就好啦。

主席樹$O(nlogn)$

//#pragma GCC optimize(2) #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("---") using namespace std;//void rd_cre() { freopen("d://dp//data.txt","w",stdout); srand(time(NULL)); } //void rd_ac() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//AC.txt","w",stdout); } //void rd_wa() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//WA.txt","w",stdout); }typedef long long LL; typedef unsigned long long ULL; typedef pair<int,int> PII;const int N=1000010,mod=1e9+7,INF=0x3f3f3f3f; const double eps=1e-6;int n,m; int root[N],tot,idx; int dfn[N],se[N],depth[N]; vector<int>v[N]; struct Node {int l,r;LL sum; }tr[N*40];void insert(int p,int &q,int l,int r,int pos,int x) {q=++tot; tr[q]=tr[p];tr[q].sum+=x;if(l==r) return;int mid=l+r>>1;if(pos<=mid) insert(tr[p].l,tr[q].l,l,mid,pos,x);else insert(tr[p].r,tr[q].r,mid+1,r,pos,x); }LL query(int p,int q,int l,int r,int ql,int qr) {if(l>=ql&&r<=qr) return tr[q].sum-tr[p].sum;LL ans=0;int mid=l+r>>1;if(ql<=mid) ans+=query(tr[p].l,tr[q].l,l,mid,ql,qr);if(qr>mid) ans+=query(tr[p].r,tr[q].r,mid+1,r,ql,qr);return ans; }void dfs1(int u,int fa) {se[u]=1; dfn[u]=++idx;depth[u]=depth[fa]+1;for(auto x:v[u]) if(x!=fa) dfs1(x,u),se[u]+=se[x]; }void dfs2(int u,int fa) {insert(root[dfn[u]-1],root[dfn[u]],1,n,depth[u],se[u]-1);for(auto x:v[u]) if(x!=fa) dfs2(x,u); }int main() { // ios::sync_with_stdio(false); // cin.tie(0);scanf("%d%d",&n,&m);for(int i=1;i<=n-1;i++){int a,b; scanf("%d%d",&a,&b);v[a].pb(b); v[b].pb(a);}dfs1(1,0); dfs2(1,0);while(m--){int p,k; scanf("%d%d",&p,&k);LL ans=1ll*min(depth[p]-1,k)*(se[p]-1);ans+=query(root[dfn[p]],root[dfn[p]+se[p]-1],1,n,depth[p]+1,min(depth[p]+k,n));printf("%lld\n",ans);}return 0; } /**/

總結(jié)

以上是生活随笔為你收集整理的P3899 [湖南集训]谈笑风生 主席树解决二维数点的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。