注意了:
最长链与最大半联通子图是等价又不完全等价的。
最大半联通子图的节点数 会 与最长链的节点数相同,
但是边数的话,最大半联通子图还是会加进 两个点之间的所有边
但是本题的主要目的是求 最长链上的节点数,所以不要紧
#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;
}