We present a fast one dimensional hybrid method to efficiently simulate extensive air showers up to the highest observed energies. Based on precalculated pion showers and a bootstrap technique, our method predicts the average shower profile, the number of muons at detector level above several energy thresholds as well as the fluctuations of the electromagnetic and hadronic components of the shower. We study the main characteristics of proton-induced air showers up to ultra-high energy, comparing the predictions of three different hadronic interaction models: SIBYLL 1.7, SIBYLL 2.1 and QGSjet98. The influence of the hadronic interaction models on the shower evolution, in particular the elongation rate, is discussed and the applicability of analytical approximations is investigated.

17 pages (Latex), 15 postscript figures, 7 tables