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+7,mod=571111;
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=1; x=0;
    while(ch!='-'&&(ch<'0'||ch>'9')) ch=getchar();
    if(ch=='-') f=-1,ch=getchar();
    for(;ch>='0'&&ch<='9';ch=getchar()) x=x*10+ch-'0'; 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
",u,v,w);
    }
    LL f[N];
    int que[N];
    void dfs(int x,int fa) {
        f[x]=0;
        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=1,qr=0;
            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
",u,v,w);
    }
    
    int R[N],f[N][20]; LL H[N];
    void lca(int &cc,int x,int y) {
        if(R[x]<R[y]) swap(x,y);
        Rep(i,18,0) if(R[f[x][i]]>=R[y])
            x=f[x][i];
        if(x==y) { a[++cc]=x; return ; }
        Rep(i,18,0) if(f[x][i]!=f[y][i])
            x=f[x][i],y=f[y][i];
        if(f[x][0]>n) {
            a[++cc]=x; a[++cc]=y;
        }
        a[++cc]=f[x][0]; return ;
    }
    
    int dfk,sz[N];
    int dfs(int x,int fa) {
        sz[x]=1;
        f[x][0]=fa;
        R[x]=R[fa]+1;
        dfn[x]=++dfk;
        For(i,1,18) f[x][i]=f[f[x][i-1]][i-1];
        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=0;
        read(cnt); sz=cnt;
        For(i,1,cnt) { read(a[i]); isnode[a[i]]=1; }
        sort(a+1,a+cnt+1,cmp);
        For(i,1,cnt-1) {
            lca(sz,a[i],a[i+1]);
        }
        sort(a+1,a+sz+1,cmp);
        cnt=unique(a+1,a+sz+1)-(a+1);
        For(i,1,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=0; V.dfs(rt,0);
        printf("%lld
",ans);
        For(i,1,cnt) { isnode[a[i]]=0; V.fir[a[i]]=0; } V.ecnt=0;
    }
}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^1]=1;
            if(!dfn[to[i]]) {
                tarjan(to[i]);
                low[x]=min(low[x],low[to[i]]);
                if(low[to[i]]>=dfn[x]) {
                    LL hc=0,nc=0,cc=0;
                    T.add(x,++tot,0);
                    Rep(j,top,1) {
                         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==1) 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=1; 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*2333+v)%mod,up=vc[hs].size();
        For(i,0,up-1) {
            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,1,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,1,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(1); 
    T.dfs(1,0);
    read(Q);
    For(i,1,Q) T.solve();
    Formylove;
}

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

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

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

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

http://uoj.ac/problem/87