md5 的偏移量是怎么确定的？

md5 规范给出了 4 轮偏移量来在 16 个 32 位块上执行循环左移。
每轮有 4 个偏移量的重复序列（ s in[abcd k s i]）：
第一轮是 7 、12 、17 和 22
第二轮为 5 、9 、14 和 20
第三轮为 4 、11 、16 和 23
第四轮为 6 、10 、15 和 21
rfc1321 上只说：
每轮的偏移量都经过了近似优化，以产生更快的“雪崩效应”。不同轮的偏移量是不同的。
他们是怎么找到这些偏移量的？什么保证了“雪崩效应”？