BZOJ 1093: [ZJOI2007]最大半连通子图

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

Sample Output

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;
}