1093: [ZJOI2007]最大半连通子图
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 ≤100000, M ≤1000000;对于100%的数据, X ≤10^8
Output
应包含两行,第一行包含一个整数K。第二行包含整数C Mod X.
Sample Input
6 6 20070603
1 2
2 1
1 3
2 4
5 6
6 4
Sample Output
3
3
【解题报告】
缩点之后发现最长路就是答案。跑的过程中记录一下方案数。
跑最长路的时候需要判重。
代码看得黄学长的。。。
代码如下:
/**************************************************************
Problem: 1093
User: onepointo
Language: C++
Result: Accepted
Time:1872 ms
Memory:39536 kb
****************************************************************/
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define N 100010
int mx,ans;
int ind,cnt,scc,top;
int n,m,X;
int last[N],last2[N];
int dfn[N],low[N],hav[N],belong[N],q[N];
int r[N],f[N],g[N],vis[N];
bool inq[N];
struct edge{int to,next;}e[2000005],ed[2000005];
void insert(int u,int v)
{
e[++cnt].to=v;e[cnt].next=last[u];last[u]=cnt;
}
void ins(int u,int v)
{
ed[++cnt].to=v;ed[cnt].next=last2[u];last2[u]=cnt;
r[v]++;
}
void tarjan(int x)
{
dfn[x]=low[x]=++ind;
q[++top]=x;inq[x]=1;
for(int i=last[x];i;i=e[i].next)
if(!dfn[e[i].to])
tarjan(e[i].to),low[x]=min(low[x],low[e[i].to]);
else if(inq[e[i].to])low[x]=min(low[x],dfn[e[i].to]);
int now=0;
if(low[x]==dfn[x])
{
scc++;
while(now!=x)
{
now=q[top];top--;
inq[now]=0;
hav[scc]++;
belong[now]=scc;
}
}
}
void rebuild()
{
cnt=0;
for(int x=1;x<=n;x++)
{
for(int i=last[x];i;i=e[i].next)
if(belong[x]!=belong[e[i].to])
ins(belong[x],belong[e[i].to]);
}
}
void dp()
{
int head=0,tail=0;
for(int i=1;i<=scc;++i)
{
if(!r[i]) q[tail++]=i;
f[i]=hav[i];g[i]=1;
}
while(head!=tail)
{
int now=q[head];head++;
for(int i=last2[now];i;i=ed[i].next)
{
r[ed[i].to]--;
if(!r[ed[i].to])q[tail++]=ed[i].to;
if(vis[ed[i].to]==now)continue;
if(f[now]+hav[ed[i].to]>f[ed[i].to])
{
f[ed[i].to]=f[now]+hav[ed[i].to];
g[ed[i].to]=g[now];
}
else if(f[now]+hav[ed[i].to]==f[ed[i].to])
g[ed[i].to]=(g[ed[i].to]+g[now])%X;
vis[ed[i].to]=now;
}
}
}
int main()
{
scanf("%d%d%d",&n,&m,&X);
for(int i=1;i<=m;++i)
{
int u,v;
scanf("%d%d",&u,&v);
insert(u,v);
}
for(int i=1;i<=n;++i) if(!dfn[i]) tarjan(i);
rebuild();
dp();
for(int i=1;i<=scc;++i)
{
if(f[i]>mx) mx=f[i],ans=g[i];
else if(f[i]==mx) ans=(ans+g[i])%X;
}
printf("%d\n%d\n",mx,ans);
return 0;
}