This paper presents fast algorithms for the computation of discrete Hartley transform (DHT). When the sequence length N = p*q, where p and q are integers and relatively prime, the one dimensional DHT can be decomposed into p length-q DHT's and q length-p discrete Fourier transforms (DFT). Compared to other reported algorithms, the proposed one has a regular computational structure, provides flexibility for composite sequence lengths and achieves substantial savings on the required number of operations.