[ZJOI2007]最大半连通子图（强联通分量+dp最长链）

1093：[ZJOI2007]最大半连通子图

时间限制: 3000 ms
内存限制: 165888 KB

题目描述

一个有向图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的余数。

输入

第一行包含两个整数N，M，X。N，M分别表示图G的点数与边数，X的意义如上文所述接下来M行，每行两个正整数a, b，表示一条有向边(a, b)。图中的每个点将编号为1,2,3…N，保证输入中同一个(a,b)不会出现两次。

N ≤100000, M ≤1000000；对于100%的数据， X ≤$10^8$

输出

应包含两行，第一行包含一个整数K。第二行包含整数C Mod X.

输入样例

6 6 20070603
1 2
2 1
1 3
2 4
5 6
6 4

输出样例

3
3

解：

这不是我认识的zjoi，怎么这么简单？
就是tarjin缩点然后给点权找最长链。
code:

#include<iostream>
#include<cstdio>
#include<cstring>
#include<stack>
using namespace std;
struct lxy{
    int to,next;
}eg[2000005];

int n,m,x,a1,a2,tim,cnt_,cnt,ansl,ansnum;
int head[200005],ti[200005],min_ti[200005];
int bl[200005],ineg[200005],num[200005],lg[200005],f[200005];
bool vis[200005];
stack <int> d;
stack <int> p;

void add(int op,int ed)
{
    eg[++cnt].next=head[op];
    head[op]=cnt;
    eg[cnt].to=ed;
}

void dfs(int u)
{
    d.push(u);
    ti[u]=++tim;
    min_ti[u]=tim;
    for(int i=head[u];i!=-1;i=eg[i].next)
      if(bl[eg[i].to]==0)
      {
         if(min_ti[eg[i].to]==0) dfs(eg[i].to);
         min_ti[u]=min(min_ti[u],min_ti[eg[i].to]);
      }
    if(ti[u]==min_ti[u]){
        cnt_++;
        while(d.top()!=u){
            bl[d.top()]=cnt_;
            num[cnt_]++;
            d.pop();
        }
        bl[d.top()]=cnt_;
        num[cnt_]++;
        d.pop();
    }
}

int main()
{
    memset(head,-1,sizeof(head));
    scanf("%d%d%d",&n,&m,&x);
    cnt_=n;
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d",&a1,&a2);
        add(a1,a2);
    }
    for(int i=1;i<=n;i++)
      if(bl[i]==0)
        dfs(i);
    for(int i=1;i<=n;i++)
      for(int j=head[i];j!=-1;j=eg[j].next)
      if(bl[i]!=bl[eg[j].to])
      {
        add(bl[i],bl[eg[j].to]);
        ineg[bl[eg[j].to]]++;
      }
    for(int i=n+1;i<=cnt_;i++)
      if(ineg[i]==0)
        d.push(i),f[i]=1,lg[i]=num[i];
    while(!d.empty()){
        int now=d.top();d.pop();
        for(int i=head[now];i!=-1;i=eg[i].next)
        {
          if(--ineg[eg[i].to]==0)
            d.push(eg[i].to);
          if(vis[eg[i].to]==0)
          {
            p.push(eg[i].to);vis[eg[i].to]=1;
            if(lg[now]+num[eg[i].to]==lg[eg[i].to])
              f[eg[i].to]=(f[eg[i].to]+f[now])%x;
            else if(lg[now]+num[eg[i].to]>lg[eg[i].to]){
              f[eg[i].to]=f[now];
              lg[eg[i].to]=lg[now]+num[eg[i].to];
            }
          }
        }
        while(!p.empty()){
            vis[p.top()]=0;
            p.pop();
        }
    }
    for(int i=n+1;i<=cnt_;i++)
    {
        if(lg[i]>ansl){
            ansl=lg[i];
            ansnum=f[i];
        }
        else if(lg[i]==ansl)
          ansnum=(ansnum+f[i])%x;
    }
    printf("%d\n%d",ansl,ansnum);
}