BZOJ1093: [ZJOI2007]最大半连通子图

题目链接

【分析】

（题目看了好半天。。但是其实很好懂得。。
缩点+dp找最长链和最长链方案数。还是很好想的。但是需要注意重构图后关于重边的处理。
为什么要去重边呢？？这个问题困扰了我半天。注意原题中一句话即可

【代码】

#include <cstdio>
#include <iostream>
#include <queue>
#include <vector>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <stack>
#define N 100005
#define M 1000005
#define INF 1000000000
using namespace std;
typedef long long ll;
typedef pair<ll,ll> pa;

int read()
{
    int x=0,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-') f=-1;ch=getchar();}
    while(isdigit(ch)){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}

int n,m,ind,scc,cnt,mod,mx,ans;
int b[M],p[N],nextedge[M];
int dfn[N],low[N],Belong[N],num[N],in[N];
int dp[N],f[N];
bool Instack[N];
stack<int>st;
vector<int>g[N];

void Add(int x,int y)
{
    cnt++;
    b[cnt]=y;
    nextedge[cnt]=p[x];
    p[x]=cnt;
}

void Input_Init()
{
    n=read(),m=read(),mod=read();
    for(int i=1;i<=m;i++)
    {
        static int x,y;
        x=read(),y=read();
        Add(x,y);
    }
}

void Tarjan(int x)
{
    dfn[x]=low[x]=++ind;
    st.push(x);Instack[x]=1;
    for(int i=p[x];i;i=nextedge[i])
    {
        int v=b[i];
        if(!dfn[v])
        {
            Tarjan(v);
            low[x]=min(low[x],low[v]);
        }
        else if(Instack[v]) low[x]=min(low[x],dfn[v]);
    }
    if(dfn[x]==low[x])
    {
        scc++;int now=-1;
        while(now!=x)
        {
            now=st.top();st.pop();
            Belong[now]=scc;
            Instack[now]=0;
            num[scc]++;
        }
    }
}

void Get_Scc(){
    for(int i=1;i<=n;i++) if(!dfn[i]) Tarjan(i);
}

void Build_Graph()
{
    for(int i=1;i<=n;i++)
    for(int j=p[i];j;j=nextedge[j])
        if(Belong[i]!=Belong[b[j]])
            g[Belong[i]].push_back(Belong[b[j]]);
    for(int i=1;i<=scc;i++)
    {
        sort(g[i].begin(),g[i].end());
        g[i].erase(unique(g[i].begin(),g[i].end()),g[i].end());
    }
    for(int i=1;i<=scc;i++)
    for(int j=0;j<g[i].size();j++)
        in[g[i][j]]++;
}

void Top_Sort()
{
    queue<int>q;
    for(int i=1;i<=scc;i++) 
    {
        if(!in[i]) q.push(i);
        dp[i]=num[i];f[i]=1;
    }
    while(!q.empty())
    {
        int k=q.front();q.pop();
        for(int i=0;i<g[k].size();i++)
        {
            int v=g[k][i];
            in[v]--;if(!in[v]) q.push(v);
            if(dp[v]==dp[k]+num[v]) f[v]+=f[k];
            else if(dp[v]<dp[k]+num[v]) f[v]=f[k],dp[v]=dp[k]+num[v];
            while(f[v]=>mod) f[v]-=mod;
        }
    }
}

void Solve()
{
    Top_Sort();
    for(int i=1;i<=scc;i++)
    {
        if(mx<dp[i])
        {
            mx=dp[i];
            ans=f[i];
        }
        else if(mx==dp[i]) ans+=f[i];
        while(ans>=mod) ans-=mod;
    }
    printf("%d\n%d\n",mx,ans);
}

int main()
{
    Input_Init();
    Get_Scc();
    Build_Graph();
    Solve();
    return 0;
}