bzoj1093[ZJOI2007]最大半连通子图 tarjan+拓补排序

7 篇文章 0 订阅
6 篇文章 0 订阅

这是一个看懂题目就能切的题。
题意:给你个图,求最长路。
怎么样,是不是被震惊到了。。。其实就这么简单。。
由于原图可能有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);
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值