【JZOJ 4784】Map

Description

upd:保证原图连通。
“不相交路径”的定义为不存在相同的边。可以存在相同的点。重边视为不同的边。

Solution

这个肯定是先tarjan缩点，
缩完点以后出来的肯定是一棵树，
推一下公式，做一个前缀和，LCA一下即可

复杂度：$O(n\log(n))$

Code

#include<iostream>
#include<cstdio>
#include<cstdlib>
#define efo(i,q) for(int i=A[q];i;i=B[i][0])
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fod(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
typedef long long LL;
const int N=500500,M=25;
int read(int &n)
{
    char ch=' ';int q=0,w=1;
    for(;(ch!='-')&&((ch<'0')||(ch>'9'));ch=getchar());
    if(ch=='-')w=-1,ch=getchar();
    for(;ch>='0' && ch<='9';ch=getchar())q=q*10+ch-48;n=q*w;return n;
}
int m,n,m1,H,w;
LL ans;
int g[N],c[N];
LL f[N],su[N],b[N];
int A[N],B[2*N][3],HE[N],B0=1;
bool z[N];
int za[N],dfn[N],low[N];
int G[N][M+1];
void link(int q,int w)
{
    if(!A[q])HE[q]=B0+1;B[++B0][0]=A[q],B[A[q]=B0][1]=w;
    if(!A[w])HE[w]=B0+1;B[++B0][0]=A[w],B[A[w]=B0][1]=q;
}
int tarjan(int q)
{
    za[dfn[q]=low[q]=++za[0]]=q;z[q]=1;
    efo(i,q)if(!B[i][2])
        if(!z[B[i][1]])B[i^1][2]=1,low[q]=min(low[q],tarjan(B[i][1])),B[i^1][2]=0;
        else low[q]=min(low[q],low[B[i][1]]);
    if(dfn[q]==low[q])
    {
        g[q]=q;z[q]=0;
        while((w=za[za[0]--])!=q)g[w]=q,z[w]=0;
    }
    return low[q];
}
void dfs(int q,int fa)
{
    G[q][0]=fa;c[q]=c[fa]+1;
    su[q]=su[fa]+b[q];f[q]=f[fa]+b[q]*su[fa];
    efo(i,q)if(B[i][1]!=fa&&!su[B[i][1]])dfs(B[i][1],q);
}
int LCAO(int q,int e)
{
    while(c[q]>e)
    {
        int I=0;
        while(c[G[q][I+1]]>e)I++;
        q=G[q][I];
    }
    return q;
}
int LCA(int q,int w)
{
    q=LCAO(q,c[w]);
    w=LCAO(w,c[q]);
    while(q!=w)
    {
        int I=0;
        while(G[q][I+1]!=G[w][I+1])I++;
        q=G[q][I];w=G[w][I];
    }
    return q;
}
int main()
{
    int q,w;
    read(n),read(m),read(m1);
    fo(i,1,m)read(q),read(w),link(q,w);
    tarjan(1);
    fo(i,1,n)
    {
        efo(j,i)B[i][1]=g[B[i][1]];
        if(g[i]!=i)B[HE[i]][0]=A[g[i]],A[g[i]]=A[i],A[i]=0;
        b[g[i]]++;
    }
    dfs(1,0);
    fo(j,1,M)fo(i,1,n)G[i][j]=G[G[i][j-1]][j-1];
    fo(i,1,m1)
    {
        read(q),read(w);
        q=g[q];w=g[w];
        if(q!=w)
        {
            int t=LCA(q,w);
            ans+=(su[q]-su[t])*(su[w]-su[t])+(f[q]-f[t]-su[G[t][0]]*(su[q]-su[t]))+(f[w]-f[t]-su[G[t][0]]*(su[w]-su[t]));
        }   
    }
    printf("%lld\n",ans*2);
    return 0;
}