P4645 [COCI2006-2007 Contest#3] BICIKLI

题目链接

题意分析

首先\(1→2\)的每一天路径上不可以存在环

所以我们先来一个\(Tarjan\)找环

然后如何确定\(1\)到\(2\)的路径上存在环呢

我们先正向边\(1→2\)跑一次

然后反向边\(2→1\)跑一次

全部都标记到的点就是会经过的点

如果有一个点位于环上 就是\(inf\)

否则的话 正常的拓扑\(dp\)就可以了

CODE:

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<cstdlib>
#include<string>
#include<queue>
#include<map>
#include<stack>
#include<list>
#include<set>
#include<deque>
#include<vector>
#include<ctime>
#define ll long long
#define inf 0x7fffffff
#define N 500008
#define IL inline
#define M 308611
#define D double
#define mod 1000000000
#define R register
using namespace std;
template<typename T>IL void read(T &_)
{
    T __=0,___=1;char ____=getchar();
    while(!isdigit(____)) {if(____=='-') ___=0;____=getchar();}
    while(isdigit(____)) {__=(__<<1)+(__<<3)+____-'0';____=getchar();}
    _=___ ? __:-__;
}
/*-------------OI使我快乐-------------*/
int n,m,tot,cnt,top;
int to[M],nex[M],head[M],sta[M];
int dfn[M],low[M],siz[M],ans[M],bel[M],in[M];
queue<int> Q;
bool vis[M],vis_cdy[M],vis_wzy[M];
vector<int> cdy[M],wzy[M];
IL void add(int x,int y)
{to[++tot]=y;nex[tot]=head[x];head[x]=tot;}
IL void dfs(int now)
{
    dfn[now]=low[now]=++cnt;vis[now]=1;
    sta[++top]=now;
    for(R int i=head[now];i;i=nex[i])
    {
        int v=to[i];
        if(!dfn[v]) {dfs(v);low[now]=min(low[now],low[v]);}
        else if(vis[v]) low[now]=min(low[now],dfn[v]);
    }
    if(dfn[now]==low[now])
    {
        ++tot;
        while(sta[top+1]!=now)
        {
            bel[sta[top]]=tot;
            siz[tot]++;
            vis[sta[top--]]=0;
        }
    }
}
IL void bfs_cdy()
{
    memset(vis_cdy,0,sizeof vis_cdy);
    Q.push(1);vis_cdy[1]=1;
    for(;!Q.empty();)
    {
        int u=Q.front();Q.pop();
        for(R int i=0;i<(int)cdy[u].size();++i)
        {
//          printf("%d -> %d\n",u,cdy[u][i]);
            in[cdy[u][i]]++;
            if(!vis_cdy[cdy[u][i]])
            {
                vis_cdy[cdy[u][i]]=1;
                
                Q.push(cdy[u][i]);
            }
        }
    }
}
IL void bfs_wzy()
{
    memset(vis_wzy,0,sizeof vis_wzy);
    Q.push(2);vis_wzy[2]=1;
    for(;!Q.empty();)
    {
        int u=Q.front();Q.pop();
        for(R int i=0;i<(int)wzy[u].size();++i)
        {
            if(!vis_wzy[wzy[u][i]])
            {
                vis_wzy[wzy[u][i]]=1;
                Q.push(wzy[u][i]);
            }
        }
    }
}
IL void Bfs_ans()
{
    Q.push(1);ans[1]=1;
    for(;!Q.empty();)
    {
        int u=Q.front();Q.pop();
        for(R int i=head[u];i;i=nex[i])
        {
            int v=to[i];
            if(vis_wzy[v])
            {
                --in[v];
                ans[v]=(ans[v]+ans[u])%mod;
                if(!in[v]) Q.push(v);
            }
        }
    }
}
int main()
{
//  freopen(".in","r",stdin);
//  freopen(".out","w",stdout);
    read(n);read(m);
    for(R int i=1,x,y;i<=m;++i)
    {
        read(x);read(y);add(x,y);
        cdy[x].push_back(y);
        wzy[y].push_back(x);
    }
    tot=0;
    for(R int i=n;i;--i) if(!dfn[i]) dfs(i);
//  for(R int i=1;i<=n;++i) printf("%d%c",bel[i],(i==n ? '\n':' '));
    bfs_cdy();bfs_wzy();
    for(R int i=1;i<=n;++i)
    if(vis_cdy[i]&&vis_wzy[i]&&siz[bel[i]]>1)
    {puts("inf");return 0;}
//  for(R int i=1;i<=n;++i) printf("%d%c",in[i],(i==n ? '\n':' '));
    Bfs_ans();
    printf("%d\n",ans[2]);
//  fclose(stdin);
//  fclose(stdout);
    return 0;
}

HEOI 2019 RP++

转载于:https://www.cnblogs.com/LovToLZX/p/10628476.html