【JZOJ4784】Map

来源：互联网发布：淘宝天猫人工投诉电话编辑：程序博客网时间：2024/06/05 03:54

Description

这里写图片描述

Solution

首先我们把边双缩成点，设每个点的权值是边双的大小，原图会变成一棵树。

对于增加一条道路，就会形成一个环。那么增加的点对个数就是：点权和的平方减去点权的平方和。

那么只要在树上倍增即可。

Code

#include<iostream>#include<cstdio>#include<cstdlib>#include<cmath>#define fo(i,j,k) for(int i=j;i<=k;i++)#define fd(i,j,k) for(int i=j;i>=k;i--)#define N 200001#define M 800001#define ll long longusing namespace std;int to[M],next[M],last[M],num=0;int t1[M],nx[M],ls[M],cnt=0;int b[N],c[N],tot=0;int d[N];int f[N][19];ll g[N][19];ll z[N][19];bool bz[M];bool vis[N];int low[N],dfn[N],st[N],top=0,dep=0;struct node{    int x,y;}e[M];void link(int x,int y){    num++;    to[num]=y;    next[num]=last[x];    last[x]=num;}void link1(int x,int y){    cnt++;    t1[cnt]=y;    nx[cnt]=ls[x];    ls[x]=cnt;}ll sqr(ll x){    return x*x;}int turn(int x){    return x%2?x+1:x-1;}void tarjan(int x){    low[x]=dfn[x]=++dep;    st[++top]=x;    vis[x]=true;    for(int i=last[x];i;i=next[i])    if(!bz[i])    {        int v=to[i];        if(!dfn[v])        {            bz[i]=bz[turn(i)]=true;            tarjan(v);            bz[i]=bz[turn(i)]=false;            low[x]=min(low[x],low[v]);        }        else if(vis[v]) low[x]=min(low[x],dfn[v]);    }    if(low[x]==dfn[x])    {        tot++;        while(st[top+1]!=x && top)        {            c[st[top]]=tot;            b[tot]++;            vis[st[top]]=false;            top--;        }    }}void find(int x){    //cout<<x<<endl;    //system("pause");    for(int i=ls[x];i;i=nx[i])    {        int v=t1[i];        if(v!=f[x][0])        {            f[v][0]=x;            g[v][0]=b[v];            z[v][0]=b[v]*1ll*b[v];            d[v]=d[x]+1;            find(v);        }    }}ll get(int u,int v){    if(d[v]>d[u]) swap(u,v);    ll tg=0,tz=0;    fd(i,18,0)    if(f[u][i] && d[f[u][i]]>=d[v]) tg+=g[u][i],tz+=z[u][i],u=f[u][i];    if(u==v)    {        tg+=b[u];        tz+=sqr(b[u]);        return sqr(tg)-tz;    }    fd(i,18,0)    if(f[u][i] && f[v][i] && f[u][i]!=f[v][i])    {        tg+=g[u][i]+g[v][i];        tz+=z[u][i]+z[v][i];        u=f[u][i];        v=f[v][i];    }    tg+=g[u][0]+g[v][0];    tz+=z[u][0]+z[v][0];    u=f[u][0];v=f[v][0];    tg+=b[u];    tz+=sqr(b[u]);    return sqr(tg)-tz;}ll ans=0;int main(){    int n,m,q;    cin>>n>>m>>q;    fo(i,1,m)    {        int x,y;        scanf("%d %d",&x,&y);        link(x,y);        link(y,x);        e[i].x=x;        e[i].y=y;    }    tarjan(1);    int p=0;    fo(i,1,m)    {        int x=c[e[i].x],y=c[e[i].y];        if(x!=y)        {            link1(x,y);            link1(y,x);        }    }    d[0]=-1;    g[c[1]][0]=b[c[1]];    z[c[1]][0]=b[c[1]]*1ll*b[c[1]];    find(c[1]);    fo(j,1,18)    fo(i,1,n)    {        f[i][j]=f[f[i][j-1]][j-1];        g[i][j]=g[i][j-1]+g[f[i][j-1]][j-1];        z[i][j]=z[i][j-1]+z[f[i][j-1]][j-1];    }    while(q--)    {        int x,y;        scanf("%d %d",&x,&y);        ans+=get(c[x],c[y]);    }    cout<<ans;}

1 0