calculate lishu - phase3