挺有意思的一道图论。

Description

  一个有向图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的余数。

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

  N ≤10000, M ≤1000000,X ≤10^8。

Solution

  拿到这题,我们首先思考半连通分量是个什么东西。

  首先我们知道,强连通分量一定是半连通分量,

  题目要求我们求最大的半连通分量,所以如果选取了一个强连通分量里的点,那么把该点所在的整个强连通分量都选进去肯定没问题,选取一个强连通分量和选取一个点是等价的。

  所以我们很自然地用tarjan缩了缩点……

  然后我们得到了一个带点权的拓扑图。

  仔细一想,我们发现拓扑图中的半连通分量是一条链,

  所以问题也就变成了找拓扑图中的最长链,并统计最长链的条数。

  (这个用DP不会做你退群吧)

  注意缩点之后要处理掉重边。时间复杂度O(n+m)。

#include <cstdio>
#include <algorithm>
#include <cstring>
#define MN 100005
#define MM 1000005
using namespace std;
struct edge{int nex,to;}e[MM];
struct bian{int x,y;}b[MM];
bool u[MN],ink[MN];
int low[MN],st[MN],hr[MN],bel[MN],w[MN],d[MN],q[MN],f1[MN],f2[MN];
int dfn,pin,tp,n,m,mod,hd,tl,ans1,ans2;
char B[<<],*SS,*TT; inline char getc() {return SS==TT&&(TT=(SS=B)+fread(B,,<<,stdin),SS==TT)?EOF:*SS++;}
inline int read()
{
register int n=; char c;
do c=getc(); while (c<'' || c>'');
do n=n*+c-'',c=getc(); while (c>='' && c<='');
return n;
} inline void ins(int x,int y) {e[++pin]=(edge){hr[x],y}; hr[x]=pin;} void tarjan(int x)
{
register int i,lt;
low[x]=lt=++dfn;
u[x]=ink[x]=true; st[++tp]=x;
for (i=hr[x];i;i=e[i].nex)
{
if (u[e[i].to]&&!ink[e[i].to]) continue;
if (!u[e[i].to]) tarjan(e[i].to);
low[x]=min(low[x],low[e[i].to]);
}
if (low[x]==lt)
for (;st[tp+]!=x;--tp) ink[st[tp]]=false,bel[st[tp]]=x,++w[x];
} int main()
{
register int i,x;
n=read(); m=read(); mod=read();
for (i=;i<=m;++i) b[i].x=read(),b[i].y=read(),ins(b[i].x,b[i].y);
for (i=;i<=n;++i) if (!u[i]) tarjan(i);
memset(hr,,sizeof(hr)); pin=;
for (i=;i<=m;++i)
if (bel[b[i].x]!=bel[b[i].y])
ins(bel[b[i].x],bel[b[i].y]),++d[bel[b[i].y]];
for (i=hd=;i<=n;++i) if (bel[i]==i&&!d[i]) q[++tl]=i,f1[i]=w[i],f2[i]=;
for (;hd<=tl;++hd)
{
for (x=q[hd],i=hr[x];i;i=e[i].nex)
{
if (!--d[e[i].to]) q[++tl]=e[i].to;
if (!u[e[i].to]) continue; else u[e[i].to]=false;
if (f1[x]+w[e[i].to]>f1[e[i].to]) f1[e[i].to]=f1[x]+w[e[i].to],f2[e[i].to]=f2[x];
else if (f1[x]+w[e[i].to]==f1[e[i].to]) f2[e[i].to]+=f2[x],f2[e[i].to]-=f2[e[i].to]>=mod?mod:;
}
for (i=hr[x];i;i=e[i].nex) u[e[i].to]=true;
}
for (i=;i<=n;++i)
if (f1[i]>ans1) ans1=f1[i],ans2=f2[i];
else if (f1[i]==ans1) ans2+=f2[i],ans2-=ans2>=mod?mod:;
printf("%d\n%d",ans1,ans2);
}

Last Word

  有向图用tarjan缩完点得到的是拓扑图,无向图缩环会变成树(森林)。

  用了n+e光速读入后立竿见影地卡到了这道题的rank2。(代码画风崩坏不可避)

[BZOJ]1093 最大半连通子图(ZJOI2007)的更多相关文章

  1. bzoj 1093 最大半连通子图 - Tarjan - 拓扑排序 - 动态规划

    一个有向图G=(V,E)称为半连通的(Semi-Connected),如果满足:?u,v∈V,满足u→v或v→u,即对于图中任意两点u,v,存在一条u到v的有向路径或者从v到u的有向路径.若G'=(V ...

  2. BZOJ 1093 最大半连通子图 题解

    1093: [ZJOI2007]最大半连通子图 Time Limit: 30 Sec  Memory Limit: 162 MBSubmit: 2767  Solved: 1095[Submit][S ...

  3. BZOJ 1093 最大半连通子图

    缩点求最长链. #include<iostream> #include<cstdio> #include<cstring> #include<algorith ...

  4. Tarjan水题系列(5):最大半连通子图 [ZJOI2007 luogu P2272]

    题目 大意: 缩点后转为求最长链的长度和最长链的个数 思路: 看懂题就会做系列 长度和个数都可以拓扑排序后DP求得 毕竟是2007年的题 代码: 如下 #include <cstdio> ...

  5. BZOJ 1093 [ZJOI2007] 最大半连通子图(强联通缩点+DP)

    题目大意 题目是图片形式的,就简要说下题意算了 一个有向图 G=(V, E) 称为半连通的(Semi-Connected),如果满足图中任意两点 u v,存在一条从 u 到 v 的路径或者从 v 到 ...

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

    1093: [ZJOI2007]最大半连通子图 Time Limit: 30 Sec  Memory Limit: 162 MBSubmit: 1986  Solved: 802[Submit][St ...

  7. bzoj 1093 [ZJOI2007]最大半连通子图(scc+DP)

    1093: [ZJOI2007]最大半连通子图 Time Limit: 30 Sec  Memory Limit: 162 MBSubmit: 2286  Solved: 897[Submit][St ...

  8. BZOJ 1093: [ZJOI2007]最大半连通子图( tarjan + dp )

    WA了好多次... 先tarjan缩点, 然后题意就是求DAG上的一条最长链. dp(u) = max{dp(v)} + totu, edge(u,v)存在. totu是scc(u)的结点数. 其实就 ...

  9. 【刷题】BZOJ 1093 [ZJOI2007]最大半连通子图

    Description 一个有向图G=(V,E)称为半连通的(Semi-Connected),如果满足:?u,v∈V,满足u→v或v→u,即对于图中任意 两点u,v,存在一条u到v的有向路径或者从v到 ...

随机推荐

  1. GPUImage实战问题解决

    在项目中遇到了使用完GPUImage以后,内存不释放的问题,翻阅官方API,找到了解决方法: deinit{ GPUImageContext.sharedImageProcessingContext( ...

  2. css中的em 简单教程 -- 转

    先附上原作的地址: https://www.w3cplus.com/css/px-to-em 习惯性的复制一遍~~~~ -------------------------------我是分界线---- ...

  3. Python之旅.第三章.函数3.30

    一.迭代器 1.什么是迭代?:迭代是一个重复的过程,并且每次重复都是基于上一次的结果而来2.要想了解迭代器到底是什么?必须先了解一个概念,即什么是可迭代的对象?可迭代的对象:在python中,但凡内置 ...

  4. Python内置函数(40)——dir

    英文文档: dir([object]) Without arguments, return the list of names in the current local scope. With an ...

  5. axure 预览"HTTP/1.1 302 Found"

    使用Axure编辑原型时,点击预览出现"HTTP/1.1 302 Found" 第一想到的就是重新安装Axure和检查原型文件是否损坏,验证后证明前Axure和.rp文件都是完好的 ...

  6. maven的使用之一简单的安装

    首先,我们知道,在传统的项目中,我们会导入一堆的jar包,那样的话,我们会发现我们的jar包的大小已经占了整个项目大小的90%以上,甚至更多,而且,我们的jar包只能自己使用,如果 其他人想用的话,还 ...

  7. 03、NetCore2.0下Web应用之搭建最小框架

    03.NetCore2.0下Web应用之搭建最小框架 这里我们不使用VS2017或者CLI命令的方式创建Asp.Net Core 2.0网页应用程序,而是完全手工的一点点搭建一个Web框架,以便更好的 ...

  8. z-index的权重是叠加的

    <!DOCTYPE html> <html lang="en"> <head> <meta charset="utf-8&quo ...

  9. 前端之旅HTML与CSS篇之清除浮动塌陷

    以下内容为转载. 方法1:给浮动的元素的上级添加高度如果一个元素要浮动,那么它的祖先元素一定要有高度.高度的盒子,才能关住浮动.只要浮动在一个有高度的盒子中,那么这个浮动就不会影响后面的浮动元素.所以 ...

  10. 前端之旅HTML与CSS篇之a便签中放入其他块元素会撑大高度的原因

    原因:a元素下有一个匿名文本,这个文本外有一个匿名行级盒子,它有的默认vertical-align是baseline的,而且往往因为上文line-height的影响,使它有个line-height,从 ...