bzoj 1023: [SHOI2008]cactus仙人掌圖【tarjan+dp+單調隊列】
阿新 • • 發佈:2019-02-28
最短距離 esp main 仙人掌 i++ clu pac tarjan 處理
本來想先求出點雙再一個一個處理結果寫了很長發現太麻煩
設f[u]為u點向下的最長鏈
就是再tarjan的過程中,先照常處理,用最長兒子鏈和次長兒子鏈更新按ans,然後處理以這個點為根的環,也就是這個點是dfs第一次到這個環訪問到的點
環用來更新ans的是兒子鏈+到根的一段,這個直接for一邊就行,還有就是一個兒子鏈+環上的一段+另一個兒子鏈,這個把環復制一遍然後單調隊列掃即可
註意距離的定義是最短距離!
#include<iostream> #include<cstdio> using namespace std; const int N=100005; int n,m,h[N],cnt,fa[N],de[N],f[N],s[N],dfn[N],low[N],dft,ans,a[N],q[N]; struct qwe { int ne,to; }e[N*10]; int read() { int r=0,f=1; char p=getchar(); while(p>'9'||p<'0') { if(p=='-') f=-1; p=getchar(); } while(p>='0'&&p<='9') { r=r*10+p-48; p=getchar(); } return r*f; } void add(int u,int v) { cnt++; e[cnt].ne=h[u]; e[cnt].to=v; h[u]=cnt; } void dfs(int u,int fat) { dfn[u]=low[u]=++dft; fa[u]=fat; de[u]=de[fat]+1; for(int i=h[u];i;i=e[i].ne) if(e[i].to!=fat) { if(!dfn[e[i].to]) { dfs(e[i].to,u); low[u]=min(low[u],low[e[i].to]); } else low[u]=min(low[u],dfn[e[i].to]); if(low[e[i].to]>dfn[u]) { ans=max(ans,f[u]+f[e[i].to]+1); f[u]=max(f[u],f[e[i].to]+1); } } for(int j=h[u];j;j=e[j].ne) if(fa[e[j].to]!=u&&dfn[e[j].to]>dfn[u]) { int tot=de[e[j].to]-de[u]+1,l=1,r=1; q[1]=1; for(int x=e[j].to;x!=u;x=fa[x]) a[tot--]=f[x]; a[tot]=f[u]; tot=de[e[j].to]-de[u]+1; for(int i=1;i<=tot;i++) a[i+tot]=a[i]; for(int i=2;i<=2*tot;i++) { while(l<=r&&i-q[l]>tot/2) l++; if(l<=r) ans=max(ans,a[i]+i+a[q[l]]-q[l]); while(l<=r&&a[q[r]]-q[r]<=a[i]-i) r--; q[++r]=i; } for(int i=2;i<=tot;i++) f[u]=max(f[u],a[i]+min(i-1,tot-i+1)); } } int main() { n=read(),m=read(); for(int i=1;i<=m;i++) { int k=read(),la=read(),x; for(int j=2;j<=k;j++) { x=read(); add(la,x),add(x,la); la=x; } } dfs(1,0); printf("%d\n",ans); return 0; }
bzoj 1023: [SHOI2008]cactus仙人掌圖【tarjan+dp+單調隊列】