Birdwatching

题目链接

题目大意

给一个有向图，和一个点T，
简单点说就是 问跟T相连（x->T）的那些点中，有哪些点是到T必须经过与T相连的那条路。

题解

刚开始的想法：建反边，如果T到那些点的路径中有大于1的就肯定不可以，但是有环这种特殊情况。不知道该怎么判断。
题解：建反边，在直接跟t相连的这些点的集合中如果可以互相到达，那么就不可以。
所以就可以是在这些点的集合中，到达一个点的方式有两次就不可以。但是怎么特判一下环那种的呢。用一个vector，因为只用两次就可以判断出来。
所以直接搜就好了。

#include <algorithm>
#include <cstdio>
#include <iostream>
#include <vector>
#include <stack>
#include <queue>
#include <map>
#include <cmath>
#include <set>
#include <cstring>
#include <string>
#include <bitset>
#include <stdlib.h>
#include <time.h> 
using namespace std;
typedef long long ll;
typedef pair<int,int> pii;
typedef unsigned long long ull;
#define st first
#define sd second
#define mkp make_pair
#define pb push_back
void wenjian(){
   freopen("concatenation.in","r",stdin);freopen("concatenation.out","w",stdout);}
void tempwj(){
   freopen("hash.in","r",stdin);freopen("hash.out","w",stdout);}
ll gcd(ll a,ll b){
   return b == 0 ? a : gcd(b,a % b);}
ll qpow(ll a,ll b,ll mod){
   a %= mod;ll ans = 1;while(b){
   if(b & 1)ans = ans * a % mod;a = a * a % mod;b >>= 1;}return ans;}
struct cmp{
   bool operator()(const pii & a, const pii & b){
   return a.second > b.second;}};
int lb(int x){
   return  x & -x;}

const int inf = 0x3f3f3f3f;
const ll INF = 0x3f3f3f3f3f3f3f3f;
const ll mod = 1e9+7;
const int maxn = 2e5+5;

std::vector<int> v;
std::vector<int> vv[maxn];
std::vector<int> vp[maxn];
std::vector<int> ans;
queue<pii> qq;
int main()
{
   
	int n,m,s;
	scanf("%d%d%d",&n,&m,&s);
	for (int i = 1; i <= m; i ++ )
	{
   
		int x,y;
		scanf("%d%d",&x,&y);
		if(y == s)
			v.pb(x);
		vv[y].pb(x);
	}
	for (int i =0 ; i < v.size();i ++ )
	{
   
		qq.push(mkp(v[i],v[i]));
		vp[v[i]].pb(v[i]);
	}
	while(!qq.empty())
	{
   
		pii t = qq.front();
		qq.pop();
		int x = t.st;
		// printf("%d\n",x);
		for (int i = 0; i < vv[x].size(); i ++ )
		{
   
			int v = vv[x][i];
			if(v == s)
				continue;
			if(vp[v].size() == 0)
			{
   
				qq.push(mkp(v,t.sd));
				vp[v].pb( t.sd);
				// printf("%d %d\n",v,t.sd);
			}
			else if(vp[v].size() == 1)
			{
   
				if(vp[v][0] != t.sd)
				{
   
					qq.push(mkp(v,t.sd));
					vp[v].pb(t.sd);
					// printf("%d %d\n",v,t.sd);
				}
			}
		}
	}
	for (int i = 0; i < v.size(); i ++ )
	{
   
		if(vp[v[i]].size() < 2)
		{
   
			ans.pb(v[i]);
		}
	}
	printf("%d\n",ans.size());
	for (int i =0 ; i < ans.size(); i ++ )
	{
   
		printf("%d\n",ans[i]);
	}
	// printf("\n");
}