【JZOJ4784】Map

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 long
using 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;
}