Description
Input
第一行包含两个整数N,M,X。N,M分别表示图G的点数与边数,X的意义如上文所述。接下来M行,每行两个正整数a, b,表示一条有向边(a, b)。图中的每个点将编号为1,2,3…N,保证输入中同一个(a,b)不会出现两次。
Output
应包含两行,第一行包含一个整数K。第二行包含整数C Mod X.
Sample Input
6 6 20070603
1 2
2 1
1 3
2 4
5 6
6 4
1 2
2 1
1 3
2 4
5 6
6 4
Sample Output
3
3
3
HINT
对于20%的数据, N ≤18;对于60%的数据, N ≤10000;对于100%的数据, N ≤100000, M ≤1000000;对于100%的数据, X ≤10^8。
题解
强连通分量缩点重构+拓扑排序dp首先看题目感觉就跟强连通有点相关。我们可以发现一个强连通分量必定满足半联通的性质,所以,一个半联通子图中可以包含(合法情况下)若干个强连通分量。所以我们可以先用强连通缩点。之后,我们考虑要使一个有向图满足半联通,通过画图可发现,除了强连通分量外,一条链就满足半联通的性质,所以可以用拓扑排序找最长链。Ps:经过这两天的省队集训感觉自己太弱了,代码风格中有按照常数优化的方法改了一些。感觉怪怪的……
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<cmath>
#include<algorithm>
using namespace std;
int n,m,mod,zz,head[100002],hd[100002];
struct bian {int to,nx;} e[1000002],tr[1000002];
int pre[100002],low[100002],sccnum[100002],ct,scc;
int stack[100002],top,have[100002];
int in[100002],q[100002],f[100002],num[100002],vis[100002];
void insert(int x,int y)
{zz++; e[zz].to=y; e[zz].nx=head[x]; head[x]=zz;}
void init()
{
int i,x,y;
scanf("%d%d%d",&n,&m,&mod);
for(i=1;i<=m;i++)
{scanf("%d%d",&x,&y); insert(x,y);}
}
void dfs(int x)
{
int i,p;
pre[x]=low[x]=++ct;
stack[++top]=x;
for(i=head[x];i;i=e[i].nx)
{p=e[i].to;
if(!pre[p]) {dfs(p); low[x]=min(low[x],low[p]);}
else if(!sccnum[p]) low[x]=min(low[x],pre[p]);
}
p=-1;
if(pre[x]==low[x])
{scc++;
while(p!=x)
{p=stack[top]; top--;
sccnum[p]=scc;
have[scc]++;
}
}
}
void tarjan()
{
int i;
for(i=1;i<=n;i++)
{if(!pre[i]) dfs(i);}
}
void rebuild()
{
int i,j,p;
zz=0;
for(i=1;i<=n;i++)
for(j=head[i];j;j=e[j].nx)
{p=e[j].to;
if(sccnum[p]!=sccnum[i])
{tr[++zz].to=sccnum[p]; tr[zz].nx=hd[sccnum[i]]; hd[sccnum[i]]=zz;
in[sccnum[p]]++;
}
}
}
void dp()//用拓扑排序找到最长链
{
int i,j,x,p,t=0,w=1;
for(i=1;i<=scc;i++)
{if(in[i]==0)
{q[w++]=i;
f[i]=have[i]; num[i]=1;
}
}
while(t<w)
{x=q[t]; t++;
for(i=hd[x];i;i=tr[i].nx)
{p=tr[i].to;
in[p]--;
if(in[p]==0) q[w++]=p;
if(vis[p]==x) continue;//重边
if(f[p]<f[x]+have[p])
{f[p]=f[x]+have[p];
num[p]=num[x];
}
else if(f[p]==f[x]+have[p])
{num[p]=(num[p]+num[x]);
if(num[p]>=mod) num[p]-=mod;
}
vis[p]=x;
}
}
}
int main()
{
init(); tarjan(); rebuild(); dp();
int i,ans=0,ann=0;
for(i=1;i<=scc;i++)
{if(f[i]>ans)
{ans=f[i]; ann=num[i];}
else if(f[i]==ans)
{ann=ann+num[i];
if(ann>=mod) ann-=mod;
}
}
printf("%d\n%d\n",ans,ann);
//system("pause");
return 0;
}