uoj#87. mx的仙人掌

 //Achen
 #include<bits/stdc++.h>
 #define For(i,a,b) for(int i=(a);i<=(b);i++)
 #define Rep(i,a,b) for(int i=(a);i>=(b);i--)
 #define Formylove return 0
 const int N=2e6+,mod=;
 typedef long long LL;
 typedef double db;
 using namespace std;
 int n,m,Q;

 template<typename T> void read(T &x) {
     char ch=getchar(); T f=; x=;
     while(ch!='-'&&(ch<''||ch>'')) ch=getchar();
     if(ch=='-') f=-,ch=getchar();
     for(;ch>=''&&ch<='';ch=getchar()) x=x*+ch-''; x*=f;
 }

 LL ans,Hdis[N];
 int isnode[N];
 struct VMtree{
     int ecnt,fir[N],nxt[N],to[N]; LL val[N];
     void add(int u,int v,LL w) {
         nxt[++ecnt]=fir[u]; fir[u]=ecnt; to[ecnt]=v; val[ecnt]=w;
         nxt[++ecnt]=fir[v]; fir[v]=ecnt; to[ecnt]=u; val[ecnt]=w;
         //printf("%d %d %lld\n",u,v,w);
     }
     LL f[N];
     int que[N];
     void dfs(int x,int fa) {
         f[x]=;
         for(int i=fir[x];i;i=nxt[i]) if(to[i]!=fa)
             dfs(to[i],x);
         if(x<=n) {
             for(int i=fir[x];i;i=nxt[i]) if(to[i]!=fa) {
                 int y=to[i];
                 ans=max(ans,f[y]+val[i]+f[x]); f[x]=max(f[x],f[y]+val[i]);
             }
         }
         else {
             int ql=,qr=;
             for(int i=fir[x];i;i=nxt[i]) if(to[i]!=fa) {
                 int y=to[i];
                 f[x]=max(f[x],f[y]+val[i]);
                 while(ql<=qr&&Hdis[y]-Hdis[que[ql]]>Hdis[x]-(Hdis[y]-Hdis[que[ql]])) ql++;
                 if(ql<=qr) ans=max(ans,Hdis[y]-Hdis[que[ql]]+f[que[ql]]+f[y]);
                 while(ql<=qr&&f[y]-Hdis[y]>=f[que[qr]]-Hdis[que[qr]]) qr--;
                 que[++qr]=y;
             }
             for(int i=fir[x];i;i=nxt[i]) if(to[i]!=fa) {
                 int y=to[i];
                 while(ql<=qr&&Hdis[x]-(Hdis[que[ql]]-Hdis[y])>Hdis[que[ql]]-Hdis[y]) ql++;
                 if(ql<=qr) ans=max(ans,Hdis[x]-(Hdis[que[ql]]-Hdis[y])+f[que[ql]]+f[y]);
             }
         } if(isnode[x]) ans=max(ans,f[x]);
     }
 }V;

 int tot,dfn[N],a[N];
 bool cmp(const int &A,const int &B) { return dfn[A]<dfn[B]; }
 struct Tree {
     int ecnt,fir[N],nxt[N],to[N]; LL val[N];
     void add(int u,int v,LL w) {
         nxt[++ecnt]=fir[u]; fir[u]=ecnt; to[ecnt]=v; val[ecnt]=w;
         nxt[++ecnt]=fir[v]; fir[v]=ecnt; to[ecnt]=u; val[ecnt]=w;
         //printf("%d %d %lld\n",u,v,w);
     }

     int R[N],f[N][]; LL H[N];
     void lca(int &cc,int x,int y) {
         if(R[x]<R[y]) swap(x,y);
         Rep(i,,) if(R[f[x][i]]>=R[y])
             x=f[x][i];
         if(x==y) { a[++cc]=x; return ; }
         Rep(i,,) if(f[x][i]!=f[y][i])
             x=f[x][i],y=f[y][i];
         if(f[x][]>n) {
             a[++cc]=x; a[++cc]=y;
         }
         a[++cc]=f[x][]; return ;
     }

     int dfk,sz[N];
     int dfs(int x,int fa) {
         sz[x]=;
         f[x][]=fa;
         R[x]=R[fa]+;
         dfn[x]=++dfk;
         For(i,,) f[x][i]=f[f[x][i-]][i-];
         for(int i=fir[x];i;i=nxt[i]) if(to[i]!=fa) {
             H[to[i]]=H[x]+val[i];
             dfs(to[i],x);
             sz[x]+=sz[to[i]];
         }
     }

     bool in(int a,int b) { return dfn[a]>=dfn[b]&&dfn[a]<dfn[b]+sz[b]; }
     int sta[N],top;
     void solve() {
         int cnt,sz,rt=;
         read(cnt); sz=cnt;
         For(i,,cnt) { read(a[i]); isnode[a[i]]=; }
         sort(a+,a+cnt+,cmp);
         For(i,,cnt-) {
             lca(sz,a[i],a[i+]);
         }
         sort(a+,a+sz+,cmp);
         cnt=unique(a+,a+sz+)-(a+);
         For(i,,cnt) {
             if(!rt||R[a[i]]<R[rt]) rt=a[i];
             while(top&&!in(a[i],sta[top])) top--;
             if(top) V.add(sta[top],a[i],H[a[i]]-H[sta[top]]);
             sta[++top]=a[i];
         } while(top) top--;
         ans=; V.dfs(rt,);
         printf("%lld\n",ans);
         For(i,,cnt) { isnode[a[i]]=; V.fir[a[i]]=; } V.ecnt=;
     }
 }T;

 struct Cactus {
     int ecnt,fir[N],nxt[N],to[N],val[N],vis[N];
     void add(int u,int v,int w) {
         nxt[++ecnt]=fir[u]; fir[u]=ecnt; to[ecnt]=v; val[ecnt]=w;
         nxt[++ecnt]=fir[v]; fir[v]=ecnt; to[ecnt]=u; val[ecnt]=w;
     }

     int dfk,dfn[N],low[N],sta[N],top;
     void tarjan(int x) {
         dfn[x]=low[x]=++dfk;
         for(int i=fir[x];i;i=nxt[i]) if(!vis[i]) {
             sta[++top]=i;
             vis[i]=vis[i^]=;
             if(!dfn[to[i]]) {
                 tarjan(to[i]);
                 low[x]=min(low[x],low[to[i]]);
                 if(low[to[i]]>=dfn[x]) {
                     LL hc=,nc=,cc=;
                     T.add(x,++tot,);
                     Rep(j,top,) {
                          hc+=val[sta[j]]; cc++;
                         if(sta[j]==i) break;
                     }
                     Hdis[tot]=hc;
                     while(top) {
                         int j=sta[top--];
                         if(to[j]!=x) {
                             if(cc==) T.add(tot,to[j],hc);
                             else T.add(tot,to[j],min(hc-nc,nc));
                             Hdis[to[j]]=nc;
                         }
                         nc+=val[j];
                         if(j==i) break;
                     }
                 }
             }
             else low[x]=min(low[x],dfn[to[i]]);
         }
     }

     void init() { ecnt=; tot=n; }
 }C;

 struct hash {
     int ecnt;
     struct edge{ int u,v,w; }e[N];
     vector<int>vc[N];
     void ins(int u,int v,int w) {
         if(u>v) swap(u,v);
         int hs=(u*+v)%mod,up=vc[hs].size();
         For(i,,up-) {
             int ec=vc[hs][i];
             if(e[ec].u==u&&e[ec].v==v) { e[ec].w=min(e[ec].w,w); return; }
         }
         e[++ecnt]=(edge){u,v,w}; vc[hs].push_back(ecnt);
     }

     void add() {
         For(i,,ecnt) C.add(e[i].u,e[i].v,e[i].w);
     }
 }H;

 int main() {
     //freopen("1.in","r",stdin);
     //freopen("1.out","w",stdout);
     read(n); read(m);
     For(i,,m) {
         int u,v,w;
         read(u); read(v); read(w);
         C.add(u,v,w);
         H.ins(u,v,w);
     }
     C.init(); H.add();
     C.tarjan();
     T.dfs(,);
     read(Q);
     For(i,,Q) T.solve();
     Formylove;
 }

写完花了1.5h找bug发现是数组开小了。

其实是因为太长实在不想看一直各种磨蹭各种骚扰其他同学。

无脑圆方树+虚树+单调队列优化dp。

如果两个人的lca是方点要把方点的两个圆儿子加进去。

http://uoj.ac/problem/87