The first quantum technologies to solve computational problems that are beyond the capabilities of classical computers are likely to be devices that exploit characteristics inherent to a particular physical system, to tackle a bespoke problem suited to those characteristics. Evidence implies that the detection of ensembles of photons, which have propagated through a linear optical circuit, is equivalent to sampling from a probability distribution that is intractable to classical simulation. However, it is probable that the complexity of this type of sampling problem means that its solution is classically unverifiable within a feasible number of trials, and the task of establishing correct operation becomes one of gathering sufficiently convincing circumstantial evidence. Here, we develop scalable methods to experimentally establish correct operation for this class of sampling algorithm, which we implement with two different types of optical circuits for 3, 4, and 5 photons, on Hilbert spaces of up to 50, 000 dimensions. With only a small number of trials, we establish a confidence > 99% that we are not sampling from a uniform distribution or a classical distribution, and we demonstrate a unitary specific witness that functions robustly for small amounts of data. Like the algorithmic operations they endorse, our methods exploit the characteristics native to the quantum system in question. Here we observe and make an application of a "bosonic clouding" phenomenon, interesting in its own right, where photons are found in local groups of modes superposed across two locations. Our broad approach is likely to be practical for all architectures for quantum technologies where formal verification methods for quantum algorithms are either intractable or unknown.The construction of a universal quantum computer, capable of implementing any quantum computation or quantum simulation, is a major long term experimental objective. However, it is expected that non-universal quantum machines, that exploit characteristics of their own physical system to solve specific problems, will outperform classical computers in the near-term [1]. Ensembles of single photons in linear optical circuits are a recently proposed example: despite being non-interacting particles, their detection statistics are described by functions that are intractable to classical computers -matrix permanents [2]. It is therefore believed that linear optics constitutes a platform for the efficient sampling of probability distributions that cannot be simulated by classical computers, with strong evidence provided in the case of circuits described by large random matrices [3][4][5][6][7].A universal quantum computer, running for example Shor's factoring algorithm [8], creates an exponentially large probability distribution with individual peaks at highly regular intervals that facilitate the solution to the factoring problem allowing efficient classical verification, as is the case for all problems in the NP complexity class [9]. Accordingly, correct operation of the...