We consider the melting of the vortex solid in highly anisotropic layered superconductors with a small concentration of random columnar pinning centers. Using large-scale numerical minimization of a free-energy functional, we find that melting of the low-temperature, nearly crystalline vortex solid (Bragg glass) into a vortex liquid occurs in two steps as the temperature increases: the Bragg glass and liquid phases are separated by an intermediate Bose glass phase. A suitably defined local melting temperature exhibits spatial variation similar to that observed in experiments.PACS numbers: 74.25. Qt,74.72.Hs,,74.25.Ha,74.78.Bz The mixed phase of type-II superconductors with random pinning constitutes an excellent test system for studies of the effects of quenched disorder on the structure and melting of crystalline solids. In systems with weak random point pinning, the existence of a lowtemperature topologically ordered Bragg glass (BrG) phase with quasi-long-range translational order is now well established [1,2]. A variety of fascinating "glassy" behavior has been experimentally observed [3] near the first-order melting transition of the BrG phase in both conventional and high-T c superconductors. It has been suggested [3,4] that these observations can be understood if it is assumed that the melting of the BrG phase occurs in two steps: the BrG first transforms into a "multidomain" glassy phase which melts into the usual vortex liquid at a slightly higher temperature.In the presence of random columnar pinning, a "strong" Bose glass (BoG) phase [5] without quasi-longrange translational order occurs at low temperatures if the concentration of pins is larger than that of vortex lines. In the opposite limit of dilute pins, one expects [6] a "weak" BoG phase at low temperatures which would melt into an interstitial liquid (IL) as the temperature is increased. In the IL phase, some of the vortices remain pinned at the strong pinning centers, while the other, interstitial ones form a liquid. A recent numerical study [7] suggests that a topologically ordered BrG phase is also possible in such systems if the pin concentration is sufficiently small. It is also found experimentally, for both point [8] and columnar [9] pinning, that the melting of the solid phase is "broadened": the local transition temperature, measured by a discontinuity of the local magnetization, is different in different regions of the sample.Here we report results of a numerical study that provides insights and explanations for some of the observations described above. From minimization of an appropriate free energy functional, we find that the vortex system in an extremely anisotropic, layered, superconductor with a random dilute array of strong columnar pins (with both pins and magnetic field perpendicular to the layers) forms a BrG phase at low temperatures. As T is increased, this phase undergoes a first order transition into a glassy phase which we identify as a polycrystalline BoG. This phase then transforms, at a slightly higher T , into th...