编程实现计算贝祖等式的算法，即对于两个正整数a,b，可计算出整数s,t使得s·a+t·b=(a,b)

可以使用扩展欧几里得算法来计算贝祖等式。具体步骤如下：

1. 初始化s1=1, s2=0, t1=0, t2=1, r1=a, r2=b。

2. 计算商q和余数r3=r1 mod r2。

3. 如果r3=0，则(a,b)=r2，s=s2，t=t2，算法结束。

4. 否则，计算新的s和t：s=s1-q*s2，t=t1-q*t2。

5. 更新r1=r2，r2=r3，s1=s2，s2=s，t1=t2，t2=t。

6. 回到步骤2。

最终得到的s和t就是贝祖等式中的整数解。