def get_DiAC_phases(cliffNum):
    """

    Returns the phases (in multiples of pi) of the three Z gates dressing the two X90
    pulses comprising the DiAC pulse correspoding to cliffNum
    e.g., get_DiAC_phases(1) returns a=0, b=1, c=1, in
    Ztheta(a) + X90 + Ztheta(b) + X90 + Ztheta(c) = Id
    """
    DiAC_table = [
    [0, 1, 1],
    [0.5, -0.5, 0.5],
    [0, 0, 0],
    [0.5, 0.5, 0.5],
    [0, -0.5, 1],
    [0, 0, 1],
    [0, 0.5, 1],
    [0, 1, -0.5],
    [0, 1, 0],
    [0, 1, 0.5],
    [0, 0, 0.5],
    [0, 0, -0.5],
    [1, -0.5, 1],
    [1, 0.5, 1],
    [0.5, -0.5, -0.5],
    [0.5, 0.5, -0.5],
    [0.5, -0.5, 1],
    [1, -0.5, -0.5],
    [0, 0.5, -0.5],
    [-0.5, -0.5, 1],
    [1, 0.5, -0.5],
    [0.5, 0.5, 1],
    [0, -0.5, -0.5],
    [-0.5, 0.5, 1]]
    return DiAC_table[cliffNum]