For planar integrable billiards, the eigenstates can be classified with respect to a quantity determined by the quantum numbers. Given the quantum numbers as m, n, the index which represents a class is c = m mod k n for a natural number, k. We show here that the entire tower of states can be generated from an initially given state by application of the operators introduced here. Thus, these operators play the same role for billiards as raising and lowering operators in angular momentum algebra.Quantum billiards are systems where a single particle is confined inside a boundary on which the eigenfunctions vanish [1]. One seeks the solutions of the time-independent Schrödinger equation, which is the same as the Helmholtz equation in the context of general wave phenomena. The solutions of this problem for an arbirarily shaped enclosure is a very challenging open problem, even when we restrict ourselves to two-dimensional cases [2,3,4,5]. There are some very interesting connections between exactly solvable models and random matrix theories, a summary may be seen in [6]. The Helmholtz operator is separable in certain coordinate systemsfor these cases, the solutions can be found [7]. The non-separable problems for which the classical dynamics is integrable have been recently studied in detail [8,9,10]. Although the solutions of these systems have been known, * Electronic address: srjain@barc.gov.in; Tel.: +912225593589 1 arXiv:1703.07587v1 [quant-ph]