有向图的强联通分量之:【求最长链】【求最长链的方案数(图论中的方案数DP)】【最长链和最大半联通子图 节点数相同】【最长链与最大半联通子图等价又不完全等价】

 注意了:

最长链与最大半联通子图是等价又不完全等价的。

最大半联通子图的节点数 会 与最长链的节点数相同,

但是边数的话,最大半联通子图还是会加进 两个点之间的所有边

但是本题的主要目的是求 最长链上的节点数,所以不要紧

 

 

#include <bits/stdc++.h>
using namespace std;
const int N=1e5+10,M=2e6+10;
typedef pair<int,int> PII;
typedef long long LL;
int n,m,mode;
int e[M],ne[M],h[N],idx,hs[N];
int dfn[N],low[N],timestamp;
int stk[N],top;
int id[N],siz[N],scc_cnt;
bool in_stk[N];
int f[N],g[N];

void add(int h[],int a,int b)
{
    e[idx]=b;
    ne[idx]=h[a];
    h[a]=idx++;
    
}


void tarjan(int u)
{
    dfn[u]=low[u]=++timestamp;
    stk[++top]=u,in_stk[u]=true;
    
    for(int i=h[u];i!=-1;i=ne[i])
    {
        int j=e[i];
        if(!dfn[j])
        {
            tarjan(j);
              low[u]=min(low[u],low[j]);
        }
        else if(in_stk[j])
        {
            low[u]=min(low[u],dfn[j]);
        }
   
    }
    
     
    if(dfn[u]==low[u])
    {
        ++scc_cnt;
        int y;
        do
        {
            y=stk[top--];
            in_stk[y]=false;
            id[y]=scc_cnt;
            siz[scc_cnt]++;
        }while(y!=u);
    }
}

int main()
{
    memset(h,-1,sizeof h);
    memset(hs,-1,sizeof hs);
    
    scanf("%d%d%d",&n,&m,&mode);
    
    
    while(m--)
    {
        int a,b;
        scanf("%d%d",&a,&b);
        add(h,a,b);
    }
    
    for(int i=1;i<=n;i++)
    {
        if(!dfn[i])
        {
            tarjan(i);
          
        }
    }
   
    
    
    //通过缩点的过程,把强连通分量这个新图(拓扑图)建立起来
    /*判边重,
    但是想一下,我们用什么来表示一个边呢?
    (用边上两端点的 编号)
    来表示,由于编号的范围在 1e5左右,所以
    我们可以用 u*1e6+v 来区分
    */
    unordered_set<LL> s;
    
    for(int i=1;i<=n;i++)
    for(int j=h[i];j!=-1;j=ne[j])  //枚举原图
    {
        int k=e[j];
        int a=id[i],b=id[k];
        LL hash=(LL)a*1000000+b;
        if(a!=b&&!s.count(hash))
        {
            add(hs,a,b);
            s.insert(hash);
        }
    }
    
    //枚举联通块,更新dp数组
    for(int i=scc_cnt;i;i--)
    {
        if(!f[i])
        {
            f[i]=siz[i];  //f[i]是以第i个点为终点的节点数量最大值
            g[i]=1;  //让f[i]取得最大值的方案数
        }
        
        for(int j=hs[i];j!=-1;j=ne[j])
        {
            int k=e[j];
            if(f[k]<f[i]+siz[k])
            {
                f[k]=f[i]+siz[k];
                g[k]=g[i];
            }
            else if(f[k]==f[i]+siz[k])
            {
                g[k]=(g[k]+g[i])%mode;
            }
        }
    }
    
    
    int maxf=0,sum=0;  //存最长值,和方案数
    for(int i=1;i<=scc_cnt;i++)
    {
        if(f[i]>maxf)
        {
            maxf=f[i];
            sum=g[i];
        }
        else if(f[i]==maxf)
        {
            sum=(sum+g[i])%mode;
        }
    }
    
    printf("%d\n",maxf);
    printf("%d\n",sum);
    
    
    return 0;
}

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值