码迷,mamicode.com
首页 > 其他好文 > 详细

【BZOJ】1093: [ZJOI2007]最大半连通子图(tarjan+拓扑序)

时间:2014-11-19 14:07:39      阅读:229      评论:0      收藏:0      [点我收藏+]

标签:des   blog   http   io   ar   os   sp   for   数据   

http://www.lydsy.com/JudgeOnline/problem.php?id=1093

两个条件综合起来加上求最大的节点数,那么很明显如果是环一定要缩点。

然后再仔细思考下就是求dag的最长路的数目啦。。。

然后wa了。。。

看了题解。。。噗!第一次注意到缩点后会有重边QAQ。。。于是。。

orz orz

然后思考了下怎么处理重边。。。很简单,每个点bfs时记录一下就行了。。

#include <cstdio>
#include <cstring>
#include <cmath>
#include <string>
#include <iostream>
#include <algorithm>
#include <queue>
#include <set>
#include <map>
using namespace std;
typedef long long ll;
#define pii pair<int, int>
#define mkpii make_pair<int, int>
#define pdi pair<double, int>
#define mkpdi make_pair<double, int>
#define pli pair<ll, int>
#define mkpli make_pair<ll, int>
#define rep(i, n) for(int i=0; i<(n); ++i)
#define for1(i,a,n) for(int i=(a);i<=(n);++i)
#define for2(i,a,n) for(int i=(a);i<(n);++i)
#define for3(i,a,n) for(int i=(a);i>=(n);--i)
#define for4(i,a,n) for(int i=(a);i>(n);--i)
#define CC(i,a) memset(i,a,sizeof(i))
#define read(a) a=getint()
#define print(a) printf("%d", a)
#define dbg(x) cout << (#x) << " = " << (x) << endl
#define error(x) (!(x)?puts("error"):0)
#define printarr2(a, b, c) for1(_, 1, b) { for1(__, 1, c) cout << a[_][__]; cout << endl; }
#define printarr1(a, b) for1(_, 1, b) cout << a[_] << ‘\t‘; cout << endl
inline const int getint() { int r=0, k=1; char c=getchar(); for(; c<‘0‘||c>‘9‘; c=getchar()) if(c==‘-‘) k=-1; for(; c>=‘0‘&&c<=‘9‘; c=getchar()) r=r*10+c-‘0‘; return k*r; }
inline const int max(const int &a, const int &b) { return a>b?a:b; }
inline const int min(const int &a, const int &b) { return a<b?a:b; }
#define rdm(x,i) for(int i=ihead[x]; i; i=e[i].next)

const int N=100005, M=1000005;
int TT, FF[N], LL[N], vis[N], s[N], top, scc, p[N], w[N], MD, ans1, ans2, in[N], front, tail;
struct Gr {
	int ihead[N], n, cnt;
	struct dat { int next, to; }e[M];
	void add(int u, int v) { e[++cnt].next=ihead[u]; ihead[u]=cnt; e[cnt].to=v; ++in[v]; }
	void dfs(int x) {
		FF[x]=LL[x]=++TT; vis[x]=1; s[++top]=x;
		rdm(x, i) {
			int y=e[i].to;
			if(!FF[y]) {
				dfs(y);
				LL[x]=min(LL[x], LL[y]);
			}
			else if(vis[y]) LL[x]=min(LL[x], FF[y]);
		}
		if(FF[x]==LL[x]) {
			++scc;
			int y;
			do {
				y=s[top--];
				vis[y]=0;
				p[y]=scc;
				w[scc]++;
			} while(x!=y);
		}
	}
	void tarjan() {
		for1(i, 1, n) if(!FF[i]) dfs(i);
	}
	void bfs() {
		int *d=FF, *f=LL, *q=s;
		for1(i, 1, n) d[i]=f[i]=vis[i]=0;
		front=tail=0;
		for1(i, 1, n) if(!in[i]) q[tail++]=i, d[i]=w[i], f[i]=1;
		while(front!=tail) {
			int u=q[front++]; if(front==N) front=0;
			rdm(u, i) {
				int v=e[i].to;
				if(!(--in[v])) { q[tail++]=v; if(tail==N) tail=0; }
				if(vis[v]) continue;
				if(d[v]==d[u]+w[v]) {
					f[v]=(f[v]+f[u])%MD;
				}
				else if(d[v]<d[u]+w[v]) {
					d[v]=d[u]+w[v];
					f[v]=f[u];
				}
				vis[v]=1;
			}
			rdm(u, i) vis[e[i].to]=0;
		}
		ans1=-1;
		for1(i, 1, n) ans1=max(ans1, d[i]);
		for1(i, 1, n) if(ans1==d[i]) ans2=(ans2+f[i])%MD;
	}
}G, g;
void build() {
	for1(i, 1, g.n) in[i]=0;
	G.n=scc;
	for1(x, 1, g.n) for(int i=g.ihead[x]; i; i=g.e[i].next) if(p[x]!=p[g.e[i].to]) G.add(p[x], p[g.e[i].to]);
}

int main() {
	read(g.n);
	int m=getint(); read(MD);
	for1(i, 1, m) { int u=getint(), v=getint(); g.add(u, v); }
	g.tarjan();
	build();
	G.bfs();
	printf("%d\n%d\n", ans1, ans2);
	return 0;
}

  

 


 

 

Description

bubuko.com,布布扣

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

 

对于100%的数据, N ≤100000, M ≤1000000;对于100%的数据, X ≤10^8。

 

Source

 

 

【BZOJ】1093: [ZJOI2007]最大半连通子图(tarjan+拓扑序)

标签:des   blog   http   io   ar   os   sp   for   数据   

原文地址:http://www.cnblogs.com/iwtwiioi/p/4108023.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!