C-[ZJOI2007]最大半连通子图_2022图论班第二章连通性例题与习题 (nowcoder.com)
题目描述t
一个有向图G=(VE)称为半连通的(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的余数。
输入描述:
第一行包含两个整数N,M,X。N,M分别表示图c的点数与边数,x的意义如上文所述接下来M行,每行两个正整数a,b,表示一条有向边(a,b)。
图中的每个点将编号为1,2,3...N,保证输入中同一个(a,b)不会出现两次。N≤100000,M≤ 1000000;对于100%的数据,x≤10^8
输出描述:
应包含两行,第一行包含一个整数K。第二行包含整数c Mod x.
示例1
输入
复制
6 6 20070603
1 2
2 1
1 3
2 4
5 6
6 4
输出
3
3
题解:
根据题意,我们发现在一个强联通分量内,点可以互相到达,肯定是满足题意的,所以我们把他考虑为一个点,用tarjin算法找到所有的强联通分量
其次我们可以发现他要我们找的图是一条链(强联通分量看作点),所以我们应该找一条最长的链,并找到有多少这样的链
进行缩点后,利用拓扑排序,我们可以找到最长的链的长度,同时进行DP
DP的思想如下:
f[i]i之前链的最大长度
ff[i]有多少种方式可以达到f[i]
如果f[u] + 此时点(不是真的点,强连通分量的大小)的数目 > f[v]
f[v] = f[u] + siz[v]
ff[v] = ff[u]
如果f[u] + 此时点(不是真的点,强连通分量的大小)的数目 = f[v]
ff[v] = ff[u] + ff[v]
最终把强连通分量中,最长链数目 = ma 的ff[i]相加即可
#include<iostream>
#include<algorithm>
#include<string>
#include<queue>
#include<vector>
#include<map>
#include<cstring>
#include<cmath>
#include<stack>
#include<set>
using namespace std;
#define int long long
typedef pair<int,int> PII;
int n,m,x;
const int N = 2e5 + 10;
vector<int> p[N];
int dfn[N];
int low[N];
int vis[N];
int idx;
stack<int> s;
int cnt;
int siz[N];
int color[N];
vector<int> ne[N];
int rd[N];
int f[N];
int ff[N];
int pre[N];
int ma;
int cd[N];
queue<int> q;
void tarjin(int u)
{
dfn[u] = low[u] = ++idx;
vis[u] = 1;
s.push(u);
for(auto v:p[u])
{
if(!dfn[v])
{
tarjin(v);
low[u] = min(low[u],low[v]);
}
else if(vis[v])
{
low[u] = min(low[u],dfn[v]);
}
}
if(low[u] == dfn[u])
{
cnt++;
while(1)
{
int t = s.top();
s.pop();
color[t] = cnt;
siz[cnt] ++;
vis[t] = 0;
if(t == u)
break;
}
}
}
void toupu()
{
while(q.size())
{
int u = q.front();
q.pop();
ma = max(ma,f[u]);
for(auto v:ne[u])
{
rd[v]--;
if(rd[v] == 0)
{
q.push(v);
}
if(pre[v] == u)
{
continue;
}
if(f[u] + siz[v] > f[v])
{
f[v] = f[u] +siz[v];
ff[v] = ff[u];
}
else if(f[u] + siz[v] == f[v])
{
ff[v] = (ff[u] + ff[v])%x;
}
pre[v] = u;
}
}
}
void solve()
{
cin >> n >> m >> x;
for(int i = 1;i <= m;i++)
{
int l,r;
cin >> l >> r;
p[l].push_back(r);
}
for(int i = 1;i <= n;i++)
{
if(!dfn[i])
tarjin(i);
}
for(int i = 1;i <= n;i++)
{
for(auto j : p[i])
{
int u = color[i];
int v = color[j];
if(u!=v)
{
ne[u].push_back(v);
cd[u]++;
rd[v]++;
}
}
}
for(int i = 1;i <= cnt;i++)
{
if(!rd[i])
{
ff[i] = 1;
f[i] = siz[i];
q.push(i);
}
}
toupu();
int ans = 0;
for(int i = 1;i <= cnt;i++)
{
if(f[i] == ma)
{
ans = (ans + ff[i])%x;
}
}
cout <<ma<<"\n"<<ans;
}
//1 -7 1 -7
signed main(){
// ios::sync_with_stdio(false);
// cin.tie(0);
// cout.tie(0);
int t = 1;
// cin >> t;
while(t--)
{
solve();
}
}
//3 3 2 2 1 0
//1 2 1 1 1 0