BZOJ1093 || 洛谷P2272 [ZJOI2007]最大半连通子图【Tarjan+DAG最长链计数】

Time Limit: 30 Sec
Memory Limit: 162 MB

Description

一个有向图G=(V,E)称为半连通的(Semi-Connected)，如果满足：?u,v∈V，满足u→v或v→u，即对于图中任意
两点u，v,存在一条u到v的有向路径或者从v到u的有向路径。若G’=(V’,E’)满足V’?V，E’是E中所有跟V’有关的边，
则称G’是G的一个导出子图。若G’是G的导出子图，且G’半连通，则称G’为G的半连通子图。若G’是G所有半连通子图
中包含节点数最多的，则称G’是G的最大半连通子图。给定一个有向图G，请求出G的最大半连通子图拥有的节点数K
，以及不同的最大半连通子图的数目C。由于C可能比较大，仅要求输出C对X的余数。

Input

第一行包含两个整数N，M，X。N，M分别表示图G的点数与边数，X的意义如上文所述接下来M行，每行两个正整
数a, b，表示一条有向边(a, b)。图中的每个点将编号为1,2,3…N，保证输入中同一个(a,b)不会出现两次。N ≤1
00000, M ≤1000000；对于100%的数据， X ≤10^8

Output

应包含两行，第一行包含一个整数K。第二行包含整数C Mod X.

---

题目分析

一开始搞错题意了还以为要有边直接相连。。。
Tarjan缩点求DAG最长链再计数即可
注意去重边

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<vector>
#include<queue>
using namespace std;
typedef long long lt;

int read()
{
    int f=1,x=0;
    char ss=getchar();
    while(ss<'0'||ss>'9'){if(ss=='-')f=-1;ss=getchar();}
    while(ss>='0'&&ss<='9'){x=x*10+ss-'0';ss=getchar();}
    return f*x;
}

const int maxn=2000010;
int n,m,mod;
struct edge{int u,v;}ed[maxn<<1];
struct node{int v,nxt;}E[maxn<<1],nE[maxn<<1];
int head[maxn],tot,nh[maxn];
int dfn[maxn],low[maxn],cnt,st[maxn],ins[maxn],top;
int col[maxn],coln,val[maxn];
int in[maxn],out[maxn];
lt dp[maxn],num[maxn];

bool cmp(edge a,edge b){ return a.u==b.u?a.v<b.v:a.u<b.u;}

void add(int u,int v,int d)
{
    if(d==0){
        E[++tot].nxt=head[u];
        E[tot].v=v;
        head[u]=tot;
    }
    else{
        nE[++tot].nxt=nh[u];
        nE[tot].v=v;
        nh[u]=tot;
    }
}

void tarjan(int u)
{
    dfn[u]=low[u]=++cnt;
    st[++top]=u; ins[u]=1;
    for(int i=head[u];i;i=E[i].nxt)
    {
        int v=E[i].v;
        if(!dfn[v]){ tarjan(v); low[u]=min(low[u],low[v]);}
        else if(ins[v]) low[u]=min(low[u],dfn[v]);
    }
    if(dfn[u]==low[u])
    {
        coln++;
        do{
            col[st[top]]=coln;
            ins[st[top]]=0;
            val[coln]++;
        }while(st[top--]!=u);
    }
}

lt DP(int u)
{
    if(!out[u]){
        num[u]=1;
        return dp[u]=val[u];
    }
    if(dp[u]) return dp[u];
    
    lt res=0;
    for(int i=nh[u];i;i=nE[i].nxt)
    {
        int v=nE[i].v;
        lt tt=DP(v);
        if(tt>res) res=tt,num[u]=num[v];
        else if(tt==res) num[u]=(num[u]+num[v])%mod;
    }
    return dp[u]=res+val[u];
}

int main()
{
    n=read();m=read();mod=read();
    for(int i=1;i<=m;++i)
    {
        int u=read(),v=read();
        add(u,v,0);
    }
    
    for(int i=1;i<=n;++i)
    if(!dfn[i]) tarjan(i);
    
    int sz=0; 
    for(int u=1;u<=n;++u)
    for(int i=head[u];i;i=E[i].nxt)
    {
        int v=E[i].v;
        if(col[u]!=col[v])
        ed[++sz]=(edge){col[u],col[v]};
    }
    
    tot=0; sort(ed+1,ed+1+sz,cmp);
    for(int i=1;i<=sz;++i)
    if(ed[i].u!=ed[i-1].u||ed[i].v!=ed[i-1].v)
    {
        add(ed[i].u,ed[i].v,1);
        ++out[ed[i].u]; ++in[ed[i].v];
    }
    
    lt ans=0,ct=0;
    for(int i=1;i<=coln;++i)
    if(!in[i]) ans=max(ans,DP(i));
    for(int i=1;i<=coln;++i)
    if(!in[i]&&dp[i]==ans) ct=(ct+num[i])%mod;
    
    printf("%lld\n%lld",ans,ct);
    return 0;
}