当前位置： 移动技术网 > IT编程>开发语言>C/C++ > BZOJ1093: [ZJOI2007]最大半连通子图(tarjan dp)

BZOJ1093: [ZJOI2007]最大半连通子图(tarjan dp)

2018年09月03日  | 移动技术网IT编程  | 我要评论

只卖来客不外送,冰寒祛痘教程,我的法宝是电脑

题意

一个有向图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;
}
/*
*/

您可能感兴趣的文章:

如对本文有疑问，请在下面进行留言讨论，广大热心网友会与你互动！！ 点击进行留言回复

相关文章:

验证码：
最近更新的文章
大家感兴趣的文章
移动技术网