BZOJ1093 [ZJOI2007]最大半连通子图

首先，我们要tarjan。。。 然后我们要缩点。。。

注意，缩点的时候两个新建的点会有重边，需要判重

正常的判重方法是bfs一边，但是我YY的比较奇葩，方法下面将。。。

缩好点就变成了一个DAG，然后就类似树形DP的方法求最大权值链

我是用记忆化搜索，当dfs某个点p时用数组vis记录一些东西：

首先vis[x] > 0代表x已经访问过了， vis[x] == p表示vis[x]上次是由p来的，于是就可以去重了。。。

然后就没有然后了，搞定yeah！

/**************************************************************
    Problem: 1093
    User: rausen
    Language: C++
    Result: Accepted
    Time:1604 ms
    Memory:31220 kb
****************************************************************/
 
#include <cstdio>
#include <algorithm>
#include <cstring>
 
using namespace std;
const int inf = (int) 1e8;
struct edges{
    int from, to, next;
}e[2000005];
 
int first[100005],tot;
int dfn[100005], low[100005], w[100005], sum[100005];
int vis[100005], s[100005], f[100005], g[100005], F[100005];
int n, m, MOD, X, Y, ans, top, cnt, num, ANS;
 
void add_edge(int x, int y){
    e[++tot].next = first[x];
    first[x] = tot;
    e[tot].from = x, e[tot].to = y;
}
 
void tarjan(int p){
    dfn[p] = low[p] = ++cnt;
    s[++top] = p;
    vis[p] = 1;
    int x, y;
    for (x = first[p]; x; x = e[x].next){
        y = e[x].to;
        if (!vis[y]) tarjan(y);
        if (vis[y] < 2)
            low[p] = min(low[p], low[y]);
    }
    if (dfn[p] == low[p]){
        ++num;
        while (s[top + 1] != p){
            int y = s[top--];
            vis[y] = 2, w[y] = num;
            ++sum[num];
        }
    }
}
 
void dfs(int p){
    if (vis[p]) return;
    vis[p] = inf;
    int x, y;
    for (int x = first[p]; x; x = e[x].next){
        y = e[x].to;
        if (vis[y] == p) continue;
        vis[y] = p;
        dfs(y);
        if (f[p] < f[y])
            f[p] = f[y], g[p] = g[y]; else
        if (f[p] == f[y]) g[p] += g[y];
    }
    g[p] %= MOD;
    f[p] += sum[p];
    if (ans < f[p])
        ans = f[p], ANS = g[p]; else
    if (ans == f[p]) ANS += g[p], ANS %= MOD;
}
 
int main(){
    scanf("%d%d%d", &n, &m, &MOD);
    for (int i = 1; i <= m; ++i)
        scanf("%d%d", &X, &Y), add_edge(X, Y);
    for (int i = 1; i <= n; ++i)
        if (!vis[i]) tarjan(i);
     
    memset(first, 0, sizeof(first));
    tot = 0;
    for (int i = 1; i <= m; ++i){
        X = e[i].from;
        Y = e[i].to;
        if (w[X] != w[Y]) add_edge(w[X], w[Y]);
    }
    memset(vis, 0, sizeof(vis));
    for (int i = 1; i <= num; ++i)
        if (!first[i]) g[i] = 1;
    for (int i = 1; i <= num; ++i)
        if (!vis[i]) dfs(i);
    printf("%d
%d
", ans, ANS % MOD);
    return 0;
}

(p.s. 各种压代码+卡内存此题status终于上了前两页。。。)