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;
}