A lattice model for dense polymer solutions and polymer mixtures in three dimensions is presented, aiming to develop a model suitable for efficient computer simulation on vector processors, with a qualitatively realistic local dynamics. It is shown that the bond fluctuation algorithm for a suitable set of allowed bond vectors has the property that due to the excluded volume constraint no crossing of bonds by local motions can occur, and entanglement restrictions thus are fully taken into account. For athermal binary (AB) symmetrical polymer mixtures, the dependence of both self-diffusion coefficient and interdiffusion coefficient on polymer density is obtained, simulating a thin film geometry where a film of polymer A is coated with a film of polymer B. For one density, the dependence of the interdiffusion coefficient on an attractive energy between unlike monomers is also studied. For weak attraction an enhancement of interdiffusion proportional to this energy occurs. For strong attraction, however, a rather immobile tightly bound AB layer forms in the interface which hampers further unmixing.