Birdwatching 【Gym - 102501K】

题目链接

---

  抗疫期间，在家读如此长的题目容易烦躁hh，于是我就帮大伙读了。

  有N个点，M条边的无向图，我们给出图P是图G的一个衍生图，图G中的点和边图P中都有，但是图P中可能存在一些多余边，怎么说呢，就是图G中有a->b->c这样的边，那么图P中可能就会新的生成a->c这样的边了，类推有效。

  现在，我们要找一个点T，要确定这样的点T'，使得T'确保在图G中有T‘->T的边。输出T‘点的数量，以及所有的T’（升序排列）。

  这问题可以看成，从起点T到其余邻接点，只能有一次到达的方案，于是我们不妨bfs直接搜，因为到达一个点一次以上都是不可的，但是有一种情况，我们得排除掉。

这种情况，到达右上角那个点是有两种方法，但是都是以它为核心的，所以特殊处理一下，用个vector存一下状态即可。

5 6 0
1 0
2 1
3 2
2 3
0 3
3 0

#include <iostream>
#include <cstdio>
#include <cmath>
#include <string>
#include <cstring>
#include <algorithm>
#include <limits>
#include <vector>
#include <stack>
#include <queue>
#include <set>
#include <map>
#include <bitset>
//#include <unordered_map>
//#include <unordered_set>
#define lowbit(x) ( x&(-x) )
#define pi 3.141592653589793
#define e 2.718281828459045
#define INF 0x3f3f3f3f
#define eps 1e-8
#define HalF (l + r)>>1
#define lsn rt<<1
#define rsn rt<<1|1
#define Lson lsn, l, mid
#define Rson rsn, mid+1, r
#define QL Lson, ql, qr
#define QR Rson, ql, qr
#define myself rt, l, r
#define MP(a, b) make_pair(a, b)
using namespace std;
typedef unsigned long long ull;
typedef unsigned int uit;
typedef long long ll;
const int maxN = 1e5 + 7;
int N, M, T, head[maxN], cnt, ans, p[maxN];
struct Eddge
{
    int nex, to;
    Eddge(int a=-1, int b=0):nex(a), to(b) {}
}edge[maxN];
inline void addEddge(int u, int v)
{
    edge[cnt] = Eddge(head[u], v);
    head[u] = cnt++;
}
int tim[maxN] = {0};
vector<int> vt[maxN];
queue<pair<int, int>> Q;
pair<int, int> now;
inline void bfs()
{
    for(int i=head[T], v; ~i; i=edge[i].nex)
    {
        v = edge[i].to;
        tim[v]++;
        vt[v].push_back(v);
        Q.push(MP(v, v));
    }
    while(!Q.empty())
    {
        now = Q.front(); Q.pop();
        int u = now.first;
        for(int i=head[u], v; ~i; i=edge[i].nex)
        {
            v = edge[i].to;
            if(v == T) continue;
            if(tim[v] < 2)
            {
                if(!tim[v])
                {
                    tim[v]++;
                    vt[v].push_back(now.second);
                    Q.push(MP(v, now.second));
                }
                else if(tim[v] == 1)
                {
                    if(vt[v][0] == now.second) continue;
                    tim[v]++;
                    vt[v].push_back(now.second);
                    Q.push(MP(v, now.second));
                }
            }
        }
    }
}
inline void init()
{
    cnt = ans = 0;
    for(int i=0; i<=N; i++) head[i] = -1;
}
int main()
{
    scanf("%d%d%d", &N, &M, &T);
    init();
    for(int i=1, u, v; i<=M; i++)
    {
        scanf("%d%d", &u, &v);
        addEddge(v, u);
    }
    tim[T] = 1;
    bfs();
    for(int i=head[T], v; ~i; i=edge[i].nex)
    {
        v = edge[i].to;
        if(tim[v] < 2) p[++ans] = v;
    }
    sort(p + 1, p + ans + 1);
    printf("%d\n", ans);
    for(int i=1; i<=ans; i++) printf("%d\n", p[i]);
    return 0;
}