Two-loop anomalous dimension matrix for ΔS = 1 weak non-leptonic decays (I). O(α2s)