BZOJ1093: [ZJOI2007]最大半连通子图(tarjan dp)
2018-09-01 05:38:25来源:博客园 阅读 ()
题意
一个有向图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的余数。
Sol
很zz的题然而我因为没判重边的缘故wa了好久qwq
首先强连通分量内的点一定是半联通图
如果任意链各个强连通分量之间有边的话,它们构成的图是半联通图
那么我们最长路dp一下就好,同时dp出方案数。
#include<cstdio> #include<algorithm> #include<stack> #include<queue> #define Pair pair<int, int> #define fi first #define se second #define MP(x, y) make_pair(x, y) //#define int long long using namespace std; const int MAXN = 2 * 1e5 + 10; inline int read() { char c = getchar(); int x = 0, f = 1; while(c < '0' || c > '9') {if(c == '-') f = -1; c = getchar();} while(c >= '0' && c <= '9') x = x * 10 + c - '0', c = getchar(); return x * f; } int N, M, mod; vector<int> v[MAXN], E[MAXN]; int low[MAXN], dfn[MAXN], col[MAXN], vis[MAXN], tot, cn, f[MAXN], g[MAXN], inder[MAXN], siz[MAXN], fuck[MAXN]; stack<int> s; void tarjan(int x) { low[x] = dfn[x] = ++tot; vis[x] = 1; s.push(x); for(int i = 0; i < v[x].size(); i++) { int to = v[x][i]; if(!dfn[to]) tarjan(to), low[x] = min(low[x], low[to]); else if(vis[to]) low[x] = min(low[x], dfn[to]); } if(dfn[x] == low[x]) { int h; cn++; do { h = s.top(); s.pop(); col[h] = cn; vis[h] = 0; siz[cn]++; }while(x != h); } } void Topsort() { queue<int> q; for(int i = 1; i <= cn; i++) { if(!inder[i]) q.push(i); f[i] = siz[i], g[i] = 1; } Pair ans = MP(0, 0); while(!q.empty()) { int p = q.front(); q.pop(); for(int i = 0; i < E[p].size(); i++) { int to = E[p][i], val = f[p] + siz[to]; inder[to]--; if(!inder[to]) q.push(to); if(fuck[to] == p) continue; if(f[to] == val) (g[to] += g[p]) %= mod; else if(val > f[to]) f[to] = val, g[to] = g[p] % mod; fuck[to] = p; } } for(int i = 1; i <= cn; i++) if(f[i] > ans.fi) ans.fi = f[i]; for(int i = 1; i <= cn; i++) if(f[i] == ans.fi) ans.se = (ans.se + g[i]) % mod; printf("%d\n%d", ans.fi, ans.se % mod); } main() { // freopen("3.in", "r", stdin); N = read(); M = read(); mod = read(); for(int i = 1; i <= M; i++) { int x = read(), y = read(); v[x].push_back(y); } for(int i = 1; i <= N; i++) if(!dfn[i]) tarjan(i); for(int x = 1; x <= N; x++) { // sort(v[x].begin(), v[x].end(), comp); // int last = -1; for(int i = 0; i < v[x].size(); i++) { int to = v[x][i]; // if(i > 0 && (col[to] == col[last])) continue; if(col[x] != col[to]) E[col[x]].push_back(col[to]), inder[col[to]]++; // last = to; } } Topsort(); return 0; } /* */
标签:
版权申明:本站文章部分自网络,如有侵权,请联系:west999com@outlook.com
特别注意:本站所有转载文章言论不代表本站观点,本站所提供的摄影照片,插画,设计作品,如需使用,请与原作者联系,版权归原作者所有
- BZOJ1059: [ZJOI2007]矩阵游戏(二分图匹配) 2018-07-09
- 【程序人生】 我的程序员之路 2018-06-18
- C#学习笔记-封装 2018-06-18
- P1169 [ZJOI2007]棋盘制作 2018-06-17
- BZOJ 1059: [ZJOI2007]矩阵游戏 2018-06-17
IDC资讯: 主机资讯 注册资讯 托管资讯 vps资讯 网站建设
网站运营: 建站经验 策划盈利 搜索优化 网站推广 免费资源
网络编程: Asp.Net编程 Asp编程 Php编程 Xml编程 Access Mssql Mysql 其它
服务器技术: Web服务器 Ftp服务器 Mail服务器 Dns服务器 安全防护
软件技巧: 其它软件 Word Excel Powerpoint Ghost Vista QQ空间 QQ FlashGet 迅雷
网页制作: FrontPages Dreamweaver Javascript css photoshop fireworks Flash