jzoj3910-Idiot的间谍网络【倍增,dfs】
生活随笔
收集整理的這篇文章主要介紹了
jzoj3910-Idiot的间谍网络【倍增,dfs】
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
正題
題目鏈接:https://jzoj.net/senior/#main/show/3910
題目大意
nnn個點的森林mmm個操作,
解題思路
這里考慮離線的算法。我們可以先將這個森林的最終狀態構建出來,每條邊的邊權表示上一個標記的編號。然后我們發現個查詢需要滿足兩個條件,
這里用樹上倍增的方法計算即可。時間復雜度O(nlog?n)O(n\log n)O(nlogn)。加上快讀可以卡過
正解好像是O(α(n)n)O(\alpha (n)n)O(α(n)n)搜索樹這里就不具體寫了。
codecodecode
#include<cstdio> #include<cstring> #include<algorithm> #include<cctype> #include<vector> using namespace std; const int N=5e5+100; struct node{int to,next,w; }a[N]; int tot,cnt,z,num,n,m,f[N][20],v[N][20],dep[N]; int ls[N],rfn[N],ed[N],p[N],ask[N],in[N]; bool ans[N]; vector<int> q[N]; void addl(int x,int y,int w) {a[++tot].to=y;a[tot].next=ls[x];ls[x]=tot;a[tot].w=w; } int read() {int x=0,f=1; char c=getchar();while(!isdigit(c)) {if(c=='-')f=-f;c=getchar();}while(isdigit(c)) x=(x<<1)+(x<<3)+c-48,c=getchar();return x*f; } void dfs(int x) {rfn[x]=++cnt;for(int i=ls[x];i;i=a[i].next){int y=a[i].to;dep[y]=dep[x]+1;f[y][0]=x;v[y][0]=a[i].w;dfs(y);}ed[x]=cnt;return; } int LCA(int x,int y) {int z=0;for(int i=19;i>=0;i--)if(dep[f[y][i]]>=dep[x])z=max(v[y][i],z),y=f[y][i];return z; } int main() {//freopen("data.in","r",stdin);//freopen("data.out","w",stdout);n=read();m=read();for(int i=1;i<=m;i++){int op,x,y;op=read();if(op==1){x=read();y=read();in[x]++;addl(y,x,z);}if(op==2)p[++z]=read();if(op==3){ask[++num]=read();q[read()].push_back(num);}}for(int i=1;i<=n;i++)if(!in[i]) dep[i]=1,dfs(i);for(int j=1;j<20;j++)for(int i=1;i<=n;i++)f[i][j]=f[f[i][j-1]][j-1],v[i][j]=max(v[i][j-1],v[f[i][j-1]][j-1]);for(int i=1;i<=z;i++){for(int j=0;j<q[i].size();j++){int x=ask[q[i][j]],y=p[i];if(rfn[x]<=rfn[y]&&ed[x]>=rfn[y]&&LCA(x,y)<i)ans[q[i][j]]=1;}}for(int i=1;i<=num;i++)if(ans[i]) puts("YES");else puts("NO"); }總結
以上是生活随笔為你收集整理的jzoj3910-Idiot的间谍网络【倍增,dfs】的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: jzoj3920-噪音【贪心,dp】
- 下一篇: jzoj4815-ksum【堆】