Codeforces Round #436 (Div. 2) F. Cities Excursions 字典序dfs+Tarjan判環

題目鏈接：Cities Excursions
題意：給你一個有向圖，問任意兩點間的字典序最小路徑（如果存在）上的第k個節點是啥。
解法：提前存儲好每個點能直接到達的點，並排序，因爲這樣遍歷就是字典序的遍歷。提前存儲好q個詢問，對q的f按從小到大排序，每次對f相同的詢問做dfs，起始點是f，提前把起始點是f的詢問用vector保存下來，這樣在做dfs的過程中，一次性把具有共同起點的這部分詢問給得出結果，這樣做節約了時間。初始化res數組爲-1，res作用保存結果，在dfs中，用Tarjan判環，棧stk存放路徑。在遍歷過程，只要能加入到棧的點，從起點必能到達；沒被加入到棧的點，從起點走字典序最小路必不能到達該點（可能是環導致的），res還是-1。能被加入到棧的點，結果取決於從起點到達該點，字典序最小路徑的長度（即當前棧的棧頂指針的值）和該查詢的k大小比較，小於關係，說明k超出路徑長度，res還是-1,；大於關係，說明詢問的是該路徑的第k個點，stk數組對應下標的值賦給res[當前詢問]。

感覺思路理得比較清楚了，還不懂看代碼吧，代碼借鑑別人的。

#include <iostream>
#include <cstdio>
#include <fstream>
#include <algorithm>
#include <cmath>
#include <deque>
#include <vector>
#include <queue>
#include <string>
#include <cstring>
#include <map>
#include <stack>
#include <set>
#define Max(a,b) a>b?a:b
#define Min(a,b) a>b?b:a
#define mem(a,b) memset(a,b,sizeof(a))
using namespace std;
typedef long long ll;
int dir[4][2]= {{1,0},{-1,0},{0,1},{0,-1}};
const double eps = 1e-6;
const double Pi = acos(-1.0);
const int INF=0x3f3f3f3f;

const int maxn = 3010;
const int maxq = 4e5+10;

vector<int> edge[maxn];
struct query{
    int id,f,t,k;
    query(int _id = -1, int _f = 0, int _t = 0, int _k = 0): id(_id), f(_f), t(_t), k(_k){}
}qs[maxq];
vector<query> Q[maxn];
int n,m,q;
int tot,top;
int res[maxq];
int dfn[maxn],low[maxn],instk[maxn],stk[maxn];
bool cmp(const query &a, const query &b){
    return a.f < b.f;
}

void dfs(int u, int ok){
    dfn[u] = ++tot;
    low[u] = INF;
    stk[top++] = u;
    instk[u] = 1;
    if(ok){
        for(int i = 0; i < Q[u].size(); i++){
            if(Q[u][i].k <= top) res[Q[u][i].id] = stk[Q[u][i].k-1];
        }
    }
    for(int i = 0; i < edge[u].size(); i++){
        int v = edge[u][i];
        if(!dfn[v]){
            dfs(v,ok && dfn[u] < low[u]);
            low[u] = min(low[u],low[v]);
        }
        else if(instk[v]){
            low[u] = min(dfn[v],low[u]);
        }
    }
    instk[u] = 0;
    --top;
}

int main(){
    while(~scanf("%d%d%d",&n,&m,&q)){
        int a,b;

        for(int i = 1; i <= n; i++){
            edge[i].clear();
        }

        for(int i = 0; i < m; i++){
            scanf("%d%d",&a,&b);
            edge[a].push_back(b);
        }
        for(int i = 1; i <= n; i++){
            sort(edge[i].begin(),edge[i].end());
        }
        for(int i = 0; i < q; i++){
            int f,t,k;
            scanf("%d%d%d",&f,&t,&k);
            qs[i] = query(i,f,t,k);
        }

        sort(qs,qs+q,cmp);
        memset(res,-1,sizeof(res));
        for(int i = 0; i < q; i++){
            Q[qs[i].t].push_back(qs[i]);
            if(qs[i].f != qs[i+1].f){
                memset(low,0,sizeof(low));
                memset(dfn,0,sizeof(dfn));
                memset(instk,0,sizeof(instk));
                tot = 0, top = 0;
                dfs(qs[i].f,1);
                for(int j = 1; j <= n; j++){
                    Q[j].clear();
                }
            }
        }
        for(int i = 0; i < q; i++){
            printf("%d\n",res[i]);
        }
    }
    return 0;
}