I have a feeling that there may be a solution on the following lines, but this is far from a complete solution.
Let S be a subset of vertices. For each vertex V in the graph, consider the set D_S (V), which I define as follows: D_S (V) = {V} if V is in S, and otherwise D_S (V) is the union of {V} with the set D_S ( W) for all direct descendants of W from V. (That is, it is “all possible descendants of V, but stop the recursion every time you click on the element V”.) The question is, can we find a collection S such that size S is equal to O (f (N)), and also D_S (V) has size O (g (N)) for all V, and where f and g are asymptotically sublinear? (Maybe we can achieve sqrt for both, for example.)
If we can find this, I propose the following strategy. For preprocessing, create a hash table for each U in S containing all the vertices that are ultimately accessible from U. This can be achieved in O (f (N) * M); which is not necessarily better than O (N ^ 2), but better than O (M * Q), at least.
Now, to answer the query “is it accessible from V?”, It is trivial if V is in S. Otherwise, we check if V = U, and in this case it is also trivial. Finally, we apply the same process to all descendants of V, recursively, returning “yes” if the answer “yes” through one of two cases is higher, but “no” only if we never find U. This recursion takes O (g (N)).
The question remains how to choose S. I think that if a graph arises from some process where the external degree follows the power law, you can simply take the vertices sqrt (N) with the highest degree. But, for example, if we have a graph with N = 2 * K vertices (i, 0) and (i, 1), with K ^ 2 edges: from each (i, 0) to each (j, 1); then there is no suitable subset of S. But maybe graphs that do not have a suitable S need some degree of homogeneity, which we can use ... Or maybe not. I dont know. Any ideas let me know!