A novel algorithm for computing the discrete Hartley transform (DHT) is presented. The proposed algorithm is based on the algebraic integer encoding scheme. With the aid of this scheme, an error-free representation of the cas function becomes possible. Furthermore, for the implementation of the algorithm, a fully pipelined systolic architecture with O(N) throughput is proposed.