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

无事闲来水一水Tarjan,结果水挂了。。。。。。

先Tarjan缩点。

第一问是DAG上的最长路径(点权)

第二问是最长路径的方案数。

都是DP。

然而第二问坑了。。。。。。

因为没判重边(很明显重新建图的时候会出现重边啊)

毕竟还是弱啊。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<stack>
using namespace std;
const int N=100000+5;
const int M=1000000+5;
typedef long long ll;
inline int read(){
	int x=0;char ch;
	while(ch<'0'||ch>'9')ch=getchar();
	while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
	return x;
}
struct Edge{
	int to,next;
}e[M],g[M];
int head[N],cnt,last[N];
int val[N],sccno[N],scc_cnt;
int pre[N],lowlink[N],dfs_clock;
int d[N];
ll f[N],mod;
stack<int>s;
int vis[N];
void ins1(int u,int v){
	cnt++;e[cnt].to=v;e[cnt].next=head[u];head[u]=cnt;
}
void ins2(int u,int v){
	cnt++;g[cnt].to=v;g[cnt].next=last[u];last[u]=cnt;
}
void tarjan(int u){
	pre[u]=lowlink[u]=++dfs_clock;
	s.push(u);
	for(int i=head[u];i;i=e[i].next){
		int v=e[i].to;
		if(!pre[v]){
			tarjan(v);
			lowlink[u]=min(lowlink[u],lowlink[v]);
		}else if(!sccno[v])lowlink[u]=min(lowlink[u],pre[v]);
	}
	if(lowlink[u]==pre[u]){
		scc_cnt++;
		while(true){
			int x=s.top();s.pop();
			sccno[x]=scc_cnt;
			val[scc_cnt]++;
			if(x==u)break;
		}
	}
}
int dp(int u){
	if(d[u])return d[u];
	d[u]=val[u];
	for(int i=last[u];i;i=g[i].next)
	d[u]=max(d[u],val[u]+dp(g[i].to));
	return d[u];
}
ll calc(int u){
	if(f[u])return f[u];
	int lim=0;
	for(int i=last[u];i;i=g[i].next)
	lim=max(lim,d[g[i].to]);
	for(int i=last[u];i;i=g[i].next)
	if(d[g[i].to]==lim){
		if(vis[g[i].to]==u)continue;
		vis[g[i].to]=u;
		f[u]=(f[u]+calc(g[i].to))%mod;
	}
	if(!f[u])f[u]=1;
	return f[u];
}
int main(){
	int n,m;n=read();m=read();mod=read();
	int u,v;
	while(m--){
		u=read();v=read();
		ins1(u,v);
	}
	for(int i=1;i<=n;i++)if(!pre[i])tarjan(i);
	cnt=0;
	for(int u=1;u<=n;u++)
	for(int i=head[u];i;i=e[i].next){
		int v=e[i].to;
		if(sccno[u]!=sccno[v])
		ins2(sccno[u],sccno[v]);
	}
	int ans=0;
	ll num=0;
	for(int i=1;i<=scc_cnt;i++)ans=max(ans,dp(i));
	printf("%d\n",ans);
	for(int i=1;i<=scc_cnt;i++)if(d[i]==ans)num=(num+calc(i))%mod;
	printf("%lld",num);
	return 0;
}


你可能感兴趣的:(1093: [ZJOI2007]最大半连通子图)