Abstract
Benzenoid chains, also known as hexagonal chains, are a class of organic compounds that consist of an arrangement of hexagonal rings fused together. In this article, we first present reduction formulas to compute the matching polynomial and independence
polynomial of any benzenoid chain by utilizing the transfer matrix technique. Subsequently, computational formulas for the Hosoya index and Merrifield-Simmons index of benzenoid chains are derived. Furthermore, the expected values of the
Hosoya index and Merrifield-Simmons index for random benzenoid chains are also obtained.