最大半连通子图 bzoj 1093
最大半连通子图
【问题描述】
一个有向图G = (V,E)称为半连通的(Semi-Connected),如果满足:∀ u, v ∈V,满足u—>v 或 v —> u,即对于图中任意两点u,v, 存在一条u到v的有向路径或者从v到u的有向路径。
若满足,则称G’是G的一个导出子图。
若G’是G的导出子图,且G’半连通,则称G’为G的半连通子图。
若G’是G所有半连通子图中包含节点数最多的,则称G’是G的最大半连通子图。
给定一个有向图G,请求出G的最大半连通子图拥有的节点数K,以及不同的最大半连通子图的数目C。由于C可能比较大,仅要求输出C对X的余数。
【输入文件】
第一行包含两个整数N,M,X。N,M分别表示图G的点数与边数,X的意义如上文所述。接下来M行,每行两个正整数a, b,表示一条有向边(a, b)。图中的每个点将编号为1,2,3…N,保证输入中同一个(a,b)不会出现两次。
【输出文件】
应包含两行,第一行包含一个整数K。第二行包含整数C Mod X.
【样例输入】
6 6 20070603
1 2
2 1
1 3
2 4
5 6
6 4
【样例输出】
3
3
【数据规模】
对于20%的数据, N ≤ 18;
对于60%的数据, N ≤ 10000;
对于100%的数据, N ≤ 100000, M ≤ 1000000;
对于100%的数据, X ≤ 108。
题解:
首先用Tarjon缩点,去重连边,得到新图,那么题目就变成了求图中最长链及最长链个数
最长链可以直接用拓扑排序
最长链个数用一个类似递推的方法
记录每一个点的方案数
那么当前点的方案数就等于连到此点且满足距离相等的点的方案数之和
最后查找距离等于最长链的点,答案为它们的方案数之和
#include<algorithm>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<cstdio>
#include<cmath>
using namespace std;
inline int Get()
{
int x = ;
char c = getchar();
while('' > c || c > '') c = getchar();
while('' <= c && c <= '')
{
x = (x << ) + (x << ) + c - '';
c = getchar();
}
return x;
}
const int me = ;
int n, m, mo;
int x[me], y[me];
int tot;
int de[me], to[me], fir[me], nex[me];
int ue[me];
int si[me];
inline void Ins(int x, int y)
{
nex[++tot] = fir[x];
fir[x] = tot;
to[tot] = y;
}
int num, top, col;
int ti[me], lo[me], st[me], co[me];
inline void Tarjan(int u)
{
ti[u] = lo[u] = ++num;
st[++top] = u;
for(int i = fir[u]; i; i = nex[i])
{
int v = to[i];
if(!ti[v])
{
Tarjan(v);
lo[u] = min(lo[u], lo[v]);
}
else
if(!co[v])
lo[u] = min(lo[u], ti[v]);
}
if(lo[u] == ti[u])
{
co[u] = ++col;
++si[col];
while(st[top] != u)
{
++si[col];
co[st[top]] = col;
--top;
}
--top;
}
}
int t, w;
int ans;
int e[me];
int dis[me];
inline bool rule(int a, int b)
{
if(x[a] != x[b]) return x[a] < x[b];
return y[a] < y[b];
}
int nu[me];
inline void Remove()
{
for(int i = ; i <= m; ++i)
{
nu[i] = i;
x[i] = co[x[i]];
y[i] = co[y[i]];
}
sort(nu + , nu + + m, rule);
}
inline void Build()
{
tot = ;
memset(fir, , sizeof(fir));
for(int i = ; i <= m; ++i)
{
int z = nu[i];
if((x[z] != y[z]) && (x[z] != x[nu[i - ]] || y[z] != y[nu[i - ]]))
{
++de[y[z]];
Ins(x[z], y[z]);
}
}
}
inline void Reset()
{
for(int i = ; i <= col; ++i)
if(!de[i])
{
ue[++w] = i;
dis[i] = si[i];
e[i] = ;
if(dis[ans] < dis[i]) ans = i;
}
}
inline void Topo()
{
while(t < w)
{
int u = ue[++t];
for(int i = fir[u]; i; i = nex[i])
{
int v = to[i];
--de[v];
if(dis[v] < dis[u] + si[v])
{
dis[v] = dis[u] + si[v];
e[v] = ;
if(dis[ans] < dis[v]) ans = v;
}
if(dis[v] == dis[u] + si[v])
e[v] = (e[v] + e[u]) % mo;
if(!de[v]) ue[++w] = v;
}
}
}
int anss;
inline void Ask()
{
for(int i = ; i <= n; ++i)
if(dis[i] == dis[ans])
anss = (anss + e[i]) % mo;
}
int main()
{
n = Get(), m = Get(), mo = Get();
for(int i = ; i <= m; ++i)
{
x[i] = Get(), y[i] = Get();
Ins(x[i], y[i]);
}
for(int i = ; i <= n; ++i)
if(!ti[i])
Tarjan(i);
Remove();
Build();
Reset();
Topo();
Ask();
printf("%d\n%d", dis[ans], anss);
}
最大半连通子图 bzoj 1093的更多相关文章
- BZOJ 1093 [ZJOI2007] 最大半连通子图(强联通缩点+DP)
题目大意 题目是图片形式的,就简要说下题意算了 一个有向图 G=(V, E) 称为半连通的(Semi-Connected),如果满足图中任意两点 u v,存在一条从 u 到 v 的路径或者从 v 到 ...
- bzoj 1093 最大半连通子图 - Tarjan - 拓扑排序 - 动态规划
一个有向图G=(V,E)称为半连通的(Semi-Connected),如果满足:?u,v∈V,满足u→v或v→u,即对于图中任意两点u,v,存在一条u到v的有向路径或者从v到u的有向路径.若G'=(V ...
- BZOJ 1093 [ZJOI2007]最大半连通子图
1093: [ZJOI2007]最大半连通子图 Time Limit: 30 Sec Memory Limit: 162 MBSubmit: 1986 Solved: 802[Submit][St ...
- bzoj 1093 [ZJOI2007]最大半连通子图(scc+DP)
1093: [ZJOI2007]最大半连通子图 Time Limit: 30 Sec Memory Limit: 162 MBSubmit: 2286 Solved: 897[Submit][St ...
- BZOJ 1093: [ZJOI2007]最大半连通子图( tarjan + dp )
WA了好多次... 先tarjan缩点, 然后题意就是求DAG上的一条最长链. dp(u) = max{dp(v)} + totu, edge(u,v)存在. totu是scc(u)的结点数. 其实就 ...
- BZOJ 1093 最大半连通子图 题解
1093: [ZJOI2007]最大半连通子图 Time Limit: 30 Sec Memory Limit: 162 MBSubmit: 2767 Solved: 1095[Submit][S ...
- [BZOJ]1093 最大半连通子图(ZJOI2007)
挺有意思的一道图论. Description 一个有向图G=(V,E)称为半连通的(Semi-Connected),如果满足:∀u,v∈V,满足u→v或v→u,即对于图中任意两点u,v,存在一条u到v ...
- 【刷题】BZOJ 1093 [ZJOI2007]最大半连通子图
Description 一个有向图G=(V,E)称为半连通的(Semi-Connected),如果满足:?u,v∈V,满足u→v或v→u,即对于图中任意 两点u,v,存在一条u到v的有向路径或者从v到 ...
- bzoj1093: [ZJOI2007]最大半连通子图 scc缩点+dag上dp
一个有向图G=(V,E)称为半连通的(Semi-Connected),如果满足:?u,v∈V,满足u→v或v→u,即对于图中任意两点u,v,存在一条u到v的有向路径或者从v到u的有向路径.若G'=(V ...
随机推荐
- Flexible 弹性盒子模型之CSS flex-basis 属性
实例 设置第二个弹性盒元素的初始长度为 80 像素: div:nth-of-type(2){flex-basis:80px;} 效果预览 浏览器支持 表格中的数字表示支持该属性的第一个浏览器的版本 ...
- class-dump 反编译私有的库和应用
一.下载并安装class-dump 下载class-dump-3.5.dmg 点击下载 下载完成以后双击.dmg的文件,将里面的class-dump拷贝到/usr/local/bin 设置权限chm ...
- 1.Hibernate简介
1.框架简介: 定义:基于java语言开发的一套ORM框架: 优点:a.方便开发; b.大大减少代码量; c.性能稍高(不能与数据库高手相比,较一般数据库使用者 ...
- JavaScript
2015-08-01 16:20 JavaScript使用时需要注意的地方 1.引入JS的位置:最好的做法是把<script>的标签放到HTML文档的最后.</body>标签之 ...
- 技术笔记:Indy控件发送邮件
工作中有个需求需要发送邮件,因为使用的delphi6,所以自然就选择了indy组件,想想这事挺简单的.实现的过程倒是简单,看着Indy的demo很快就完了,毕竟也不是很复杂的功能. 功能要求: 1.压 ...
- 微软开源全新的文档生成工具DocFX
微软放弃Sandcastle有些年头了,微软最近开源了全新的文档生成工具DocFX,目前支持C#和VB,类似JSDoc或Sphinx,可以从源代码中提取注释生成文档之外,而且还有语法支持你加入其他的文 ...
- Web 项目杂记(一)
1.Tomcat 多实例部署 在Tomcat下多实例部署后,发现如下问题,采用etc/init.d/tomcat start方式无法启动,而需要采用startup.sh.查找原因发现,是因为多实例部署 ...
- Atitit 图像处理 调用opencv 通过java api attilax总结
Atitit 图像处理 调用opencv 通过java api attilax总结 1.1. Opencv java api的支持 opencv2.4.2 就有了对java api的支持1 1. ...
- 一步步构造自己的vue2.0+webpack环境
前面vue2.0和webpack都已经有接触了些(vue.js入门,webpack入门之简单例子跑起来),现在开始学习如何构造自己的vue2.0+webpack环境. 1.首先新建一个目录vue-wk ...
- Git(远程仓库:git@oschina)-V2.0
1.注册git@osc(也就是“码云”) 这里会提示注册密码==push密码,反正一定要记住的东西. 2.安装git 这里要设置个人信息 git config --list //查看git信息 g ...