bzoj1093[ZJOI2007]最大半連通子圖 tarjan+拓補排序

這是一個看懂題目就能切的題。
題意：給你個圖，求最長路。
怎麼樣，是不是被震驚到了。。。其實就這麼簡單。。
由於原圖可能有scc，我們需要tarjan處理掉，然後就是最長路了，由於並不固定起點，我們拓補排序的同時去更新最長路，也可說是DP吧，記錄一下到達每個點的方法數，這個也很好記錄。
1A，爽。

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<queue>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
const int N=1e6+5;
int mo;
int n,m,head[N],next[N],go[N];
int scc,sz,tot,low[N],dfn[N],pre[N];
int dis[N],sum[N],bl[N],size[N];
int sta[N],top,d[N];
struct node
{
    int x,y;
}e[N];
queue<int> q;
bool vis[N];
inline void add(int x,int y)
{
    go[++tot]=y;
    next[tot]=head[x];
    head[x]=tot;
}
inline void tarjan(int x)
{
    dfn[x]=low[x]=++sz;
    sta[++top]=x;
    vis[x]=1;
    for(int i=head[x];i;i=next[i])
    {
        int v=go[i];
        if (!dfn[v])
        {
            tarjan(v);
            low[x]=min(low[x],low[v]);
        }
        else if (vis[v])low[x]=min(low[x],dfn[v]);
    }
    if (dfn[x]==low[x])
    {
        ++scc;
        int y=0;
        while (y!=x)
        {
            y=sta[top--];
            bl[y]=scc;
            size[scc]++;
            vis[y]=0;
        } 
    }
}
inline void topsort()
{
    memset(dis,0,sizeof(dis));
    fo(i,1,n)
    if (!d[i])
    {
        q.push(i);
        dis[i]=size[i];
        sum[i]=1;
    }
    while (!q.empty())
    {
        int x=q.front();
        q.pop();
        for(int i=head[x];i;i=next[i])
        {
            int v=go[i];
            d[v]--;
            if (!d[v])q.push(v);
            if (x==pre[v])continue;
            pre[v]=x;
            if (dis[v]==dis[x]+size[v])
            (sum[v]+=sum[x])%=mo;
            if (dis[v]<dis[x]+size[v])
            {
                dis[v]=dis[x]+size[v];
                sum[v]=sum[x];
            }
        }
    }
}
int main()
{
    scanf("%d%d%d",&n,&m,&mo);
    fo(i,1,m)
    {
        scanf("%d%d",&e[i].x,&e[i].y);
        add(e[i].x,e[i].y);
    }
    fo(i,1,n)if (!dfn[i])tarjan(i);
    tot=0;
    memset(head,0,sizeof(head));
    fo(i,1,m)
    {
        int x=e[i].x,y=e[i].y;
        if (bl[x]!=bl[y])add(bl[x],bl[y]),d[bl[y]]++;
    }
    //fo(i,1,n)printf("%d ",bl[i]);printf("\n");
    topsort();
    int ans1=0,ans2=0;
    fo(i,1,scc)ans1=max(ans1,dis[i]);
    fo(i,1,scc)
    if (ans1==dis[i])(ans2+=sum[i])%=mo;
    printf("%d\n%d\n",ans1,ans2);
}