【bzoj1093】最大半连通子图 强连通分量&记忆化搜索

       由于要使点最多，所以同一个强连通分量中的点要么都选，要么都不选。然后tarjan缩点，新点的权值为该强连通分量中点的个数。然后在新的DAG上求一条最长的链并统计最长链的个数即可。记忆化搜索解决（当然也可以bfs顺推然而我懒-_-。

下附AC代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 100005
#define M 1000005
using namespace std;
 
int n,m,mod,tot,cnt,fst[N],pnt[M],nxt[M];
int dfsclk,tp,stk[N],pos[N],low[N],scc[N],sum[N],f[N],g[N];
bool bo[N];
struct edg{ int x,y; }a[M];
bool cmp(edg aa,edg bb){
    return aa.x<bb.x || (aa.x==bb.x && aa.y<bb.y);
}
int read(){
    int x=0; char ch=getchar();
    while (ch<'0' || ch>'9') ch=getchar();
    while (ch>='0' && ch<='9'){ x=x*10+ch-'0'; ch=getchar(); }
    return x;
}
void add(int aa,int bb){
    pnt[++tot]=bb; nxt[tot]=fst[aa]; fst[aa]=tot;
}
void dfs(int x){
    pos[x]=low[x]=++dfsclk; int p;
    stk[++tp]=x;
    for (p=fst[x]; p; p=nxt[p]){
        int y=pnt[p];
        if (!pos[y]){
            dfs(y); low[x]=min(low[x],low[y]);
        } else if (!scc[y]) low[x]=min(low[x],pos[y]);
    }
    if (low[x]==pos[x]){
        cnt++; scc[x]=cnt;
        for (sum[cnt]=1; stk[tp]!=x; tp--){
            scc[stk[tp]]=cnt; sum[cnt]++;
        }
        tp--;
    }
}
void dp(int x){
    if (f[x]) return; int p;
    for (p=fst[x]; p; p=nxt[p]){
        int y=pnt[p]; dp(y);
        if (f[y]>f[x]){ f[x]=f[y]; g[x]=g[y]; } else
        if (f[y]==f[x]) g[x]=(g[x]+g[y])%mod;
    }
    if (!f[x]) g[x]=1; f[x]+=sum[x];
}
int main(){
    n=read(); m=read(); mod=read(); int i,p;
    for (i=1; i<=m; i++){
        int x=read(),y=read(); add(x,y);
    }
    for (i=1; i<=n; i++) if (!pos[i]) dfs(i); 
    m=0;
    for (i=1; i<=n; i++)
        for (p=fst[i]; p; p=nxt[p]){
            int j=pnt[p];
            if (scc[i]!=scc[j]){
                m++; a[m].x=scc[i]; a[m].y=scc[j];
            }
        }
    sort(a+1,a+m+1,cmp);
    tot=0; memset(fst,0,sizeof(fst));
    for (i=1; i<=m; i++)
        if (i==1 || a[i].x!=a[i-1].x || a[i].y!=a[i-1].y){
            add(a[i].x,a[i].y); bo[a[i].y]=1;
        }
    for (i=1; i<=cnt; i++) if (!bo[i]) add(0,i);
    dp(0); printf("%d\n%d\n",f[0],g[0]);
    return 0;
}

by lych

2015.12.4