calculate lishu - phase1